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ABSTRACT 


The  purpose  of  thi^  study  is  to  examine  the  probieir.  of 

/ 

.determiniiig  rise  and  /set  tir.es  for  visibility'  periods  of 

/ 

multiple  satellite/-'satellite  and  sateliitevground  station 


communicat ions . 


Additionally,  the  ability  to  maintain 


constant  communication  links  is  explored. 

This  research  examines  three  means  to  solve  the 
problem:  basic  iteration,  the  Lavton  Method,  and  an  improved 
method  involving  parabolic  blending  and  a  closed  fc’^i  root 
solving  technique.  The  greatest  concentration  is  placed  on 
the  final  solution  and  its  application. 

The  result  of  this  investigation  is  a  computer  program 
;:i  ac  te:  mining  rise  anc  set  times  for  as  many  as 
four  •iatciiites  and  four  ground  stations.  In  addition,  the 
program  performs  a  check  for  continuous  communications 


betv.cen  the  indicated  satellites  and  earth  based  resources. 
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CHAPTER  1 


INTRC^'CTION 


Iii  space  mission  architecture,  much  rests  on  a 
satellite’s  ability  to  communicate  with  either  an  earth 
based  tracking  station  cr  another  satellite.  Data  relays  as 
well  as  tracking,  telemetry,  and  command  directions  depend 
OP  line  of  sight  opportunities.  Addit  ionall'',  availability 


oi 


communicat ions 


determine 


satellite 


data 


rate 


requirements,  data  storage  capabilities,  and  on  orbit 
processing.  All  these  factors  affect  satellite  complenity, 
orbit  selection  a.nd  satellite  constellation  size.  for  these 
rcasor.s,  it  is  impcrta.-it  to  possess  the  ability  to  rapidly 
oc :  ermine  ’.'.sioility  r.sv  nna  s:t  t-mci. 

Three  methoiis  of  determining  the  desired  rise  and  set 
timoo  arc  of  ii.tc-rcst  to  this  study.  The  first  propagaTes  a 
hoslc  •.  isibii.ty  functic.n  then  requires  linear  interpolation 
to  find  rise, ''set  times,.  The  second  takes  advantage  ;f  tiic 


fast  Fourier  transform 


define  characteristics 


of  the 


isibility  cur\e.  This  method  makc-s  use  of  the  Re  gala  Falsi 
itc-iati'.  e  techiiique  to  determine  rise,' set  times.  The  final 
method  reconstructs  the  visibility  cur'- c  using  parabolic 
blending,  then  determines  risc.'set  times  using  a  closed  form 
root  solver.  The  process  for  determining  visibility  periods 
by  this  third  method  is  det*»iled.  The  final  result  is  a 
computer  program  capable  of  handling  multiple  satellites  and 
ground  stations,  determining  whether  or  not  continuous 


ccmmunicav ions  among  the  satellites  and  with  the  ground 
possible,  and  completing  the  task  in  a  timely  manner. 


CHAPTLR 


'  VISIEUITY  FUNCTION 
^.1  General  Discussion 

The  concept  of  the  visibility  function  is  rei.tral  to 
any  mpthod  of  determinin?  rise/set  times  for  satcliito- 
satellite  or  satellite-ground  station  in-view  periods.  The 
visibility  function  involves  a  line  of  sight  calculation 
between  the  two  objects  of  interest.  The  value  cf  the 
function  varies  with  time  in  accordance  to  the  changing 
positions  of  the  two  objects. 

For  two  satellites  in  orbit  around  the  earth,  line  cf 
sight  is  poss.'.  ie  only  when  beth  objects  are  above  a  plane 
tangent  to  the  earth’s  surface.  The  '  most  e:;tr-.ne  r-csitions 
for  the ' sate  11 ites  which  still  allow  for  line  of  sight  occur 
when  both  lie  in  the  tangent  piano.  These  posit  iciiS  are 
defined  ty  the  fclicwir.g  cc^uot ic ns' : 

a*  *  cos 

g:  =  cos*  '  (  Rr  /  re  I 

Ihe  angle  between  thv  two  position  vectors  at  any  point  is 
defined  as; 

o  =  ■  !  n  •  r:  .) 

*r*i**T*"* 

•  i  i  I  I  A  *  I 

The  I isibiiity  function  is  then  defined  by'; 

,r  -  ai  a:  -  o  , 

When  the  value  thi  vi-ibility  function  is  positive,  the 
satellites  hav’e  line  of  sight.  Similarly,  a  negative  value 
indicates  no  line  of  sight.  A  value  of  zero  is  either  a 


rise  or  set  time  for  visibility  (See  Figure 
visibility  function  applies  to  a  satellite 
station. 


Orbital  Geometry 


-  OC,  *  0(:  -  Cp 


Figure  1 

2.2  Basic  Rise\Set  Determination' 

The  most  direct  method  of  determining  visibility 
opportunities  is  small  time  step  progression  of  the 
visibility  function.  At  each  step,  the  value  of  the 
visibility  function  xs  calculated.  When  plotted  versus 
time,  the  visibility  function  yields  a  visibility  waveform. 
When  a  sign  change  occurs  in  the  value  of  the  function  and 
the  wave  crosses  the  x  axis,  a  rise  or  set  time  is 
indicated.  Linear  interpolation  is  used  to  determine  the 
exact  time  of  the  crossing.  Although  this  method  is 
legitimate,  it  is  cumbersome  and  time  consuming.  The  small 
time  step  necessary  to  insure  the  best  solutions  using 


1  ) .  The  same 
and  a  .ground 


linear  interpolation  leads  to  large  data  sets  and  an 
extremely  inefficient  number  of  calculations  by  the 
computer.  In  order  to  ceter":ine  the  rise/set  time  for  two 
satellite  over  a  period  of  twenty-four  hours  with  a  time 
step  of  five  seconds,  more  than  seventeen- thousand 
calculations  of  the  visibility  function  alone  must  be 
computed.  As  a  method  for  real  time  solutions,  this  process 
is  simply  not  realistic. 


CHAPTER  3 


THE  LAWTON  METHOD  • 

3.1  The  Fast  Fourier  Transform 
In  his  paper,  J,  A.  Lawton  presents  a  method  for 
determining  visibility  opportunities  using  a  fast  Fourier 
transform,  FFT,  to  estimate  the  period  of  the  visibility 
waveform^ .  The  use  of  the  FFT  is  in  response  to  the  need  to 
process  large  quantities  of  data  points  with  fewer 

calculations  while  still  maintaining  the  accuracy  of  the 
interpolating  trigonometric  polynomial^  .  '  In  his 

presentation,  Lawton  uses  the  FFT  to ■  select  the  most 
prorinent  modulus  of  the  waveform  from  the  equation* 
a 

K  ^  1 

where 

ritkl  =  the  ’.isibility  function  at  time  tt 
i  =  >'-l 

n  =  number  of  data  points 

The  value  of  Cj  with  the  largest  modulus,  Cjina::,  is  then 
used  to  determine  the  value  of  the  period  of  the  waveform 
from  the ' following  equations* . 

n  ■  • 

P  =  2  n 


Once  the  period  of  the  vaveform  is  determine,  it  is 
used  in  conjunction  with  a  curve  searching  technique  to 
determine  the  rise/set  times  of  the  objects  in  question. 


3.2  The  Regula-Falsi  Method  of  Zero  Finding 


The  visibility  waveform  is  first  sectioned  by  fractions 
of  the  period  determined  by  the  FFT .  The  slope  of  the 
visibility  function  is  determined  in  each  progressive 
section  of  the  wave  until  a  local  maximum  value  of  the 
function  is  isolated  between  a  positive  and  negative  slope 
’.aluc-.  '  If  the  local  maximum  value  is  found  to  be  positive, 
.  n.:i  -  ^  at :  line  of  sight  is  possio-c.  ti.s  g-rt..;.  f 

tl.e  cur',  e  in  that  particular  seexion  is  searched  for  rise 
and  set  times.  If,  the  value  of  .the  isibility  function  is 
ncgati'.c  at  ti.o  local  maximum,  tnc  scarciiir.g  process  resumes 
to  isolate  the  next  local  maximum  •alue^ . 

Carve  searching  in  this  process  is  cccmplishcd  using 
t.ho  RegaJa  Falsi  method  of  zero  finding.  The  equations  is': 


Co  =  Uo  f !  bo  )  -  bo  f  {  ao  ) 
f (  bo  )  -  f  {  ao  ) 


An  illustration  of  the  process  is  shown  below  in  Figure  2. 


Figure  2 


If  the  value  of  f(co)  is  zero,  the  search  is  complete. 
Otherwise,  an  iterative  technique  begins  with  reassignment 
of  the  variables  occurring  as  follows^  : 

if  f(a«)f(^o)  <  0«  then  ai  »  ao  and  bi  s  co 
if  f(a«)f(co)  >  Oy  then  ai  s  co  and  bi  s  b^ 

The  above  reassignment  continues  and  the  value  of  c  is 
recalculated  until  f(co)  is  equal  to  zero.  In  the  initial 
process  of  the  curve  examination,  this  method  is  used  to 
determine  the  time  associated  with  a  local  maximum  in  order 
to  determine  if  the  value  of  the  visibility  function  at  the 
local  ma.ximum  is  positive  or  negative.  In  this  ease,  no 
crossing  of  the  x  axis  takes  place  and  the  iteration  of  the 
Regula  Falsi  method  continues  until  a«  «  b«  or  until  the 
difference  between  ao  and  b»  is  sufficiently  small. 


Although  the  Lawton  process  is  sound  in  the 
determination  of  rise/set  times,  it  has  the  disadvantage  of 
being  limited  to  orbits  with  eccentricities  less  than 
C.IO  Tiiis  is  necessary  to  ensure  a  relatively  periodic 

waveform.  Very  low  eccentricity  or  nearly  circular*  orbits 
produce  nearly  sinusoidal  waves  while  highly  elliptical 
orbits  produce  wave  forms  much  more  difficult  to  predict. 
Evidence  of  this  can  be  found  in  the  figures  included  with 
the  results  to  this  work.  Aperiodic  waveforms  may  result  in 
inaccurate  determination  of  the  period  and  make  curve 
searching  a  more  difficult  task.  Additionally,  this  method 
is  only  applicable  to  e\enly  spaced  data  points* ,  Sparse 
data  or  unevenly  spaced  points  ca.nnot  be  Used  with  the  fast 
F-.  V,  j  ic  i  t  ransl'v.  1  .-s  '. 

In  his  work.  Lawton  states  this  method  accounted  for 
19.9ii'7\  of  the  rise/set  times  determined  by  the  iteration 
method  described  previously  using  a  spherical-Eart  *-  riavity 
field  and  99.992%  of  the  rise/set  times  in  an  oblate-Earth 
gravity  field.'  He  fails  to  state  the  accuracy  of  the 
lise/set  times  found  by  this  method  as  compared  to  the  basic 
iteia^cn  techni<iue.  The  work  in  this  study  did  not  include 
reproduction  of  his  results  nor  investigation  of  the 
accuracy  of  this  method. 


CHAPTER  4 


THE  PARABOLIC  BLENDING  TECHNIQUE 
Parabolic  blending  is  a  method  of  curve  fitting  which 
maintains  continuity  of  two  parabolas  formed  from  four  data 
points.  The  first  parabola,  p(r)  is  defined  by  points  1,  2 
and  3,  while  the  second,  q(s)  includes  points  2,  3,  and  4. 
The  parabolically  blended  curve  is  formed  between  points  2 
and  3  (See  Figure  3)  and  is  defined  by^ 

C{t)  =  (1  -  t)p(r)  +  tq(s) 
where  0  <  t  <  1 


C(t)  -  (1  ~  t)p(r)  ♦  tq(s) 

Figure  3 


At  t  *  0,  C(t)  has  the  same  slope  as  p(r)  while  at  t  -  1, 
q(s)  ill  the  Influential  curve.  Thb  final  result  of  the 
blending  is  a  cubic  maintaining  characteristics  of  both 


parabolas . 

The 

equa 

tion 

of 

the 

curv 

e  i 

s  written 

C( 

t  )  = 

03  t 

3  + 

02  t^ 

+ 

Q1  t  +  oc 

where 

Uo 

= 

P2 

ai 

=  -  . 

5Pi 

+  *  L* 

P3 

02 

= 

Pi 

—  o 

5P2 

+ 

2  Pa 

-  .5P4 

aa 

-  - . 

5Pi 

+  1. 

5P2 

-  1 . 

5Pa 

+  ■ .  5P4 

This  technique  differs  from  the  cubic  spline  since  the 
equation  of  the  blended  curve  is  independent  of  all  other 
data  points  except  for  the  four  points  used  to  define  the 
two  parabolas^ .  Also,  parabolic  blending  works  well  with 
both  equally  and  unequally  space  data  points*  . 


CHAPTER  5 


FINAL  ALGORITHM 

The  final  result  of  this  study  is  a  computer  program 
designed  to  give  the  user  the  option  of  determining  the 
visibility  opportunities  between  a  combination  of  up  to  four 
satellites  and  four  ground  stations.  More  than  four 
satellites  or  ground  stations  requires  significantly  more 
calculations  and  data  storage.  The  step-by-step  process  of 
rise/set  determination  is  detailed  in  the  sections  which 
follow. 

5.1  .  Satellite  Position  Determination 

To  begin  the  process  of  determining  the  rise/set  times 
for  satellite-satellite  or  satellite-ground  station 
\isibility,  it  is  necessary  to  calculate  the  position  vector 
of  the  satellite  at  the  start  of  the  time  of  interest.  This 
requires  that  si?;  critical  orbital  elements  be.  known  for 
tliis  calculation.  They  are  the  mean  motion  (n),  inclination 
ti  j  ,  eccentricity  (e>,  longitude  of  the  ascending  node  ( D ) , 
argument,  of  periapsis  (w),  and  the  mean  anomaly  (M). 
Although  the  position  vector  can  be  determined  from  a 
variety  of  elements,  these  si.x  were  chosen  since  they  are 
readily  known  from  an  ephemeris  of  satellite  information. 
These  elements  and  the  position  of  the  satellite  are 
determined  with  respect  to  the  geocentric  reference  frame. 

The  first  step  in  the  determination  of  the  position 
vector  is  the  calculation  of  the  eccentric  anomaly.  This  is 


completed  through  the  use  of  the  Newton-Rhapson  Iteration. 
The  equations  are^  : 

Mn  =  En  -  e  SIN  (Ell) 

En  *1  -  En  +  '  'hi 

1  -  e  CO'S(En  ) 

For  the  first  iteration,  M  is  used  in  the  place  of  En . 
When  En t 3  is  calculated,  a  comparison  is  made  between  En  and 
En  ♦  1  .  When  the  difference  is  sufficiently  small,  the 
iteration  is  stopped.  For  the  final  computer  program,  the 
tolerance  was  set  for  a  difference  less  than  0.000001. 
Convergence  of  the  iteration  may  be  of  concern  for  orbits 
with  eccentricities  approaching  1.0  Once  determined,  the 

final  value  of  En  is  used  in  the  calculation  of  the  true 
anomaly  as  folloiv's^  : 

COS ( V  )  =  o  -  COS ( E ) 
e  COS(E)  -  1 

SIXIv)  =  as'hl_,-.e2.1.  •  SINJE) 

P 

P  =  a( 1  -  eCOS(E) ) 

Once  the  proper  quadrant  of  v  is  determined,  the  position 
vector  is  calculated  in  the  F'QW  or  perifocal  reference 
f  rame^  . 

r  =  rC0S(v»P  +  rSI.\Mv)Q 

r  =  a_U.r-e3-i 
1  +  eCOS(v) 

A  coordinate  transformation  is  then  accomplished  to  place 
the  position  vector  in  the  geocentric  framed *1®, 


The  above  process  determines  the  position  vector  of  the 
satellite  at  a  given  time.  ''ace  this  is  accomplished,  the 
position  of  the  satellite  must  be  propagated  through  the 
orbit  by  seme  time  step  dictated  by  the  user.  This  is 
accomplished  through  the  use  of  the  following  equations^ : 
n  =  no  { 1  +  1 . 5 J2 { V 1-e^  )  ( 1  -  1 . oSIN* (  i  )  )  ] 

Q  =  -( 1.5J2C0S( i 1 )n 
K  =  (1.5J2(2  -  2.5S,INM  i  )  ]  )n 

p2  • 

These  changes  in  the  values  of  the  longitude  of  the 
ascending  node  and  the  argument  of  periapsis  represent 
secular  perturbations  resulting  from  an  aspherical  earth.® 
Once  these  elements  are  propagated  forward,  the  position 
vector  is  recalculated.  In  order  to  decrease  compounding 
error,  the  position  vector  is  propagated  each  iteration  from 
the  original  values  of  G  and  w. 

It  is  important  to  note  other  perturbing  forces 
including  atmospheric  drag  were  not  taken  into  account  for 
this  study.,  Had  drag  been  included,  secular  changes  in  the 
serai-major  axis,  eccentricity  and  the  inclination  would  have 
been  taken  into  account.® 

5.2  Ground  Station  Position  Determination 

For  the  determination  of  satellite-ground  station 
visibility  periods,  the  position  of  the  ground  station  must 
be  calculated.  This  requires  the  user  know  the  following 
information:  the  year,  month,  day,  hour,  minute,  and  second 


the  simulation  is  to  begin,  the  east  longitude  and  geodetic 
latitude  of  the  site,  and  the  altitude  of  the  site  above 
mean  sea  level.  This  information  is  used  in  the 

determination  of  the  Julian  Date,  the  local  sidereal  time  of 
the  site,  and  the  final  calculation  of  the  position  vector. 

This  process  begins  with  the  determination  of  the 
Julian  Date  given  the  year,  month,  date,  and  time.  The 
Julian  Date  is  referenced  to  noon,  1  January  4713  B.G.  The 
equation  for  the  Julian  Date  used  in  this  program  only 
applies  to  dates  between  1  March  1900  and  28  February 
2100^’*°.  The  equation  is^ : 

7 'Year  +  int ( Month+9 ) ) 

JD  =  367*Year  -  int( _  _ 12  _  )  + 

\ 

!i Second  +  Mrn;  ! 

;  „  60. _  +  Hour! 

int ( 275*Month )  +  Day  +  1721013.5  *  1  GO  _  I 

9  24 

The  Julian  Date  is  then  used  to  calculate  the  Greenwich 
Sidereal  Time  at  0  hours  UT^ : 

9,0  =  100. ‘1606184  +  36000.  •  77004Tu  <•  '  . OO038793Tu 2 

where  Tu  =  Julian  Date  -  2451545.0 

36525.0 

The  Greenwich  Sidereal  Time  is  then  calculated  according  to 
the  rotation  of  the  earth  and  the  universal  time  in  terms  of 
angles  by2 ; 

0j  =  Gjo  +  w#  ( t  ~  to  ) 


1 


Finally,  the  Local  Sidereal  Time  for  the  site  is  determined 
by  in  terms  of  angles* : 

0  •=  eg  +  6e 

From  this  point,  the  position  vector  cT  the  site  in  the 
geocentric  frame  is  ,  calculated  from^ : 

X  =  : _ a.  +  H  :  COS(L) 

!  -  e.2SIN‘-.L)  : 

Z  =  1  a.  ( 1  -  e« M  +  H  1  SINML  ) 

;  -  e.^SINMLJ  ; 

r  =  X  C0S(9!I  X  JINtejJ  +  2  K 

In  a  manner  similar  to  the  satellite  position  vector 
propagation,  the  site  position'  v,;ctor  is  propagated  ahead 
v.'ith  respect  to  the-  rotation  of  ti.e  earth  arid  tiic  time  step 
indicated.  It  should  be  noted  that  both  the  satellite 
position  vector  and  the  ground  site  position  vector  may  be 
calculated  using  either  a  spherical  earth  or  an  oblate  earth 
in  accordance  to  the  users  choice.  If  an  oblate  earth  is 
chosen,  the  K  component  of  both  position  vectors  is  scaled 
by  a  factor  of  -  e«  ^  1 

5.3  Rise\Set  Time  Determination 

Once  the  position  vectors  of  the  tvo  objects  of 
interest  have  been  determined,  the  visibility  function  is 
calculated  at  each  time  step.  Although  any  time  step  mjay  be 
specified  by  the  program  user,  it  is  suggested  23  second 
steps  be  used  to  allow  for  computer  efficiency. ♦  Extremely 


large  time  steps  may  result  in  missed  rise/set  tiacs.  In 
order  to  iiisare  complete  evaluation  of  the  entire  time  of 
interest,  the  position  vectors  and  visibility  funetim  are 
calculated  for  times  one  step  earlier  and  one  step  later 
than  specified  by  the  user.  Starting  at  time  equal  to  zero, 
parabolic  blending  is  applied  to  the  first  four  visibility 
function  data  points.  The  equation  of  the  curve  bet«ce:i 
points  2  and  3  is  calculated.  The  real  roots  of  the 
resulting  cubic  equation  are  computed  using  a  closed  form 
root  solver* .  The  equations  to  compute  the  real  roots  of  a 
quadratic  are  included  in  the  program  in  the  event  the 
blended  curve  is  a  second  degree  polynomial  or  a  straight 
line.  All  compieM  roots  are  ignored  for  this  As 

stated  pro'.' i'-usly ,  the  blended  curve  is  dci^nec 

C( t  )  =  (1  -  t  )P( r )  +  tq{ s ) 

vhere  t  may  take  an'  value  between  0  and  1.  For  this 
reason,  only  roots  between  and  including  0  and  1  are 
considered. 

When  a  real  root  meeting  the  abov  ■■  requirements  is 
determined.  the  coefficient:;  ao  ,  qi  ,  a:,,  and  qj  are 
recalculated  using  the  times  corresponding  to  points  1,  2,  3 
and  1.  li  e  roots  of  the  original  blended  equation  are  then 
substituted  into  this  second  blended  equation.  The 

resulting  value  of  the  blended  equation  is  a  rise  or  set 
time.  The  root  is  then  increased  by  a  small  increment  and 
the  value  of  the  second  equation  is  recalculated.  If  the 
value  is  positive,  then  the  root  in  question  represents  a 


T'ise  time. 


Similarly,  a  negative  value  represents  a  set 

time. 

The  next  blended  curve  is  calculated  between  points  3 
and  I  and  then  checked  for  roots.  This,  process  continues 
until  the  entire  visibility  waveform  has  been  evaluated. 

The  result  is  a  rise/set  history  for  two  objects  over  a 
given  time  period. 

5.4  Continuous  Visibility  Check 
As  previously  mentioned,  the  final  computer  program 
offers  the  user  the  option  of  selecting  up  to  fodr 


satellites  and 

four 

ground 

stations 

for  evaluation  of 

rise/set  times. 

It 

is  the 

purpose 

of  the  continuous 

visibility  check  to  evaluate  the  rise/set  history  in  order 
to  determine  if  constant  communication  is  possible  among  the 
given  satellites  or  between  the  ground  and  a  constellation 
of  satellites. 

There  are  basically  two  typed  of  communication  links 
for  satellite  communications.  They  are  forward/reverse 
links  which  occur  between  satellites,  and  uplink/downi ink 
occurring  between  a  satellite  and  the  ground**.  For  two 
satellites,  forward/reverse  continuous  communication  is  not 
possible  if  a  rise  or  set  occurs  or  if  the  satellites  were 
not  in  constant  communication  from  the  beginning  of  the 
evaluation.  For  a  three  satellite  constellation,  at  least 
two  paths  of  visibility  out  of  a  possible  three  must  be 
maintained  at  all  times  for  constant  communication  to  be 


possible.  Similarly,  four  satellites  require  three  paths  of 


visibility  out  of  a  possible  six.  Additionally,  checks  are 
made  to  determine  if  each  ground  station  is  in  view  of  at 
least  one  satellite  at  all  times  for  the  purpose  of 
uplink/downlink  capabilities. 

The  continuous  visibility  check  is  accomplished  by 
constructing  a  table  containing  every  rise  and  set  time  for 
each  possible  communications  path.  For  example,  in  a  three 
satellite  constellation,  rise/set  times  are  determined 
between  satellites  si  and  ?2,  *1  and  “3,  and  s2  and  #3.  The 
table  for  this  constellation  would  included  all  the  rise  and 
set  times  from  these  satellites.  These  rise/set 

opportunities  represent  critical  times  since  they  are  the 
only  times  during  the  simuiaticri  when  the  status  of 
visibility  can  change.  In  order  to  accomplish  the 
c_nt:nuous  visibility  check,  the  status  of  each  visibility 
path  must  be  determined  at  every  rise  set  time.  For  the 
program,  a  (II  was  used  if  two  satellites  were  vis..ble  or 
experiencing  a  rise  at  a  given  time.  Similarly,  a  10)  was 
used  for  a  state  of  no  visibility  or  a  set.  An  example  of  a 
three  satellite  table  is  shown  below. 

Time  Sat’s  ?1  &  =2  Sat’s  =1  L  =3  Sat’s  =2  L  ~2 

i  sec  ) 

0  .  ,  1  0  1  . 
128  10  0 

356  0  0  1 

700  0  1  ,  .0 

825  0  1  0 

967  1  10 


As  mentioned  before,  for  constant  communication  to  be 
maintained  for  a  three  satellite  constellation,  at  least  two 
paths  of  visibility  must  be  available  at  all  times.  That 
means  the  sum  of  the  (l)’s  and  (0)’s  at  a  given  time  in  the 
table  must  be  greater  than  or  equal  to  two.  In  the  example 
above,  communication  is  broken  at  128  seconds  and  is  not  re¬ 
established  until  967  seconds. 

The  same  process  is  used  for  a  four  satellite 
constellation.  In  such  a  case,  there  are  six  possible 
communication  paths,  and  the  sum  at  any  given  rise  or  set 
time  must  be  greater  than  or  equal  to  three.  For  ground 
station  to  satellite  communication,  the  ground  station  must 
be  able  to  communicate  with  at  least  cne  satellite  at  any 
tin.e.  Therefore,  the  sum  at  each  rise  and  set  time  niMst  be 
greater  than  or  equal  to  one  for  constant  communication  to 


be  maintained. 


CHAPTER  6 


RESULTS 

Attached  are  the  results  of  three  different  runs  of  the 
final  computer  program  attached  in  Appendix  D.  Run 
represents  a  validation  of  the  rise/set  algorithm  of  the 
program.  The  results  of  this  program  were  compared  to  the 
results  contained  in  the  research  published  in  reference 
one.  The  initial  orbital  elements  of  the  four  satellites 
used  irt'the  program  are  found  preceding  the  rise/set  times. 
This  run  was  achieved  over  a  time  of  86,400  seconds  and  a 
250  second  time  step  with  a  spherical  earth.  These  results 
were  almost  identical  to  those  achieved  by  the  separate 
p’Ograir.  created  by  Aifano  and  Negron*,  Also  included  with 
the  first  run  results  are  the  rise/set  values  calculated 
using  the  basic  iteration  method  mentioned  earlier  in  the 
research.  For  this  basic  method,  a  10,000  second  time  span 
was  evaluated  in  5  second  intervals.  The  greatest  variance 
in  the  compared  results  was  less  than  1  second.  It  should 
be  noted  the  computer  time  to  generate  the  visibility  curve 
for  the  10,000  second  basic  method  between  two  satellites 
was  12  minutes.,  while  only  3  minutes  were  required  using  the 
parabolic  blending  method  over  a  86,100  seconds  with  a  230 
second  time  step.'  Run  «1  also  included  a  results  using  a 
ground  ststion  and  one  satellite.  This  run  involved  a  125 
second  time  step  over  86,400  seconds  and  an  oblate  earth. 
These  results  are  also  found  in  Appendix  A  and  compare 
almost  identically  to  those  achieved  by  Aifano  and  Negron* . 


The  Dasic  iteration  method  was  applied  to  this  data  as  well 
and  the  results  can  be  found  in  the  adjoining  colvimn.  The 
greatest  variation  between  the  parabolic  blending  method. and 
the  basic  iteration  method  in  this  case  was  approximately  3 
seconds. 

The  results  of  the  second  run  are  included  in  Appendix 
B.  These  i^esults  use  the  same  satellites  and  the  ground 
station  from  run  ?1  but  all  calculations  involved  the  use  of. 
an  oblate  earth.  As  expected,  the  oblate  earth  resulted  in 
earlier  rise  times  and  later  set  times^  .  Additionally.,  the 
continuous  visibility  check  was  employed  in  this  run.  Plots 
of  the  visibility  function  of  the  ground  station  \ersus  each 
of  the  four  satellites  is  included  in  this  Appendix. 
.;ddit  ior.aliy ,  the  combined  plot  of  all  four  satellite-ground 
station  functions  are  included  as  well  as  the  six  possible 
visibility  combinations  of  the  four  satellites.  These 
combination  plots  can  be  used  to  verify  the  results  of  the 
continuous  visibility  check. 

Appendix  C  contains  data  resulting  from  the  third  rxin 
of  the  program.  This  run  was  intended  to  check  the  complete 
capabilities  of  the  program  by  using  four  different 
satellites  and  four  different  ground  stations.  The  run  was 
accomplished  over  a  time  span  of  86,400  seconds  at  a  time 
step  of  125  seconds.  Data  for  the  satellites  was  chosen 
randomly  from  an  orbital  element  ephemeris  of  actual 
satellites  provided  by  the  United  States  Space  Command.  The 
four  ground  stations  chosen  are  Space  Shuttle  tracking 


sites.  They  include^:  (=1)  Kaena  Point,  Hawaii;  ls2)  Mahe 
Island,  Indian  Ocean;  (=3)  Thule,  Greenland;  { ?-!  )  Vandenburg 
Air  Force  Base,  California.  Again,  the  plots  of  the 
visibility  funct  ons  are  included  and  may  be  used  to  verify 
the  results  of  the  visibility  check.  The  total  run  time  for 
this  case  was  aporoximately  twenty  minutes. 


CHAPTER  7 


COKCLUSIONS 

While  the  basic  iterative  technique  and  Lawton’s  method 
are  viable  solutions  to  the  rise  and  set  time  problem,  the 
use  of  parabolic  blending  and  a  closed  form  solution  has 
many  advantages.  First,  parabolic  blending  is  a  flexible 
means  to  i  :construct  the  visibility  Curve.  It  can  be  used 
with  sparse  or  unequally  spaced  data  and  is  not  restricted 
to  near  circular  orbits.  Second,  a  closed  form  root  solver 
allows  rise/set  times  to  be  calculated  directly  rather  than 
by  an  iterative  technique  thereby  reducing  computer  run 
time.  The  resulting  timely  determination  of  rise/set  times 
makes  this  method  a  useful  tool  for  use  with  multiple 
satellites  and  ground  stations.  Additionally,  the  inclusion 
of  a  continuous  visibility  check  makes  this  program  a  useful 
tool  for  the  investigation  of  satellite  requirements  and 
space  mission  design. 
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Appendix  Ar  Program  Run  ?1 


Longitude  of  the  Ascending  Node  : 
Argument  of  Periapsis  : 


0.00000000  DEG 
O.OOOCOOOO  DEG 


INITIAL  ORBITAL  ELEMENTS  OF  SATELLITE  w2; 

Eccentricity  :  0.9363060000 

Inclination  :  64.98740000  DEG 

Mean  Anomoly  :  0.00000000  DEG 

Longitude  of  the  Ascending  Node  :  0.00000000  DEG 

Argument  of  Periapsis  :  0.00000000  DEG 


INITIAL  ORBITAL  ELEMENTS  OF  SATELLITE  =3: 

Eccentricity  :  0.0078742000 

Inclination  :  82.87090000  DEG 

Mean  Anomoly  :  0.00000000  DEG 

Longitude  of  the  Ascending  Node  :  0.00000000  DEG 

Argument  of  Periapsis  ;  0.00000000  DEG 


INITIAL  ORBITAL  ELEMENTS  OF  SATELLITE  ?4: 

Eccentricity  :  0.0048964000 

Inclination  :  144.64140000  DEG 

Mean  Anomoly  :  0.00000000  DEG 

Longitude  of  the  Ascending  Node  :  0.00000000  DEG 

Argument  of  Periapsis  :  0.00000000  DEG 


Rise/Set  Times, between  Satellite  and  Satellite  ?2: 


83679.13091 

SEC 

— 

SET 

Rise/Set  Times  between  Satellite 

41  and  Satellite 

1450.33207 

SEC 

- 

SET 

1450.222816 

SEC 

3965.41739 

SEC 

- 

RISE 

3965 .454545 

SEC 

6889.49771 

SEC 

- 

SET 

6889 . 599589 

SEC 

9371.13568 

SEC 

- 

RISE 

9371.062356 

SEC 

12399.31335 

SEC 

- 

SET 

1  1774.02357 

SEC 

- 

RISE 

18419.90807 

SEC 

- 

SET 

20218.05526 

SEC 

- 

RISE 

25546 . 76967 

SEC 

- 

SET 

27789.37858 

SEC 

- 

RISE 

30949.66914 

SEC 

- 

SET 

33360.02044 

SEC 

- 

RISE 

36352.15032 

SEC 

- 

SET 

38814.44584 

SEC 

- 

RISE 

41760.11003 

SEC 

- 

SET 

44235.69347 

SEC 

RISE 

47181.56457 

SEC- 

- 

SET 

49643.34461 

SEC 

- 

RISE 

52636.48892 

SEC 

- 

SET 

55045.12369 

SEC 

- 

RISE 

58209.92219 

SEC 

- 

SET 

6C445 . 55653 

SEC 

- 

RISE 

63807.11706 

SEC 

mm 

SET 

67619. 78663 

SEC 

- 

RISE 

71227.16594 

SEC 

- 

SET 

73604.36321 

SEC 

- 

RISE 

76628.56161 

SEC 

- 

SET 

79111.00284 

SEC 

- 

RISE 

22033.93906 

SEC 

- 

SET 

94549.30712 

SEC 

- 

RISE 

13412.02208 

SEC  - 

RISE 

17536.54531 

SEC  - 

SET 

19673.27556 

SEC  - 

RISE 

23787.20492 

SEC  - 

SET 

25923.63498 

SEC  - 

RISE 

30032.62420 

SEC  - 

SET 

32169.13711 

SEC  - 

RISE 

36275.30607 

SEC  - 

SET 

38411.96192 

SEC  - 

RISE 

42516.36090 

SEC  - 

SET 

44.653.14374 

SEC  - 

RISE 

48756.32667 

SEC  - 

SET 

50893.22948 

SEC  - 

RISE 

54995.50864 

SEC  - 

SET 

57132 . 55012 

SEC  - 

RISE 

61234.12410 

SEC  - 

SET 

63371.29703 

SEC  - 

RISE 

67472.30883 

SEC  - 

SET 

69609.58623 

SEC  - 

RISE 

73710.14286 

SEC  - 

SET 

75847.56839 

SEC  - 

RISE 

79947.69212 

SEC  - 

SET 

82085.26989 

SEC  - 

RISE 

86185.16424 

SEC  - 

SET 

Rise/Set  Times  between  Satellite 

t3  and  Satellite 

728.00260 

SEC  - 

SET 

728.253012  SEC 

2144.97084 

SEC  - 

RISE 

2145.755968  SEC 

3508.89011 

SEC  - 

SET 

3508.638444  SEC 

34217.02408 

SEC  - 

RISE 

35246.56198 

SEC  - 

SET 

, 

36889.15759 

SEC  - 

RISE 

38274.52867 

SEC  - 

SET 

39624.21934 

SEC  - 

RISE 

41174.78140 

SEC  - 

SET 

42839.99372 

SEC  - 

RISE 

43670.74866 

SEC  - 

SET 

72091.70409 

SEC  ~ 

RISE 

72474.06899 

SEC  - 

SET 

74570.33847 

SEC  - 

RISE 

75735 . 76905 

SEC  - 

SET 

77179. 52589 

SEC  - 

RISE 

78783.75035 

SEC  - 

SET 

80215.31302 

SEC  - 

RISE 

81418.26446 

SEC  - 

SET 

83415.21845 

SEC  - 

RISE 

83937.59649 

SEC  - 

SET 

THE  LOCATION  OF  GROUND  STATION  =1  . 
east  LONGITUDE:  0.000000 

GEODETIC  LATITUDE:  39.000000 

ALTITUDE  ABOVE  THE  EARTH  ELLIPSOID:  3.90OC0O 

initial  ORBITAL  ELEMENTS  OF  SATELLITE 
Eccentricity  :  0.0078742000 
Inclination  :  82.87090000  DET 

Mean  Anoooly  :  0.00000000  DEG 

Longitude  of  ti.e  Ascending  Node  :  0.00000000  DEC 

Argument  of  Periapsis  :  0.00000000  DEG 


RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  ?1  AND  SATELLITE  =1 

360.33816  SEC  -  RISE  -366.4128596  SEC 

793.74516  SEC  -  SET  793.3902878  SEC 

5973.86027  SEC  -  RISE  5972.4175820  SEC 

6041.56945  SEC  -  SET  6044.6470590  SEC 

39460.44479  SEC, -  RISE 
39949.95608  SEC  -  SET 

44907.37075  SEC  -  RISE 

452S2. 93512  SEC  -  SET 

86C61. 28564  SEC  -  RISE 


Appendix  B:  Program  Run  ?2 


Rise/Set  Ti»es  between  Satellite  s 
1160.52812  SEC  -  SET 
3954. 531U  .SEC  -  RISE 
6901.04371  SEC  -  SET 
9357.44428  EEC  -  RISE 
12414. 551C2  SEC  -  SET 
11750.81501  SEC  -  RISE 
18432.45195  3FC  •  ?7T 
20059.68746  SEC  '  RISE 
25578.36179  SEC  -  SET 
27770.72291  SEC  -  RISE 


Rise/Set  Ti«es  between  Satellite  ? 
1807.91575  SEC  -  SET 
4069.57257  SEC  ■  RISE 
7695.00886  SEC  -  SET 
9900.04560  SEC  -  RISE 
13549.39377  SEC  -  SET 
15662.26623  SEC  -  RISE 
19335.91349  SEC  -  SET 
21381.74668  SEC  -  RISE 
25066.47341  SEC  -  SET 
27132.76382  SEC  -  RISE 


Rise/Set  Tiaes  between  Satellite  = 
2997.97685  SEC  -  SET 
,  5924.02061  SEC  -  RISE 
8914.66759  SEC  -  SET 
11127.18527  SEC  -  RISE 
SEC  -  SET 

168C8. 98772  SEC  -  RISE 
19911.62120  SEC  -  SET 
22238.62982  SEC  -  RISE 
25373.12701  SEC  -  SET 
27696.99631  SEC  -  .RISE 


Rise/Set  Times  between  Satellite 
1932.12624  SEC  -  SET 
7104.92257  SEC  -  RISE 
11273.15055  SEC  -  SET 
13407.13572  SEC  -  RISE 
17537.31159  SEC  -  SET 
19669.08990  SEC  -  RISE 
23788.08220  SEC  -  SET 
25919.79045  SEC  -  RISE 


Rise/Set  Times  between  Satellite 
738.80937  SEC  -  SET 
2137.92179  SEC  -  RISE 
3516.84185  SEC  -  SET 


and  Satellite  *3: 


1  and  Satellite  ?4: 


2  and  Satellite' =3: 


2  and  Satellite  ?l: 


3  .and  Satellite  r4: 


fNOTE:  SO  ENTRIES  OR  MESSAGE  BELOW  INDICATES  CONTINUOUS 
Communications  broken  at  3516.84185  SEC 
Communications  re-established  at  3954.59411  SEC 
Communications  broken  at  8944.66759  SEC 
Communications  re-established  at  9357.44428  SEC 
Communications  broken  at  14479.19061  SEC 
Communications  re-established  at  14750.81501  SEC 
Communications  broken  at  19335.91349  SEC 
Tommunications  re-established  at  19669.08990  SEC 
Communications  broken  at  19941.62420  SEC 
Communications  re-established  at  20069.68746  SEC 
Communications  broken  at  25375.13704  SEC 
Communications  re-established  at  27132.76382  SEC 


COMMUNICATIONS. ) 


KiSK:  AND  SET  TIMES  BETWEEN,  GROUND  STATION  #1  AND  SATELLITE  ?!: 

RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  Pi  AND  SATELLITE  ?2: 

200.76093  SEC  -  RISE 
3060.0,9662  SEC  -  SET 
13147.6k  3  SEC  -  RISE 

RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  ?l  .AND  SATELLITE  =3: 
366.33816  SEC  -  RISE 
793.74516  SEC  -  SET 
5973.86027  SEC  -  RISE 
6041.56945  SEC  -  SET 
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RISE  .A.ND  SET  TIMES  BETWEEN  GROUND  STATION  si  .AND  S.ATELLITE  s4: 

CO.m’NIC.ATrONS  ST.ATUS  BETWEEN  GROUND  ST.ATION  si  .AND  THE  FOUR  GIVEN  S.ATELLITES: 
(NOTE:  NO  ENTRIES  OR  .MESS.AGE  BELOW  INDICATES  CONTINUOUS  COMMUNIC.ATIONS.  ) 


f 


Fii' 


TIME  FOR  BEGINNING  THE  SIMULATION:  17:30.45 

THE  LOCATION  OF  GROUND  STATION  #1 
EAST  LONGITUDE:  201.735000 
GEODETIC  UTITUDE:  21.571500 

\  ALTITUDE  ABOVE  THE  EARTH  ELLIPSOID:  0.000000 

THE  LOCATION  OF  GROUND  STATION  #2 

EAST  LONGITUDE:  55.480000 

GEODETIC  LATITUDE:  -4.669900 

ALTITUDE  .ABOVE  THE  EARTH  ELLIPSOID:  0.000000 

THE  LOCATION  OF  GROUND  ST.ATION  ?t3 

LAST  LONGITUDE:  291.400000 

GEODETIC  LATITUDE:  76.508200 

ALTITUDE  ABOVE  THE  EARTH  ELLIPSOID:  0.000000 

THE  LOCATION  OF  GROUND  STATION  #4 

LAST  LONGITUDE:  239.498300 

GEODETIC  UTITUDE:  34.823100 

ALTITUDE  .ABOVE  THE  EARTH  ELLIPSOID:  0.000000 

INITIAL  ORBITAL  ELEMENTS  OF  SATELLITE  #1: 

Eccentricity  :  0.0668128000 

Seai-BAjor  .Axis  :  117819.33914830  EM 

Inclination  :  37.35000000  DEG 

Mean  .Anoaoly  :  274.64810000  DEG 

Longitude  of  the  Ascending  Node  :  65.63070000  DEG 

Arguaent  of  Periapsis  :  79.11000000  DEG 

IN'ITI.AL  ORBITAL  ELEMENTS  OF  SATELLITE  *2: 

Eccentricity  :  0.0003109000 

Seai-aajor  Axis  :  42164.58832806  LM 

Inclination  :  0.00990000  DEG 

Mean  .Anoaoly  :  132.80340000  DEG 

Longitude  of  the  .Ascending  .Node  :  227.28640000  DEG 

Arguaent  of  Periapsis  :  359.92590000  DEG 

INITI.AL  ORBITAL  ELEMENTS  OF  S.ATELLITE  ?3: 

Eccentricity  ;  0.0145072000 

Seai-r ajor  Axis  :  7496.14095141  LM 

IncliL ition  :  90.26190000  DEG 

Mean  Aioaoly  :  246.05610000  DEG 

Longit'ide  of  the  Ascending  Node  :  107.70380000  DEG 

Arguaent  of  Periapsis  :  115.56590000  DEG 

INITIAL  Or.BITAL  ELEMENTS  OF  SATELLITE  *4: 

Eccentricity  :  0.0531098000 

Seai-aajor  A.xis  :  7210.27162153  KM 

Inclination  :  66.05630000  DEG 

Mean  A.iOBoly  245.56390000  DEG 

Longii  .'e  of  the  Ascending  Node  :  108.27480000  DEG 

Arguaent  of  Periapsis  :  119.97980000  DEG 


Rise/Set  Tiaes  between  Satellite  #1  and  Satellite  #3: 
2096.74276  SEC  -  SET 
3749.77870  SEC  -  RISE 
8657.22758  SEC  -  SET 
10396.52711  SEC  -  RISE 
15206.50054  SEC  -  SET 
17033.58950  SEC  -  RISE 
21747.79647  SEC  -  SET 
23658.24122  SEC  -  RISE 
28283.42317  SEC  -  SET 


Rise/Set  Tiaes  between  Satellite  #1  and  Satellite  ?4: 
1307.23847  SEC  -  SET 
3272.45823  SEC  -  RISE 
V 507. 23028  SEC  -  SET 
9472.56593  SEC  -  RISE 
13707.25247  SEC  -  SET 
15679.29431  SEC  -  RISE 
19906.64430  SEC  -  SET 
21891.90738  SEC  -  RISE 
26104.35521  SEC  -  SET 
28109.37531  SEC  -  RISE 


Rise/Set  Tiaes  between  Satellite  #2  and  Satellite  #3: 
9589.46045  SEC  -  SET 
9952.64793  SEC  -  RISE 
15090.00549  SEC  -  SET 
17190.80137  SEC  -  RISE 
21406.02460  SEC  -  SET 
23759.34470  SEC  -  RISE 
27853.51599  SEC  -  SET 


Rise/Set  Tiaes  between  Satellite,  52  and  Satellite  #4; 

146.81313  SEC  -  SET 
2048.36595  SEC  -  RISE 
7176.32283  SEC  -  SET 
9093.61752  SEC  -  RISE 
13570.46290  SEC  -  SET 
15790.10429  SEC  -  RISE 
19826.60084  SEC  -  SET 
22187.89728  SEC  -  RISE 
26069.57330  SEC  -  SET 
28469.81596  SEC  -  RISE 


Rise/Set  Tiaes  between  Satellite  #3  and  Satellite '#4: 
13706.79928  SEC  -  SET 
16641.32575  SEC  -  RISE 
18745.74240  SEC  -  SET 

COMMUNICATIONS  STATUS  BETVEEN  SATELLITE  #1,  #2,  #3,  .4N1) 
(NOTE:  NO  ENTRIES  08  MESSAGE  BELOW  INDICATES  CONTINUOUS 
Coaaunications  broken  at  15090.00549  SEC 
Coaaunications  re-established  at  15790.10429  SEC 
Coaaanications  broken  at  21406.02460  SEC 
Coaaunications  re-estabiished  at  22187.89728  SEC 
Conaunications  broken  at  27853.51599  SEC 
CoBBunications  re-established  at  28109.37531  SEC 
Coaaunications  broken  at  28283.42317  SEC 
CoMonications  re-established  at  28469.81596  SEC 

s 


#4: 

COMMUNICATIONS. 
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27189.14692  SEC  "  sil  . . 

RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  il  AND  SATELLITE  #2: 

RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  ?1  AND  SATELLITE  #3: 

6659.66079  SEC  -  RISE 
7205.64965  SEC  -  SET 
12742.74680  SEC  -  RISE 
13876.25363  SEC  -  SET 
19342.09950  SEC  -  RISE 
20205.80586  SEC  -  SET 

RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  #1  AND  SATELLITE  #4: 

12011.00325  SEC  -  RISE 
12888.92281  SEC  -  SET 
18219.91633  SEC  -  RISE 
19170.28556  SEC  -  SET 

COMMUNICATIONS  STATUS  BETWEEN  GROUND  STATION  AND  THE  FOLTl  GIVE.N  SATELLITES: 

(NOTE:  NO  ENTRIES  OR  MESSAGE  BELOW  INDICATES  CONTINUOUS  CO.MMUNICATIONS. ) 


RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  52  AND  SAitLLiie.  43 

18845.80416  SEC  -  RISE 

RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  #2  AND  SATELLITE  52: 

RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  f2  AND  SATELLITE  #3: 

2776.71326  SEC  -  RISE 
3809.60547  SEC  -  SET 
9232.43874  SEC  -  RISE 
10218.11934  SEC-  SET 

RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  =2  .\ND  SATELLITE  54: 

2575.55488  SEC  -  RISE 
3233.90819  SEC  -  SET 
8752.25542  SEC  -  RISE 
9560.57871  SEC  -  SET 

COMMUNICATIONS  STATUS  BETWEEN  GROUND  STATION  s2  AND  THE  FOUR  GIVEN  SATELLITES: 

(NOTE:  NO  ENTRIES  OR  MESSAGE  BELOW  INDICATES  CONTINUOUS  COMMUNICATIONS.) 

Communications  broken  at  0.000000000000000 

Communications  re-established  at  2575.55488485389 

Communications  broken  at  3809.60546671516 

Communications  re-established  at  8752.25541685105 

Communications  broken  at  10218.1193431246 

Communications  re-established  at  18845.8041619931 

k  ■ 


I 


M  •••*  *  W*1 


44 


RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  #3  AND  SATELLITE  ff2: 

RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  #3  AND  SATELLITE  ?3: 

863.15-109  SEC  -  RISE 
1929.46674  SEC  -  SET 
7439.15083  SEC  -  RISE 
8459.80178  SEC  -  SET 
14032.54055  SEC  -  RISE 
15028.79748  SEC  -  SET 
20602.42738  SEC  -  RISE 
21618.43278  SEC  -  SET 
27132.37091  SEC  -  RISE 
28188.78822  SEC  -  SET 

RISE  AND' SET  TIMES  BETWEEN  GROUND  STATION  #3  AND  SATELLITE  #4: 

994.31344  SEC  -  RISE 
1550.61962  SEC  -  SET 
7128.82589  SEC  -  RISE 
7759.49949  SEC  -  SET 
13314.48646  SEC  -  RISE 
13955.25665  SEC  -  SET 
19524.67234  SEC  -  RISE 
20130.86382  SEC  -  SET 
25748.10806  SEC  -  RISE 
26262.61775  SEC  -  SET 

COMMUNICATIONS  STATUS  BETWEEN  GROUND  STATION  *3  AND  THE  FOUR  GIVEN  SATELLITES: 

(NOTE:  NO  ENTRIES  OR  MESSAGE  BELOW  INDICATES  CONTINUOUS  COMMILNICATIONS. ) 

Communications  broken  at  0.000000000000000 

Communications  re-established  at  863.154088411520 

Communications  broken  at  1929.46674256308 

Communications  re-established  at  7128.82589091173 

Communications  broken  at  8459.80177798542 

Communications  re-established  at  13314.4864571971 

Communications  broken  at  13955.2566531969 

Communications  re-established  at  14032.5405513214 

Pommunicaticns  broken  at  15028.7974810107 

Communications  re-established  at  19524.6723377486 

Communications  broken  at  20130.8638220675 

Communications  re-established  at  20602.4273755799 

Communications  broken  at  21618.4327784578 

Communications  re-established  at.  25748.1080570079 

Communications  broken  at  26262.6177511727 

Communications  re-established  at  27132.3709121300 

Communications  broken  at  28188.7882164166 

♦ 


RISE  AND  SET  TIMES 

BETWEEN 

GROUND 

STATION’ #4 

AND 

SATELLITE 

#1: 

15864.12078 

SEC  - 

SET 

RISE  .AND  SET 

TI.MES, 

BETWEEN 

GROUND 

STATION  #4 

.AND 

S.ATELLITE 

t'li 

RISE  .AND  SET 

TIMES 

BETWEEN 

GROUND 

ST.ATION  fi 

AND 

SATELLITE 

#3: 

160.74120 

SEC  - 

RISE 

1145.77540 

SEC  - 

SET 

6518.52238 

SEC  - 

RISE 

7628.58607 

SEC  - 

SET 

13356.74253 

SEC  - 

RISE 

13863.55253 

SEC  - 

SET 

RISE  AMD  SET 

TIMES 

BETWEEN 

GROUND 

STATION  »4 

.AND 

S.ATELLITE 

54: 

6188.38578 

SEC  - 

RISE 

. 

7079.50299 

SEC  - 

SET 

12477.73425 

SEC  - 

RISE 

13319.74884 

SEC  - 

SET 

19231.12443 

SEC  - 

RISE, 

19321.51031 

SEC  - 

SET 

COMMUNICATIONS  STATUS  BETWEEN  GROUND  STATION  #4  .AND  THE  FOUR  GIVEN  SATELLITES: 
(NOTE:  NO  ENTRIES  OR  MESSAGE  BELOW  INDICATES  CONTINUOUS  COMMUNICATIONS.) 
Communications  broken  at  13864.1207815273 

(Tommunications  re-established  at  19231.1244274431 
Communications  broken  at  19321.3103139370 


45 


Or'd  Station  *1  and  Sort 


Appendix  D:  Program  Listing 


PROGRAM  TREE 


CO 

Q 

z 

< 

0 


m  I-  _  CM  "-t 

b  CO  fc  o  o 

H  2  2  o  o  o 

^  CC  b  CO  >  CO 

T  0  w  C5  O  0 


*  PURPOSE:  This  progreiB  determine':  rise  set  times  for  inviev  ♦ 

*  periods  for  satellite  to  satellite  or  satellite  to  * 

*  ground  station  viewing.  This  program  is  designed  * 

*  to  accommodate  up  to  four  satellites  emd  as  many  * 

*  ground  station.  The  following  information  is  * 

*  necessary  to  be  able  to  use  the  progran:  * 

*  Satellite  Orbit:  1)  Eccentricity  * 

*  -2)*Length  of  Semi  major  Axis  ♦ 

*  3)  Inclination  * 

*  4)  Mean  Anomoly  ♦ 

*  5)  Longitude  of  the  Ascending  * 

*  Sode  * 

*  6)  Argument  of  Periapsis  * 

*  Ground  Station:  1)  Julian  Date  in  Question  * 

*  2)  East  Longitude  * 

*  3)  Height  Above  Sea  Level  * 

*  Other  required  information  includes  the  amount  of  time  * 

*  of  interest  and  the  desired  time  step  for  increments-  * 

*  tion.  The  program  is  not  limited  to  any  particular  ♦ 

"  type  of  orbit.  -b  A»cAn  (vtortien  * 


VAEI.ABLES: 

CHOICE 

ELIP2 

Kl'M 

S2TU 

SIZE 

STEPS 

TIME 

TORN 

S080B 


the  number  of  satellites  indicated  by  the  user  * 
square  of  the  eccentricity  of  the  earth  * 
the  number  of  ground  stations  indicated  by  the  user  * 
converts  SEC  to  TU;  includes  the  steps  size  * 
step  size  (SEC)  * 
number  of  steps'  * 
total  time  of  interest  (MIS)  * 
yes  or  nc  response  ♦ 
spherical  or  oblate  earth  * 


aEAL‘8  ELIP2vTIME,SI2E,S2TU 
I STEGER  CHC I CE , SUM , SOROB . STEPS , YORK 

OPES;LSIT=25,FILt=’BISESET.DAT’ . STATl'S=  ’ SEK ' ) 

««**»«*«* A*«««y**»«****»*»«»«««*«*»a*W*«*****«**»»«******«»***«****«f 

•  USERS  CHOICE  *, 

y*»*ay«««a»««*a*»*«e*****««««****««*»«»«*»**«**«**««**«»«***»»**«»*«e 

VRITEl*,*)'  This  program  is  designed  ^o  determine  rise  and’ 
KF.ITE{ ’ set  times  for  satellite  to  satellite  viewing  as  well’ 

,  WRITE{*,*) ’as  satellite  to  ground  station  viewing.  If  you  choose’ 
VTilTEI •,• ) ’to  enter  ground  stations,  satellite  to  satellite  rise’ 
KRITE( *,* } ’and  set  times  will  be  computed  as  well  as  satellite’ 
VR1TE{*,* )’to  ground  station  rise  and  set  times.’ 

KRITE(*,*) 

1  fcRITE{*,*)’  Are  you  interested  in  satellite  to  ground  static  a’ 
k-BITE(*,*)’viewing?  DYes  2)No  ’ 

WRITE!*,*)’  »«=*>  ’ 

RE.AD(  *,*,  ERRS  1)  YORK. 


IF  (Y0RM.GT.2.0R.YPRS.LT.1)  THEN 

WRITE!*,*) ’Your  selection  was  not  a  valid  response.  Please’ 
WHITE!*,*) ’select  again.' 

GOTO  5 
ESDIF 

2  WRITE!* ,*) ’Are  yon  Interested  in  modelihg  a  spherical  earth’ 
WRITE!*,*) ’or  an  oblate  earth?  DSpherical  2)0blate' 
WRITE!*,*)’  •••**>  • 

READ(*,*,EER»2)S0R0B 


VaWA/U 

ELSE  IF  (S0R0B.EQ.2)  THEN 
EL1P2  =  0.O0669437999012D0 
E-SDIF 

♦  The  folIoKing  prsapts  the  user  to  enter  the  values  of  the  * 

*  tine  of  interest  and  the  step  size.  • 

■$**ti1f:****tt*********t***t***********************t******************** 

5  KKITE(*,*) ’W'hat  is  the  total  tiae  of  interest  (in  ainutes)  ==>  ’ 
READ!*,*,ERR=5)TIME 

VRITE( *,♦ ) 'The  tiae  entered  above  will  be  evaluated  in  steps.' 
VrRITE{*,«.' 'The  choice  of  step  size  is  up  to  you;  hovever,  it  ’ 
VRITE(*,*) 'should  be  noted  that  step  sizes  too  large  aay  result’ 
KEITE{*,*)’in  aissed  rise  and  set  tiaes.  A  suggested  step  size* 
VRITE(*,*)’is  250  seconds.’ 

6  VrRITE(*,*)’What  is  the  step  size  of  interest  (in  secs)  ==>  ' 
READ(*,*,ERR=6)SIZE 


TIME  =  TIME*60.0D0 
STEPS  =  DINT( TIME/SIZE) 

S2TI-  =  SIZE/806. 8103818D0 

t**tt*t**t**************t*******************************t************ 

*  The  following  proapts  the  user  to  enter  the  nuaber  of  ground  * 

♦  stations  and/or  satellites  he/she  desires  ♦ 

t***t:l^t*******9*********t******************************************** 


IF  (YORS.EO.l)  THEN 

'■*RITE(  ) 'Hok  many  ground  stations  do  you  desire. 
VRITEI*,*) 'make  your  selection  froa  below.’ 


WRITE! *,♦)’  1)  1 

WRITE! *,»)’  25  2 

WRITE! *,♦)’  3)  3 

WRITE!*, *1’  4)  4 

WRITE!*, •)’=*==>  ’ 
RE.'Jj(*,*,ERR=8)NUM 


Ground 

Ground 

Ground 

Ground 


Station’ 

Stations’ 

Stations’ 

Stations’ 


Please’ 


IF  iyUM.LT.l.OR.NL‘M.GT.4)  THEN 

WRITE!*,*) 'Your  entry  was  not  valid.  Please  reenter  the’ 
WRITE!*,* ) ’nuaber  of  ground  stations.’ 

GOTO  8 
ENDIF 


9  WRITE! *,*) 'How  aany  satellites  do  you  wish  to  enter?  Please’ 
WRITE!*,*) 'sake  your  selection  froa  below.’ 

WRITE!*, ♦)’  1)  1  Satellite  ’ 

WRITE!*,*)’  2)  2  Satellites’ 

WRITE!*,*)’  3)  3  Satellites’ 

WRITE!*,*)’  4)4  Satellites’ 

WRITE! ♦,*)’==«>  ' 

READ !*,*,ERRa9) CHOICE 

CALL  GA.NDS!CHOICE,ELIP2,NL’M,S2TU, SIZE, STEPS) 


ELSE 

10  WRITE!*,.*) ’Please  aake  a  selection  for  the  type  of  viewing’ 

WRITE!*,*) ’you  prefer;’ 

WRITE!*, •)’  1)  2  Satellite  Viewing’ 

WRITE!*,*)’  2)  3  satellite  Viewing’ 

WRITE!*,*)’  3)  4  Satellite  Viewing’ 

WRITE!*, •)’=«*>  ’ 

BEAD!  *,*,EBRwlO)CI10ICE 
NUM  *  0 


IF  (CHOICE.LT.l.OR.CHOICE.GT.S)  THE.*! 

WRITE!*,*) ’Your  choice  was  not  valid.  Please  re-enter' 


CALL  SATS I  CHOI CE , EL 1 P2 , S2TU , S I Z E , STEPS ) 


LNDIF 


CLOSE (UNIT=10) 
CLOSEiUKIT=ll) 
CLOSE ( UN I T= 12) 
CLOSE{UNIT=13) 
CLOSE(UNIT=14) 
CLeSEiUNIT=15) 
CLOSE(UNIT=25) 
CLOSE(UNIT=31) 
CLOSE { UNI T=32) 
CLOSE (UNIT=33) 
CLOSE(UNIT=34) 


END 


* 

* 

PURF-OSE: 

The  purpose  of  this  subroutine  is  to  coordinate  the 

* 

* 

deteraination  of  rise/set  times  between  satellites 

y 

* 

y 

and  ground  stations. 

y 

y 

* 

V.4RI.4BLES: 

y 

* 

INPUT: 

y 

* 

CHOICE  -  the  nuaber  of  satellites  indicated  by  the  user 

y 

* 

SL?! 

-  the  nuaber  of  ground  stations  indicated  by  the  user 

y 

y 

S2TI- 

-  converts  SEC  to  TU 

y 

y 

SIZE 

-  step  size  (SEC) 

* 

y 

y 

STEPS 

<■  nuaber  of  steps 

* 

y. 

y 

LOCAL: 

* 

y 

DAY 

-  1  to  31 

* 

y 

ELIP2 

-  square  of  the  eccentricity  of  the  earth 

* 

* 

ELONG 

-  east  longitude  position  of  the  ground  station 

* 

* 

GLAT 

-  geodetic  latitude  of  the  ground  station 

* 

* 

H 

-  altitude  above  the  earth  eliipsoxd 

t 

y 

HOUR 

-  0  to  23 

* 

t 

JD 

'  Julian  Date 

* 

* 

KIND 

-  deteraines  whether  continuous  visibility  involves 

* 

* 

a  ground  station  and  three  satellites  (0)  or  three 

y 

y 

satellites  (1) 

y 

y 

MIS 

-  0  to  59 

y 

y 

MON 

-  1  to  12 

y 

y 

SEC 

-  0  to  59 

y 

y 

STATE 

-  variable  that  defines  file  number 

y 

» 

T 

-  tine  corresponding  to  a  gi'-^n  step 

y 

* 

\AS 

-  variable  that  defines  file  number 

s 

y 

VIS 

-  value  of  the  visibility  function  at  a  given  step  , 

y 

y 

YEAS 

-  1900  to  2100 

y 

*ft««x****«*«3k*»««*«*«**«****«*****«*«*««y***«««'^»***«»*«*****«***»«*» 


SUBROITIXE  GAKDS{CHOICE,ELir2,XUM,S2Tl’, SIZE, STEPS) 

REAL*8  DAV , ELI P2 , ELOSG , GLAT , H , HOUR , JD , MI S ,MOS , SSTf 

EEAL*8  SEC, SIZE, T<0: 2005 >.V1S(0: 2005 l.VE^R 

INTEGER  CHOICE, COLNT,  IN,  I NTURN,KIND,.\TM. OUT, STATE, STEPS, VAR 

0PEN{UXIT=21,FILE=’GSP0S1.DAT’ ,STATUS=’NEV’ ) 

OPES ( USI T=32 .files’ GSP0S2 . DAT ’ , STATIS= ’ NEK ’ ) 

OPES  ( USI  T=2 3 ,  FI  LE=  ’  GSP0S3 .  DAT  ’ ,  STA7US=  ’  NEK'  ’ ) 
OPES{LSIT=34,FILE=’GSPOS4.DAT’ , STATUS^ ’ SEK ’ ) 

*y««*«*ai*«»«^«»'*«««»*««t****«***s**«*«********«*****4E*»**»*M****«*** 

«  la  order  to  accurately  deteruine  the  position  vectors  of  the  ♦ 

*  given  ground  station,  it  is  necessary  to  know  the  exact  tiae  .  * 

*  (Universal  Tiae)  the  siauiation  is  to  start.  The  following  * 

*  lines  proapt  the  user  to  enter  the  year,  aonth,  day,  hour,  ain-  * 

*  ate,  and  second  the  siauiation  is  to  start.  The  Julian  Date  * 

*  will  be  calculated  froa  this  Inforaation.  * 

*s«ssssssss«sssssssssss*ssssssssss*ss«*ssssssss*«**«s««ssss«ss*ss»e*s 

KllltE(*,*)’In  order  to  deterainc  position  vectors  of  the  earth’ 
KBITE(*,*) ’ground  stationls),  it  is  necessary  to  know  the  year,’ 
KRITE(*,4f)’Bonth,  day,  and  tiae  when  you  wish  the  siauiation  to  ’ 
KBITE(*,*) ’begin.  The  tiae  of  interest  is  I'niTersal  Tiae.’ 
KRITE(*,*) ’Please  enter  the  following  inforaation: ’ 
sssssssssssssssssssssssssssssssssssssssesssssassssssssssssssssessssss 

*  TEAS  * 

ssessnsssssssssasseesssssssssssssessssssssssssssssssssssssssssssssssss 

5  WBITEI*,*)’  What  is  the  year  of  interest,  between  1900  and  2100? 
WIITE(*.190) 

REaD(«.*.EBR*5)YEAR 


GOTO  5 
ESDIF 

********iti**t********t********t**************************t************ 

*  MONTH  ♦ 

*********t********************if***************t*t******************** 

10  WRITE!*,*)’  What  is  the  month  of  interest?  Please  enter  the  ’ 
WRITE!*,*)’  number  of  the  month  between  1  and  12.* 

WRITE! *,190) 

READ!*,*,ERR=10)MON , 

IF  !M0N.LT.1.0R.M0N.GT.12)  THEN 

WRITE!*,*)’  The  month  you  entered  was  not  a,  valid  response.’ 
WRITE!*,*)’  Please  re-enter  the  value.* 

GOTO  10 
ENDIF 

*if****************************************************************t** 

*  DAY  * 

*******li***ii*****tt*****%***t*%t*******************'*t**************** 

15  WRITE!*,*)’  What  is  the  day  of  interest?  Please  enter  a  value’ 
SUITE!*,*)’  between  1  and  31.’ 

STIITE!*,190) 

EEAD!*,*,ERR=15)DAY 
IF  !DAV.LT.1.0R.DAY.GT.31)  THEN 

SUITE!*,*)’  The  day  you  entered  was  not  valid.  Please’ 
STilTE!*,*)’  re-enter  the  value’ 

GOTO  15 

ELSE  IF  ! MON. EQ. 2. AND. DAY. GT. 29)  THEN 

SPRITE!*,*) ’February  does  not  have  more  than  29  days'.  Please’ 
WRITE!*,*) ’reenter  the  month.’ 

GOTO  15 

ELSE  IF  (DAY. EG. 31)  THEN 

IF  (MON.EQ.4.0E.MON.EQ.6.CR.MON.£Q.9,OE.MON.EO.H)  THEN- 
WRITE!*,*)’  The  month  you  entered  does  not  have  more  than  30 
S-RITE!*,*)’  days.  Please  re-enter  the  month.* 

GOTO  15 
ENDIF 
ENDIF 

«***4:4:***4i****«**4:**4>**:*««*«**4:***««**Y*y:tc*!t:i:««*x*«:fckic4i!|[4:**iic  *««*****« 

*  HOUR  * 

tit*i>c**t*****************************i********it***t*t****************** 

STJITE!*,* ) ’Next  it  is  necessary  to  enter  the  time  of  interest.’ 
WRITE!*,*) ’This  will  be  done  by  hour,  minute,  2ind  second.  ’ 
WRITE!*,*) ’Please  respond  at  the  prompt.’ 

WRITE!*,*) 

20  WRITE!*,*)’  What  is  the  hour  you  wish  the  simulation  to  begin?’ 
’'■RITE!*,*) ’  The  value  must  be  between  1  and  23.’ 

WRITE!*, 190) 

READ!*,*,ERR=20)HOUR 
IF  !HOUR.LT.O.OB.HOUR.GT.23)  THEN 

WRITE!*,*)’  The  hour  you  entered  was  not  valid.  Please  ’ 
WRITE!*,*)’  re-enter  the  hour. ’ 

GOTO  20 
ENDIF 

****************4*t***********ti************************************t* 

*  MINUTE  * 

.  e«*«**e*****««»*«***«*«*«*e«««*»««««*e»»«***«*«««««««i|iik*****«i«4i«*«**« 

25  WRITE!*,*)’  At  what  minute  in  the  hour, do  you  wish  to  begin.’ 
WRITE!*,*)’  Please  enter  a  number  between  1  and  39.’ 

WRITE!*, 190) 

READ! * , *, ERRs25 )MIN 
IP  (MIN.LT.O.OR.MIN.GT.59)  THEN  ' 

WRITE!*,*)’  The  minute  you  entered  was  not  a  valid  response.’ 
WRITE!*,*) *  Please  re-enter  the  value.’ 

GOTO  25 


^^^**^*m^m*w**********t********************************************** 

30  KRITE(*,*)’  At  what  second  do  you  wish  to  begin.  Please  enter’ 
K8ITE(*,*)’  a  nunber  betKecn  1  and  59.’ 

VE!TE(*,190) 

RE.«)(*,*,ERR=30)SEC 

IF  (SEC  LT.O.OR.SEC.GT.59)  THEN 

W8ITE(*,*)’  The  second  you  entered  was  not  a  valid  response.’ 
KRITEI*,*)’  Please  re-enter  the  value.* 

GOTO  30 
ENDIF 

HTII TE  ( 2  5 , 1 50 )  MON ,  DAY ,  YE.AR 
H'RITE(25,160)HOUR,MIN,SEC 
WRITEf25,*) 

***t**************t************************************************** 

*  The  following  line  calls  the  subroutine  JULDATE  where  the  Julian  ♦ 

*  Date  is  calculated.  * 

CALL  JULDATE ( DAY, HOLTl, MIN, MON, SEC, YLAR.JD) 

*  The  following  DO  loop  call  for  the  user  to  enter  the  east  ♦ 

*  longitude  and  the  geodetic  latitude.  The  ground  station  position  * 

*  vectors  are  then  determined  in  subroutine  GRSDST.  * 

****t***********:t***************************************************** 

DO  40  COUNT  =  l.NXM 

33  KEITE(*,*) ’Please  enter  the  east  longitude  value  for  the  ’ 
^»'RITE(*,200)COUNT 
WRITE! ♦,190) 

RE.\D(*,*,ERR=33)ELONG 

IF  (ELONC,LT.-3eO.ODO.OR.ELONC.GT.360.0DO)  THEN 

'•■RITE(*,*  1 'The  value  of  longitude  entered  was  not  valid.’ 
VRITE(*,*) ’Please  re-enter  the  value.’ 

GOTO  33 
ENDIF 

VBITE(25,220)C01NT 

WT?rTE{25,230)ELONG 

35  w^RiTEf*,*) 'Please  enter  the  e.odetic  latitude  for  the  ’ 

WRITE (*,200) COUNT 

KRITE(*,*) ’The  value  shoul'  be  between  -90  and  90  degrees.’ 
WTJITE(*,190) 

RE.AD(*,*,ER8=35)GLAt 

IF  (GLAT.LT.(-90.0DOJ.OR.GL4T.GT.{SO.ODO;)  THEN 

WRITE!*,* ) ’The  geodetic  latitude  value  you  entered  was  not’ 
VfRITE{*,*) ’acceptable.  Please  re-enter  the  value.’ 

GOTO  35 
ENDIF 

WRITE! 25, 240 )GLAT 

38  WRITE!*,*) ’Please  enter  the  altitude  of  the  ground  station  ’ 
WRITE!*,*) ’above  sea  level  (in  KM).’ 

WRITE!*, 190) 

READ!*,*,EES=38)H 
WRITE! 25, 250 )H 
WRITE!25,*) 

^  ^i********************************************************************** 

*  Next,  the  subroutine  GRNDST  is  call  where  the  position  vectors  of  * 

*  ground  station  dre  determined.  * 

**»*a*««a«««aaaa»a«**«a*a**«****«»«a**sa*««a*«*«*****a*»a*a**a*«**¥a*«* 

CALL  GENDST  ( COL’NT ,  ELI  P2 ,  ELONG ,  GUT ,  H ,  JD ,  S2TU ,  STEPS ) 

40  CONTINUE 

CL0SE(UNIT*31 ) 

CLOSE! UNIT=32) 


•*.  / 


itre  aubiuutiiie  oHis  IS  called  in  oruer  to  aeteraine  tne  rise  set  * 

*  times  between  the  satellites  entered  in  this  portion  of  the  * 

*  program.  * 

*^*t*t********i.**if*************************t:*************t***t******** 

'  CHOICE  =  CHOICE  -  1 

CALL  SATS { CHOICE , ELI P2 , S2TU .SIZE, STEPS ) 

CHOICE  =  CHOICE  +  1 

«**^4:4;4:4:»***«4:**4;*«:|::t:**4:****:|cv-*«*»**4::ti***t***4:4:«*«4:**«***4;*«****4:*«!>;»t 

*  The  following  opens  the  necessary  files,  then  records  the  vis-  * 

*  ibiiity  data  to  send  to  the  subroutine  BLEND.  * 

*t**ttt*tt**************-**************************Hi:t********t*******ii:* 

DO  70  OUT  =  1,NUM 

OPES(UNIT=60,FILE=’TRACKER.DAT* ,STATUS=’USKNOWK’ ) 
0PES(LSIT=41,FILE=’STAT1.DAT’ ,STATUS=’:;EK' ) 

OPEN ( UNI T= 4 2 , FI LE= ’ STAT2 . DAT ’ , STATUS= ' SEU ’ ) 

OPEN ( UNI T=4 3, FI LE=’ STAT3.DAT’ ,STATUS=’NEK’ ) 

OPEN  I  UN I T= 4  4 , FI LE= ’ ST AT  i . DAT ’ , STATUS= ’ NEV ’ ) 

V.AR  =  1 

DO  60  IN  =  1, CHOICE 

IF  (OUl.EQ.l.AND.IN.EQ.l)  THEN 
OPEN ( UNI T=3 1 , FI LE= ’ GSPOSl . DAT ’ , STATUS= ’ OLD ’ ) 

OPLN  ( UNIT= 1 , FI LE= ’ POSl . DAT ’ , STATUS= ‘ OLD ’ ) 

ELSE  IF  {0UT.EQ.1..4ND.IN.EQ.2)  THEN 
OPEN I UNIT=31 , FI LE= ‘ GSPOSl . DAT ’ , STATUS= ’ OLD ’ ) 

OPEN ( UNI T=2 , FI LE= ’ POS2 . DAT ’ , STATUS= ' OLD ’ ) 

ELSE  IF  {OUT.EQ.l.AND.IN.EQ.3>  THEN 
0PEN(UNIT=31,FILE=’GSP0S1.DAT’ ,STATUS=’OLD’ ) 

OPEN { UN I T= 3. FILES ’POS3.DAT’ ,STATUS=’CLD’ ) 

ELSE  IF  (OlT.EQ.l.A.SD.IN.EQ.4)  THEN 
0PEN(UNIT=31,FILE=’GSP0S1.DAT’ ,STATUS=’OLD’  ) 
0PEN(UNIT=4,FILEs’P0S4.DAT» , STATUS= ’ OLD ’ ) 


ELSE  IF  {0UT.E0.2..AND.IN.EQ.1)  THEN 
OPEN ( UN I T=32, FILES ’GSPOS2.DAT’ ,STATUS=’OLD’ ) 
OPEN ( UN I Ts 1 , FI LEs ’ POSl . DAT  * , STATUSs ’ OLD ’ ) 
ELSE  IF  (0UT.EQ.2,.AND.IN.EQ.2)  THEN 

OPEN ( UNI Ts32 , FI LEs ’ GSPOS2 . DAT ’ , STATUSs ’ OLD " ) 
OPEN(UNITs2,FILEs’POS2.DAT’ .STATUSs ’OLD’ ) 
ELSE  IF  (OUT.EQ.2..AND.IN.EQ.3)  THEN 
OPEN(UNITs32,FILE=’GSPOS2.DAT’ .STATUSs’OLD’ ) 
OPEN ( UNI Ts3, FI LEs ’ POS3 . DAT ’, STATUSs ’ OLD ’ ) 
ELSE  IF  (OUT.EQ.2.AND.IN.EQ,4)  THEN 
OPEN ( UNlTs32 , FI LEs ’ GSPOS2 . DAT ’ , STATUSs ’ OLD ’ > 
OPEN ( UNI Ts4 , FI LEs ’ POS4 . DAT STATUSs ’ OLD ’ )• 

ELSE  IF  <On.EQ.3.AND.IN.EQ.l)  THEN 
OPEN! UNI Ts33 , FILES ’ GSP0S3. DAT STATUSs ’ OLD ’ ) 
OPLNIUNITsl, FILEs’POSl. DAT’ .STATUSs’OLD’ ) 
ELSE  IF  (OUT.EQ.3.A.ND.IN.EQ.2)  THEN 
OPEN{UNITs33.FILE=’GSPOS3.DAT’ .STATUSs’OLD’ ) 
OPEN(UNITs2. FILEs’POS2. DAT’ .STATUSs’OLD’ ) 
ELSE  IF  {OUT.EQ.3.A.ND.IN.EQ.3)  THEN 
OPEN(UNITs33.FILEs’GSPOS3.DAT’ .STATUSs’OLD’ ) 
)  OPEN(UNITs3,FILEs’POS3.DAT’. STATUSs’OLD’) 

ELSE  IF  (OLT.EQ.3.ASD.IN.EQ.4)  THEN 
0PEN(UNITs33.FILEs’GSP0S3.DAT’ .STATUSs’OLD’ ) 
OPEN ( USI Tsi , FI LEs * POS4 . DAT’ , STATUSs » OLD ’ ) 

ELSE  IP  {OUT.E0.4.AND.IS.EQ.1)  THEN 
OPEN(UNITs34. FILEs’GSP0S4. DAT’ .STATUSs’OLD’ ) 
OPES{ UNITsl , files’ POSl .DAT* .STATUSs’OLD’ ) 


VbV  t 

ELSE  IF  (OUT.E().4.AND.IN.EQ.3)  THEN 
OPEN(CNIT=34,FILE-’GSPOS4.DAT’  .STATl'i^’OLD’ ) 

OPEN I  UNI T=3 , FI LE= ’ P0S3 . DAT ’ , STATU?= ' OLD ’ ) 

ELSE  IF  (OUT.EQ.4.AND.IN.EQ.4)  T;ii.N 
OPEN{USIT=34,FILE=’GSPOS4  OAf’ ,STATUS=’OLD’ ) 

OPEN  ( UNI T=4 ,  FI  LE=  ’  P0S4.  CAT  ’ ,  STATIS=  ’  OLD ’ ) 

ENDIF 

CALL  SATGRN ( I N , OUT , S I Z  E , STEPS ) 

CLOSE ( UN’ T=0LT+30) 

CLOSE{lNIT=IN) 

OPEN(UNIT=40,FILE=’GSVIS.DAT’ , STATUS= ' OLD ’ ) 

DO  50  INTl'RN  =  1,STEPS+1 

READ( 40, ♦IT! INTURN), VIS( INTURN) 

IF  (IN.EQ.l)  THEN 

OPEN{UNIT=10,FILE=’VIS1.DAT’  ,STATUS=’ UNKNOWN’ ) 
WRITE 1 10 , ♦ ) T ( INTURN ) , VI S ( I NTURN ) 

CLOSE! UNITslO) 

ELSE  IF  (IN.EQ.2)  THEN 

OPEN  ( UNI  T= 1 1 ,  FI  LE=  ’VIS2.DAT’.  STATUS=  ’  LTiKNOWK  ’  ) 
WHITE! 11 ,*)T! INTUBN) ,VIS! INTURS) 

CLOSE ! UNI T=ll) 

ELSE  IF  !IN.EQ.3)  THEN 

OPEN!UNIT=12,fILE=’VIS3.DAT’ , ST 4TUS=’ UNKNOWN’ ) 
fc-RI TE ! 12 , * ) T ! I NTURN ) . VI S ! INTURN ) 

CLOSE! UNI T=12) 

ELSE  IF  !IN.EQ.4)  THEN 

OPEN!UNIT=13,FILE=’VIS4.DAT’ , STATUS= ’ UNXNOWN ’ ) 
VR I  tE  !  1 3 ,  * )  T !  I  NTl'RN ) ,  V I S  ( I  NTURN ) 

CLOSE! UNIT=13) 

ENDIF 

50  CONTINUE 

H-RITE(25,*) 

WRITE! 25, 260 )OUT, IN 
CLOSE!UNIT  =  40) 

STATE  =  VAR+40  I 

CALL  BLEND ! STATE , STEPS , T , V] S ) 

CLOSE! UNI T=STATE) 

VAR=VAB+1 

60  CONTINUE 

IF  !CHOICE.EQ.2)  THEN 
WRITE! 25,*) 

WRITE! 25, 265 )OUT 
WBITE(25,285J 
STATE  =41 
CALL  GSSn'C2! STATE) 

ELSE  IF  !CHOICE.EQ.3)  THEN 
WRITE! 25,*) 

WRITE! 25, 270 )OUT 
WRITE! 25, 285) 

KIND  =0 
STATE  =  41 

CALL  CVC3!KIND,STATE) 

ELSE  IF  !CIK}ICE.EQ.4)  THEN 
■  HBITE!25,*) 

WRITE!  25, 280 )OUT 
WR1TE(25,285)  ' 

STATE  *  41 

___  CALL  G8SCVC4! STATE) 


70 


CONTINUE 


RETURN 

150  FORMAT! ’DATE  FOR  BEGINNING  THE  SIMUUTION:  ’ ,  12, ’/’ ,12, ,  11 ) 
160  FORMATCTIME  FOR  BEGINNING  THE  SIMULATION:  ’ ,  12,  ’ :  ’ ,  12,  ’ .  ’ ,  12 ) 
190  FORMAT (2X, ’====>  ’) 

200  FORMAT!’  location  of  ground  station  f’,Il,’.’) 

220  FORM.AT!’THE  LOCATION  OF  GROUND  STATION  ?’,I1) 

230  FOR>LAT!’ EAST  LONGITUDE:  ’,F12.6) 

240  FORMAT! ’GEODETIC  LATITUDE:  ’,F12.6) 

250  FORM.AT! ’ALTITUDE  ABOVE  THE  EARTH  ELLIPSOID:  ’,F12.6) 

260  FORMAT! ’RISE  AND  SET  TIMES  BETWEEN  GROUND  STATION  if’, II,’  .AND 
+  S.ATELLITE  ?f’,Il,’:’) 

265  FORMAT! ’COMMUNICATIONS  ST.ATUS  BETWEEN  GROUND  ST.ATION  #’,11,’  .A.\D 
+  THE  TWO  GIVEN  SATELLITES:’) 

270  FOR.MAT! ’COMMUNICATIONS  ST.ATUS  BETWEEN  GROUND  STATION  #’,I1,’  AND 
+  THE  THREE  GIVEN  S.ATELLITES:  ’ ) 

280  FORMAT! ’COMMUNICATIONS  ST.ATUS  BETWEEN  GROUND  ST.ATION  s’, II,’  .A.ND 
+  THE  FOUR  GIVEN  S.ATELLITES:  ’ ) 

285  F0R.MAT!’1N0TE:  NO  E-NTRIES  OR  MESSAGE  BELOW  INDICATES  CONTI. NUOUS 
+  COMMUNICATIONS.)’) 

END 


♦  PURPOSE:  Th?  purpose  of  this  subrcutinp  is  to  deteraine  the 

•  visibility  function  between  a  given  satellite  and  a 

♦  given  ground  station. 

* 


-  do  loop  counter  representing  the  nuaber  of  sat's 
"  do  loop  counter  representing  the  nuaber  of  ground 

stations 

-  step  size 

-  total  nuaber  of  steps 


-  do  loop  counter  * 

-  angle  between  the  position  vector  to  the  satellite  * 

and  the  position  vector  to  the  ground  station  * 

-  position  veetbr  to  the  ground  station  « 

-  position  vector  to  the  satellite  ♦ 

-  value  of  the  visibility  function  between  satellite  * 

and  ground  station  * 

Lawson,  32  * 

Alfano,  Negron,  and  Moore,  2  * 

**S*«««*««t««**«*«»*««*x»««**««»«**«***«*«»**4t*»*»*«**«***e «***«**«»»*» 

SUBROUTINE  SATGRN(IN,OUT. SIZE, STEPS) 

REALMS  PHI,R1(4),R2{4),SIZE,VIS 
I NTEGER  COUNT , I N , OUT , STEPS 

*  GSVIS.DAT  contains  visibility  data  between  the  given  satellite  and  * 

*  ground  station.  This  file  will  be  overwritten  each  tiae  SATGRN  is  * 

*  called.  * 

OPEN(UNIT=40,FILE=’GSVIS.DAT’ ,STATUS=’UXKNOKN’ ) 

DO  20  COUNT  =  0, STEPS 

RE.\D(0UT+30,*)R1(1),R1(2),R1{3),R1(1) 
RE.4D(IN,*)R2(1),R2(2),R2(3),R2{  1) 

»*«»4t*«***;»»****»»4;***««»«S«*S«««««**«4r***«««»***««**«*SS*«:t«**««*:»«»** 

*  The  next  lines  define  the  visibility  function  between  the  satellite* 

*  and  the  ground  station.  * 

f  Reference:  Lt  Col  Salvatore  Alfano  '  • 

*:*»*««*»**«*S««»»*«««**S*S«««S««S«»«***«****S*«*«*»«*»*S**«»*«*«*SSSS«« 

PHI  =  ((R1{1)*R2(1))+(R1(2)*R2(2))+{R1(3)*R2(3)))/(R1(4)*B2(1)) 

IF  (PHI.GT.d.ODO))  THEN 
PHI  =  l.ODO 

aSE  IF  {PHI.LT.(-l.ODO))  THEN 
PHI  s  -l.ODO 
'  ENDIF 

VIS  =  DAC0S({B1(4))/B2{4))  -  DACOS(PHI) 

)  KRITEI40,110)COUNT*SIZE,VIS 

20  CONTINUE 
CLOSE (UNITa40) 


*  V.ARI.ABLES: 

*  INPUT: 

*  IN 

*  OUT 

* 

*  SIZE 

*  STEPS 

* 

*  LOCAL: 

*  COUNT 

»  PHI 

* 

*  R1 

*  R2 

*  VIS 

* 

» 

*  REFERENCES: 

* 


RETUBN 


SUBROUTINE  JULDATE 


* 
» 

PURPOSE:  The  purpose  of  this  subroutine  is  to  convert  convention-* 

al  tiae  to  Julian  Date.  * 

* 


V.ARIABLES: 

* 

INPUT: 

* 

DAY 

-  1  to  31 

* 

HR 

-  0  to  23 

* 

MIS 

-  0  to  59 

* 

MON 

-  1  to  12 

* 

SEC 

-  0  to  59 

* 

YLAR 

-  1900  to  2100 

♦ 

OUTPUT: 

* 

JD 

-  Julian  Date 

* 

* 

Reference: 

Valado,  1 

* 

as******************************************************************** 

SUBROniNE  JULDATEIDAV, HOUR, MIS, MON, SEC, YEAR, JD) 

BEAL*8  DAY, HOUR, MIN, MON, SEC, T*EAB,JD 

JD  =  367 . ODO* ( YEAR ) -DI NT ( ( 7 .  ODO* { YEAR+DINT { ( MOS+9 . ODD ) /12 . ODO ) ) ) 

♦  /4.0DO)+DIST((273.0D0*MON)/9.OD0)fDAY+1721013.5D0+ 

+  (|((SEC/60.0D0)+MIN)/60.0D0)+HOUR)/24.0D0 


RETURN 

END 


*■ 

* 


SUBROUTINE  GRNDST 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

t 

* 

* 

* 

* 

* 

*■ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


PURPOSE:  The  purpose  of  this  subroutine  is  to  deternine  the 

geocentric  position  of  a  ground  station  given  Julian 
Date,  east  longitude,  and  geodetic  latitude. 

VARIABLES: 

INPUT:  , 

ELIP2  -  square  of  the  eccentricity  of  the  earth 

ELONG  -  east  longitude  of  the  station  of  interest 

GLAT  -  geodetic  longitude  of  the  station  of  interest 

H  -  altitude  of  the  station  above  the  earth  ellipsoid 

JD  -  Julifu)  Date 

S2TU  -  converts  SEC  to  TU 

STEPS  -  nuaber  of  steps 

LOCAL  V.4RIABLES: 

COSLAT  cosine  of  the  geodetic  latitude 

DEGRAD  -  conversion  fro«  degrees  to  radians 

DENOM  -  denoainator  of  the  equation  to  determine  G1  and  G2 

DT  -  change  in  time 

GLAT  -  geodetic  latitude  of  the  ground  station 
GSR(l-4)  -  position  vectors  of  the  ground  station 
KM2DU  -  conversion  froa  ka  to  DU 

REP  -  counter  for  do  loop 

SINL\T  -  sine  of  the  geodetic  latitude 

TGO  -  Greenwich  Sidereal  Tiae  at  0  hr  Universal  Time 
THETA  -  observers  local  sidereal  time 
THETDOT  -  time  rate  of  change  of  the  local  sidereal  time 
based  on  the, rotation  of  the  earth  (rad/TU) 

THETG  -  sidereal  time  of  the  Greenwich  Meridian 

TU  -  tiae  measured  in  centuries 

TVOPI  -  two  times  pi 

X  -  variable  for  the  determination  of  s^ticn  position 

vector 


« 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

« 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

t 

* 

* 


*  Reference:  Valado,  2-4  * 

*  Bates,  Mueller,  and  White,  98-103  * 

*  Escobal,  20-21  * 

*  Alfano,  Negron,  and  Moore,  6  .  * 


****t***t******t*.***1f:****'**********t*****1f********t*******'$*if*ittt*ii:**lft 

SUBROUTINE  GRNDST ( COUNT , ELI P2 , ELONG , GLAT , H , JD , S2TU , STEPS ) 

REALMS  COSLAT, DEGRAD, DENC»!,DT, ELI F2, ELONG, GLAT 
REALMS  GSR(4),H,JD,KM2DU,S2TU,SINLAT,TGO 
REALMS  THETA, THETDOT, THETG, TU,TV0PI,X 
INTEGER  COUNT, REP, STEPS 

t**********4;***4:*«***»«»«*»*4i*«««»«»«*«*4;****«**4:*:r»*«*»***«*s|[***«*«««it[ 

*  CONSTANTS  * 

*s»t*:»*4:*«:|:*«****4c*a*«:^e»****«*««»**««S**««««*«**«**«*««**«*«:(t****«**««* 

TWOPI  =  2.0DO*DACOS(-1.0DO) 

DEGRAD  =  TWOPI /360.0D0 
LM2DU  =  {1.0D0/6378.137D0) 

ELONG  =  ELONG*DEGRAD 
GLAT  =  GLAT*DEGRAD 
H  =  H»KM2DU 

THEn)6T\=  .058833590688786D0 

*  The  following  lines  calculate  the  Greenwich  Sidereal  Tiae  and  * 

*  .  to  be  used  in  the  deteraination  of  the  Local  Sidereal  Tiae.  * 

********t***1tt***************************9******n******$tn*tf********* 


TU  =  (DINT(JD)  ♦  .5D0  -  2451545. ODO)/36525.0DO 

IGO  »  1.753368559D0+628.3319705D0*Tl+ 

♦  .0C00067707O8127DO*(TU**2.0DO)+(6. 30038809866574)* 

♦  DELE (JD-DIKT(JD) -.500) 

TGO  =  DMOD(TGO,TWOPI) 

IF  (TGO.LT.(O.ODO))  THEN 
TGO  =  TWOPI+TGO 
EVDIF 


=  DSQRT(1.0D0-(ELIP2*{SINL4T*»2.0D0))) 

X  =  (|1.0DO/DENOM)+H)*COSLAT 

GSR(3)  =  (((1.0DO-ELIP2)/DENOM)+H)*SINLAT*(1.0D0/ 

+  DS0RT(1.0D0-ELIP2)) 

GSR(4)  =  DSQRT((X«2.0D0)+{GSR{3)*»2.0DO)) 

TGO  =  O.ODO 

*****t******tt****t******t*t******tt***t**********t*t********tt******* 

*  The  folloKing  lines  calculate  the  x  and  z  %Mriables  for  the 

*  position  of  the  ground  station.  These  rariables  are  then  used 

*  to  calculate  the  position  vector  in  the  UK  or  geocentric  coor- 

*  dinate  systea. 

********************************************************************** 

DO  50  REP  =  0, STEPS 
DT  =  REP*S2TU 
THETG  s  TGO  ♦  THETDOT*DT 
THETA  a  THETG  ♦  ELOKG 

GSR(l)  a  X*DCOS( THETA) 

GSEI2)  a  X*DSIN(THnA) 

WRITE<COI’NT+30,100)GSR(1),GSR(2),GSR{3).GSR14) 

50  CONTI SUE 
RETURN 

ICO  FORH.4T(F12.5,5.X,F12.5,5X,F12.5,5X,F12.5) 

END 


* 

SUBROUTINE  SATS 

* 

* 

* 

t 

a 

R’RPOSE: 

* 

VARIABLES: 

¥ 

* 

* 

ORBITAL  ELEMENTS: 

* 

* 

ARGP 

-  arguaent  of  periapsis  for  @  orbit 

(DEG) 

* 

♦  _ 

AXIS 

-  seai-aajor  axis  for  9  orbit 

(KM) 

* 

* 

ECC 

-  eccentricity  for  9  orbit 

* 

* 

INC 

-  inclination  for  9  orbit 

(DEG) 

* 

* 

MA 

-  aeM  anoBoly  for  9  orbit 

(DEG) 

* 

* 

a 

NODE 

-  longitude  of  the  ascending  node  for  @  orbit 

(DEG) 

* 

* 

OTHER  VARI.4BLES: 

* 

* 

AKM 

-  value  of  the  seai-aajor  axis  in  KM 

» 

* 

CHOICE 

-  reflects  the  users  choice  for  operation 

♦ 

* 

COUNT 

-  DO  loop  counter 

« 

t 

D2RA0 

-  converts  DEG  to  RAD 

* 

* 

DU2KM 

-  converts  DUs  to  KM 

« 

* 

ELIP2 

-  square  of  the  eccentricity  of  the  earth 

* 

KIND 

-  differentiates  between  a  ground  station  and  three 

* 

* 

satellites  (0)  and  three  satellites  (1) 

* 

* 

N 

-  aean  notion  of  the  satellites 

(DU/TU) 

* 

* 

S2TU 

-  converts  SEC  to  TU 

« 

* 

SIZE 

-  step  size  (SEC) 

* 

* 

STATE 

-  file  nuaber 

* 

* 

STEPS 

-  nuaber  of  steps 

* 

* 

T 

-  chronological  tine  over  tine  span 

* 

* 

TWOPI 

-  twice  the  value  of  pi 

* 

* 

VIS 

-  value  of  the  visibility  function 

* 

*  * 
*#y««»«:*4:***4:4!*««**4:*«*»*****:»**4:**4:*******4c*it:«*****4:  :)[*«****»**«****** 

SUBROUTINE  SATS (CHOICE, ELI P2,S2TU, SIZE, STEPS) 

REALMS  AK.M ,  ARGP,  AXI S ,  D2RAD ,  DU2KM ,  ECC ,  ELI  P2 , 1 NC ,  MA ,  NODE 
REAL*8  S2TU,SIZE,T(0:500),TVOPI,VIS(0:500),N 
INTEGER  CHOICE, COUNT, KIND, STATE, STEPS 

OPEN(UNIT=l,FILE=’POSl.DAT’ „STATUS=’NEK’ ) 
OPEN(UNIT=2,FILE=’POS2.DAT’ ,STATUS=’SEW’ ) 
OPEN(UNIT=3,FILE=’POS3.DAT’ ,STATUS=’NEW’ ) 
OPEN(UNIT=4,FILE=’POS4.DAT’ ,STATUS=’SEK* ) 

OPEN  ( UNI Ts60 ,  FI  LE=  ’  TR.4CKER .  DAT ! , STATUS=  ’  NEK ’ ) 
*«««:»*»4;»*4c:|c*,**«4:t********»***««*«***«*4:*4c4:******«*»********4:****4:*** 

*  CONSTANTS  * 

TWOPI  s  2 . ODO*DACOS { -1 . ODO ) 

D2R.AD  a  (THOPI/360.0D0) 

DU2KM  s  6378.137D0 

*««*«*«*«««*««*«««««««»«»«»«*«««*«««»««««*««*****«***«»««***«***«**** 

*  The  following  proapts  the  user  to  enter  the  values  of  the  * 

*■  orbital  eleaents  for  the  satellites  of  interest.  * 

*»«s«»*«s«s«*«*ss**a*«s****«**ss4t**s*#**««s*a************4c**a*****s** 

DO  20  COUNT  a  l,CHOICE+l 
ECC  a  O.ODO 
AXIS  a  O.ODO 
INC  a  O.ODO 
MA  a  O.ODO 
NODE  a  O.ODO 
ARGP  a  O.ODO 
N  a  O.ODO 

_ TO^ITEI  25 , 510 JCOUNT 


WRITE(25,520)ECC 

12  ViRITE{*i*)’Mean  Motion  (rev/day)  ==>  ’ 

READ{*,*,ERR=12)N 

N  =  N».05867301522d0 

AXIS  =  {1.0D0/(S**2.0D0))**(.333333D0) 

AKM  =  AXIS*DU2KM 
WRITE! 25, 530) ARM 

13  WRITE!*, *)’Inclination  ==>  ’ 

READ! *,*,ERR=13) INC 
W^ITE!25,540)INC 

14  WRITE!*,*)’Mean  Anonoly  !DEG)  ==>  ' 

READ!*,*,ERR=14)MA 

WRITE!  25, 550  )M.A 

15  WRITE!*,*) ’Longitude  of  the  Ascending  Node  !DEG)  ==>  * 
RL4D!*,*,ERR=15)NODE 

WRITE! 25, 560) NODE 

16  WRITE!*,*) 'Argument  of  Periapsis  !DEG)  ==>  ’ 
REA0!*,*,ERR=16)ARGP 

WRITE! 25, 570 )ARGP 
WRITEI25,*) 

WRITE!*,*) 

WRITE!*,*) ’Please  stemd-by.  Calculations  are  in  progress.’ 
WRITE!*,*) 

*  The  following  converts  the  values  entered  from  KM  to  DU,  or  from  * 

*  DEG  to  RAD.  * 

INC  =  INC*D2R.AD 
M.-\  =  MA*D2RAD 

NODE  =  N0DE*D2RAD 
ARGP  =  ARGP*D2R.4D 

CALL  POS !  COUNT ,  EL  I P2 ,  S2TU ,  STEPS ,  .4X I S ,  ECC ,  I NC ,  N ,  NODE ,  .ARGP ,  M.A ) 
CLOSE!  ITiilT  =  COUNT) 

20  CONTINUE 

IF  !CHOICE.EQ,0)  THEN 
GOTO  400 
ENDIF 

*S:t!4:*^S*SS«'S*S*S**S*S*S4:*****4:***S*««»****4:**«:^«*****4i**«*****»«»««#* 

*  In  the  following  if  loop,  the  computer  calls  the  subroutine  * 

*  necessary  to  check  all  options  of  visibility  depending  on  the  * 

*  number  of  satellites  in  consideration.  * 

***«SS**S*SSSSS*S*«****«*«**«*«**S*S*«**S««**«***««************SS«.*** 

*  2  SATELLITES  ,  * 

***SS4;S******S*y«***«**«S««*»***S**«**»**S*««»*«****»*«*S*««S««***S»* 

WRITE!*,*) ’The  computer  is  in  the  process  of  determining  •  ise’ 
WRITE!*,*) ’and  set  times.  Please  stand  by.’ 

IF  ICHOICE.EO.l)  THEN 
CALL  TWOVIS! STEPS, SIZE) 

OPEN!USITsl0,nLE=’VISl.DAT’ ,STATlSs’OLD’ ) 

DO  30  COUNT  s  1,STEPS+1 

READ! 10 ,* )T! COUNT ), VIS! COUNT ) 

30  CONTINUE 

CLOSE!UNITslO) 

WHITE!25,*> 

W8ITE!25,600) 

STATE  *  41 

OPEN(UNITs41,FILE=’STATl.DAT’ ,STATUS=’NEW’ ) 

_  CALL  BLEND!STATE, STEPS, T, VIS)  ■ 


ELSE  IF  ( CHOICE. EQ. 2)  THEN 
C.UL  THRVIS{ STEPS, SIZE) 

OPEN { UNI T= 10 , FI LE= ’ VI S 1 . DAT ’ , STATUS= ’ OLD ’ ) 

DO  40  COUNT  -  1,STEPS+1 

READ( 10 ,♦ )T( COUNT ), VI S { COUNT ) 

40  CONTINUE 

CLOSE(UNIT=10) 

WRITE(25,*) 

WRITE{25,600) 

STATE  =41 

OPEN(UNIT=41,FILE=’STATI.DAT’ ,STATLS=*NEW‘ ) 

CALL  BLEND (ST ATE, STEPS, T, VIS) 

CLOSE (UNIT=41) 

OPEIi  ( UNI T=ll ,  FI  LE=  ’  VI S2 .  DAT  ’ ,  STATUS=  ’ OLD ’  ) 

DO  50  COUNT  =  1,STEPS+1 

BEAD ( 1 1 , ♦ ) T ( COUNT ) , V I S ( COUNT ) 

50  CONTINUE 

CLOSE{UNIT=ll) 

WRITE(25,*) 

WRITE (25, 630) 

STATE=42 

OPEN(UNIT=42,FILE=’STAT2.DAT* ,STATIS=’NEV’ ) 

CALL  BLEND ( STATE , STEPS , T , VI S ) 

CLOSE (UNIT=42) 

OPEN(UNIT=12,FILE=’VIS3.DAT’ ,STA1US=’0LD’ ) 

DO  60  COUNT  =  1,STEPS+1 

RE.AD  ( 1 2 ,  ♦ )  T  ( COUNT ) ,  V I S  ( COUNT ) 

60  CONTINUE 

CLOSE(UNIT=12) 

WRITE(25,*) 

WRITE (25, 620) 

STATE=43 

OPEN ( UNIT=43 , FI LE= ’ STAT3 . DAT ’ ,STATIS= ’ NEW ’ ) 

CALL  BLEND( STATE, STEPS, T, VIS) 

CLOSE(UNIT=43) 

CLOSE(UNIT=60) 

KIND  =1  - 

STATE=41 

WR1TE(25,*) 

WRITE(25,660) 

WRITE(25,663) 

CALL  CVC3(KIND,STATE) 

'  *««*«««*«*««**»«*««*«»«*«««««««*«*«***«««*«***««**«***«***««««****«*« 

*  4  SATELLITES  * 

ELSE  IF  (CHOICE.EQ.3)  THEN 
CALL  FOUEVIS(STEPS,SIZE), 

OPEN(UNIT=10,FILEs’VISliDAT’ ,STATUS**OLD’ ) 

DO  70  COUNT  *  1,STEPS+1 

BEAD( 10, *)T( COUNT ) ,VIS( COUNT ) 

70  CONTINUE 

CLOSE (UNIT*10) 

WRITE{2S,*) 

WBITE(25,600) 

STATE*41 


0PEN(UHIT=11,FILE=’VIS2.DAT* , STATUS= ’ OLD ’ ) 
DO  80  COUNT  =  1,STEPS+1 

BEAD ( 1 1 , * ) T ( COUNT ) , V I S ( COUNT ) 

80  CONTINUE 

CLOSE(UNIT=ll) 

KRITE{25,*) 

KRITE(25,610) 

STATE=42 

OPEN ( UNI T=42 , FI LE= ’ STAT2 . DAT ’ , STATUS= ’ NEK ’  ) 
CALL  BLEND (STATE, STEPS, t, VIS) 

CLOSE(UNIT=42) 

OPEN{UNIT=12,FILE='VIS3.DAT' , STATUS= ’ OLD ’ ) 
DO  90  COUNT  =  1,STEPS+1 

BEAD( 12 ,*) T ( COUNT ), VI S ( COUNT ) 

90  CONTINUE 

CLOSE {UNIT= 12) 

KRITE{25,*) 

KBITE(25,620) 

STATE=43 

OPEN(UNIT=43,FILE=’STAT3.DAT’  ,'STATUS=  ’ NEW ‘ ) 
CALL  BLEND ( STATE , STEPS, T, VIS) 

CLOSE! UNIT=43) 

OPEN ( UNI Ts 13 . FI LE= ’ V I S4 . DAT ' , STATUS=  *  OLD ’ ) 
DO  100  COUNT  =  1,STEPS+1 

RLAD ( 1 3 , * ) T ( COUNT ) , V I S ( COUNT ) 

100  CONTINUE 

CLOSE! UNIT=13) 

KRITE!25,*) 

WRITE (25, 630) 


STATE=44 

OPEN ! UNIT=44 , FI LE= ’ STAT4 . DAT ’ , STATUS= ' NEK ' ) 
CALL  BLEND (STATE, STEPS,!, VIS) 

CLOSE(UNIT=44) 

OPEN ( UN I T= 1 4 , FI LE= ’ V I So . DAT ’ , STATUS= ’ OLD ’ ) 
DO  110  COUNT  =  1,STEPS+1 

BE.AD  ( 14 ,  * )  T  ( COUNT ) ,  V I S  ( COUNT ) 

110  CONTINUE 

CLOoE(UNIT=14) 

WRITE(25,*) 

K^ITE(25,640) 

STATE=45 

OPEN ( UNI T=4 5, files’ STATS. OAT* , STATUS= ’ NEK ’ ) 
CALL  BLEND(3TATE, STEPS,!, VIS) 

CLOSE(UNITs45) 

0PEN(UNITsl5,FILE=’VIS6.DAT’,STATUS=’OLD’ ) 
DO  120  COUNT  s  1,STEPS+1 

READ ( 1 5 ,♦) T ( COUNT ), V I S ( COUNT ) 

120  CONTINUE 

CLOSE(UNITslS) 

WSITE(25,«) 

KRITE(2S,650) 

STATEs46 

_ OPEN(yNITs46.FILEs*STAT6.DAT* .STATUSs’WEK’ ) 


STATE=41 

WRITE(25,») 

KRITE(25,670) 

WRITE (25, 665) 

CALL  CVC4( STATE) 

ENDIF 

CLOSE(UNIT=60,STATLS=’DELETE’ ) 

400  CONTINUE 
RETURN 

500  FORMAT!'  elements  of  satellite  #*,I1,*:’) 

510  F0RM.4T{’ INITIAL  ORBITAL  ELEMENTS  OF  SATELLITE  *’,11,’:’) 

520  FORMAT! 3X, ’Eccentricity  :  ’,F12.10) 

530  FORMAT !3X,’Seai-major  Axis  :  ’,F15.8,’  KM’) 

540  FORMAT! 3X, ’Inclination  :  ’,F15.8,’  DEG’) 

550  FORMAT! 3X, ’Mean  Anomoly  •  ’,F15.8,’  DEG’) 

560  FORMAT !3X, ’Longitude  of  the  Ascending  Node  :  ’,F15.8,'  DEG’) 

570  FORM.AT!3X, ’Argument  of  Periapsis  :  ’,F15.8,’  DEG’) 

600  FORMAT! ’Rise/Set  Times  between  Satellite  #1  and  Satellite  ?2:’) 

610  FORMAT! ’Rise/Set  Times  between  Satellite  #1  and  Satellite  #3:’) 

620  FORMAT! 'Rise/Set  Times,  between  Satellite  #1  and  Satellite  ?4;’) 

630  FORMAT! ’Rise/Set  Times  between  Satellite  #2  and  Satellite  #3:*) 

640  FORMAT! ’Rise/Set  Times  between  Satellite  >2  and  Satellite  54:’) 

650  F0RM.AT! ’Rise/Set  Times  between  Satellite  #3  and  Satellite  ?4:’) 

660  FORMAT! ’COm'NI  CAT  IONS  STATUS  BETV.’EEX  SATELLITE  =1,  =2,  AND  #3:’) 
665  FORMAT! ’(NOTE:  NO  ENTRIES  OR  MESSAGE  BELOW  INDICATES  CONTINUOUS 
+  CO'WUNl CATIONS.)’) 

670  FORMAT! ’COMMUNICATIONS  STATUS  BETWEEN  SATELLITE  #1,  ?2,  =3,  AND 
+  *4:’) 


SUBSOL'TINE  KJS 


PURPOSE:  The  purpose  of  this  subroutine  is  to  deteraine  the 
position  of  a  satellite  given  the  classical  orbital 
eleaents. 


* 

)  * 

'•-r  ♦ 

« 

♦ 


VARIABLES: 

INR'T: 

.  ARGP  -  arguaent  of  periapsis  for  §  orbit  (R.AD)  * 

AXIS  -  seai-aajor  axis  for  9  orbit  (DU)  ♦ 

COUNT  -  DO  loop  counter;  indicates  the  satellite  number  ♦ 

ECC  -  eccentricity  for  @  orbit  ♦ 

ELIP2  -  square  of  the  eccentricity  of  the  earth  * 

ISC  •'  inclination  for  @  orbit  (RAD)  * 

MA  -  aean  anoaoly  for  §  orbit  (RAD)  ♦ 

NODE  -  longitude  of  the  ascending  node  for  %  orbit  (RAD)  * 

S2TU  -  time  step  (TU)  * 

STEPS  -  nuaber  of  steps  in  the  requested  interval  ♦ 

♦ 

LOCAL  VARIABLES:  » 


COSARGP 

COSFISALE 

COSINC 

COSSODE 

COSNU 

DENOM 

DIF 

ECOSFE 


ENPLUSONE  - 


FINALE 


cosine  of  the  arguaent  of  periapsis  * 
final  value  of  the  cosine  of  the  eccentric  anoaoly  * 
cosine  of  the  inclination  * 
cosine  of  the  longitude  of  the  ascending  node  * 
cosine  of  the  true  anoaoly  * 
variable  for  the  denoainator  of  a  given  quotient  * 
difference  between  two  value  of  the  eccentric  * 
anoaoly  during  Sewton-Rhapson  Itteration  * 
eccentricity  aultiplied  by  the  cosine  of  the  ♦ 
finale  value  of  the  eccentric  anoaoly  * 
first  value  of  the  eccentric  anoaoly  for  use  in  * 
the  Sewton-Rhapson  Itteration  * 
second  value  of  the  eccentric  anoaoly  as  * 
calculated  in  the  Sewton-Rhapson  Itteration  * 
final  value  of  the  eccentric  anoaoly  reached  by  * 
convergence  of  the  Sewton-Rhapson  Itteration  * 


*  M  -  value  of  the  mean  anoaoly  at  each  time  step  * 

*  MN  -  value  of  the  mean  anoaoly  calculated  using  a  * 

*  given  value  for  the  eccentric  anoaoly  (used  in  the  ♦ 

*  Sewton-Rhapson  Itteration)  * 

*  MULT  -  matrix  used  to  rotate  the  position  vectors  from  * 

*  the  perifocal  reference  frame  to  the  geocentric  * 

*  reference  fraae  » 

*  N  -  aean  motion  of  the  satellite  * 

*  SANOM  -  anoaolistic  aean  motion  * 

*  KEH’ARGP  -  propogated  arguaent  of  periapsis  * 

*  NEKNODE  -  propogated  longitude  of  the  ascending  node  ^ 

*  P  -  calculated  value  of  the  seai-latus  recttia  for  the  • 

*  orbit  ♦ 

*  B  •  -  calculated  position  vectors  of  the  satellite  ♦ 

*  REP  -  DO  loop  counter  * 

*  SPQVf  -  position  vectors  of  the  satellite  in  the  perifocal  ♦ 

*  frame  ♦ 

*  SINARGP  -  sine  of  the  argxxaent  of  periapsis  * 

*  SINFINE  -  sine  of  the  final  value  of  the  eccentric  anoaoly  * 

*  SININC  -  sine  of  the  inclination  of  the  orbit  * 

*  SINNODE  -  sine  of  the  longitude  of  the  ascending  node  * 

*  SINNU  -  sine  of  the  true  anoaoly  * 

*  TWOPI  -  constant;  two  multiplied  by  pi  * 

*  ■  '  ♦ 

*  REFERENCES:  Bates,  Mueller,  and  White,  187,  219-222  * 

****a*a««a*«a****«***«*«*»»aa*«**«*s«*«****«*«*«*«******««**««***»«»*e 

_ SUBROLTINE  POS(C(XJNT,ELIP2,S2TU.STEPS,AXIS,ECC,INC,N,NODE,ARGP,MA) 


f  f  «'*f  4v»t«'44«  f  4^v.te  A  4%/t£tVl9f  iXCtiliUlUtf  f  r(£iH?(OUCt 

REAL*8  NODEDOT,  NODE,  SU,P,R(  4  ),RPQK(-1),S2TU,SINARGP 
REALMS  SINFISE,SISINC,SINNODE,SINNU,TWOPI 
INTEGER  COUNT, REP, STEPS 

*  The  following  computes  values  of  variables  Khich  will  be  used  * 

'  ♦  multiple  times  within  the  program.  This  is  done  to  decrease  the  * 

*  amount  of  calculation  necessary  within  the  subroutine.  * 

it******************************************************************** 

COSINX  =  DCOS(INC) 

SININC  =  DSIN(INC) 

J2  =  .00108262D0 

TWOPI  =  2.0D0*DACOS(-1.0D0) 

P  =  AXIS*n.0D0  -  (ECC»*2.0D0)) 

N  =  N*S2TU 

*4aa»**««««*««»'*4:*«**»a**»««**aa****«***4**«*4;**4:***4:«*t:)[*!):*;«***«4:«** 

*  The  following  calculates  the  values  necessary  to  propogate  the  * 

*  satellite  forward  in  its  orbit.  * 

4;4E4;***:|E«**4:4:*a.*«**a4:»a44***a****4!»«««a*«a*»»»«*»»«***«t***»«**:tt****«* 

SANOM  =  N*(1.0D0+((1.5D0*J2*{DSQRT(1.0D0-(ECC»*2.0D0)))/ 

+  ( P»*2 . ODD ) ) * ( 1 . ODO- ( 1 . 5D0* ( SININC**2 . ODD ) ) ) ) ) 

NODEDOT  =  -(1.5DO*(J2/P**2.0DO)*COSINC)*NANOM 
ARGPDOT  =  |1.5D0*{J2/P**2.0D0)*(2.0D0-(2.5D0* 

+  (SININC«2.0D0))))*NANOH 

DO  50  REP  s  0, STEPS 

NEWNODE  =  DMOD{(XODE  +  ( NODEDOT ) *REP )  ,TV(0P1 ) 

NEKARGP  =  DMOD((.ARGP  +  (.tRGPDOT)*REP)  ,TWOPI ) 

COSNODE  =  DCOS(NEh^N‘ODE) 

SINNODE  =  DSIN(NEKNODE) 

COSARCP  =  DCOS(NEWARGP) 

SINARGP  =  DSIXIXEK.4RGP) 

t*************1;1tt**1ft****t******************************^************ 

*  The  following  3x3  matrix  is  used  to  transform  the  position  vector* 

*  from  the  PQK  (perifocal)  frame  to  the  UK  (geocentric)  frame,  * 

MULT  (1,1)  =  C0SN'0DE*C0SARGP  -  SINXODE*SrNARGP*COSINC 

MULT(1,2)  =  -C0SN0DE*SINARGP  -  SINNODE*COSARGP*COSINX 

MULT (1,3)=  SI KN0DE*S I N I  NX 

MULT(2,1)  =  SIXNODE*COS.tRGP  +  COSNODE*SIK.4HGP*COSI>X 
MULT  (2, 2)  =  -SISNODF*SIN.ARGP  +  XOSXODE*COS.\RGP*COSINC 
>aLT(2,3)  =  -COSXODE*SININC 

MULT(3,1)  =  SIX.ARGP*SIXIXC 

MULT(3,2)  =  COS.tRGP*SININX 

MULT(3,3)  s  COSISC  . 

*  In  the  next  loop,  the  Xewton-Rhapson  iteration  is  carried  out  in  * 

*  order  to  determine  the  value  of  tne  eccentric  anomoly  at  every  * 

*  time  step  in  the  orbit  of  the  satellite.  * 

*«»«4:«»**»**«»4:««4:««*«***«*»y«««**«**«*«««*******»*»««»**«»*«**:t:.*«**»* 

M  =  DMOD((M.\  +  {XAN0M*REP)),TK0PI) 

EN  =  M 

10  MN  =  EN  -  ECC*(DSIS{EN)) 

ENPLUSONE  s  (EN  +  {  (M  -  MN)/(1.0D0  -  {EGC*{DCOS(EN) ) ) ) ) ) 

DIF  =  DABS (EN  -ENPLUSONE) 

IF  (DIF.GT.(.OOOOOIDO))  THEN 
EN  a  ENPLUSONE 

(X)TO  10  ;  ^ 

else; 

FINALE  =  ENPLUSONE  . 


ECOSFE  =  ECC»COSFINE 

DEKOM  =  AXIS*(1.0D0  -  (ECWOSFIKE) ) 

COSXU  =  (ECC  -  COSFISE)/! ECOSFE  -  l.ODO) 

SINKU  =  ((AXIS*{DSQRT( l.ODO  -  ECC*»2.0D0) ) )*SL\FINE)/DESOM 

NU  =  DATAN2(SINSL,COSNT) 

COSNT  =  DCOS(NU) 

SISNU  =  DSIN(NT) 

*tt*************'^f***********t*******************tt****ttt************* 

♦  In  the  next  three  lines,  the  position  vector  in  the  perifocal  * 

*  coordinate  systea  is  calculated.  * 

««4:«1i4!4:4:«*:^*t*«*»«*:tt****«**«*»*»«*a*4:**»*****'************************* 

RPQK ( 1 )  =  ( P*C0SKU ) / ( 1 . ODO  +  ( ECC*COSNT ) ) 

RPQK(2)  =(P»SISNU)/( l.ODO  +  (ECC»COSSU)) 

RPQK(3)  =  O.ODO 

♦  The  next  lines  transfora  the  RPQV  vector  to  the  geocentric  * 

*  coordinate  fraae.  ♦ 

R(l)  =  RPQW(1)*MLLT(1,1)  +  RPQW  ( 2 ) ''MULT  ( 1 , 2 ) 

R{2)  =  RPQW(1)*MLLT(2,1)  +  RPqK(2)»MULT(2,2) 

R(3)  =  (RPQW{1)*MULT(3,1)  +  RPQW(2)»MULT(3,2) )*(1.0D0/DS(3RT(1.0D0 
+  -  ELIP2)) 

R(4)  =  DSQRT(tR(l)«2.0D0)  +  (R(2)**2.0D0)  +  (R(3)*»2.0D0) ) 

KRITE(COUNT,100)R|1),R{2),R(3),R(4) 

50  CONTINUE 
RETURN 

100  FOR>LAT(F12.5,5X,F12.5,5X,F12.5,5X,F12.5) 


* 


END 


SLiu(UL"ii»e;  iwuvis 


« 

*  ,  * 

*  PURPOSE:  This  subroutine  exaaines  the  visibility  opportunities  ♦ 

*  between  two  satellites.  * 

*  * 

*  VARIABLES:  * 

*  INPUT:  ♦ 

*  STEPS  -  number  of  steps  in  the  requested  interval  * 

*  SIZE  -  step  size  (SEC)  * 

*  ♦ 

»  LOCAL  V.ARIABLES:  * 

*  COUNT  -  DO  loop  counter  * 

*  FILNL^  -  file  number  * 

*  PHI  -  tingle  between  position  vectors  R1  and  R2  ♦ 

*  R1  -  position  vector  to  satellite  #1  * 

*  R2  -  position  vector  to  satellite  ?2  * 

*  VIS  -  visibility  function  between  satellites  #1  and  f2  * 

*  * 

*  REFERENCES:  Larsor,-  32  ♦ 

♦  Alfano,  Negron,  and  Moore,  2  * 


*  (These  references  apply  to  THEVIS  and  FOLT.VIS  also)  * 

SUBROUTINE  TWOVIS( STEPS, SIZE) 

INTEGER  COUNT, STEPS 

REALMS  PHI,R1(4),R2(4),SIZE,VIS 

OPEN ( UNI T= 1 , FI LE= ’ POSl . DAT ’ , STATUS= ’ OLD ’ ) 

OPEN ( UNI Ts2, FI LE=’ P0S2.DAT’ ,STATUS=’OLD’ ) 

:(ft!|;***1!*»|!»**i(!***»*#*************»********»*******»****»************» 

*  VIS. DAT  CONTAINS  VISIBILITY  DATA  BETVEEN  SATELLITES  ONE  .AND  TTO  » 

ti'»c*t*******t********************************************************* 

OPEN(UNIT=10,FILEs’VISl.DAT’ ,STATUS=’NEK’ ) 

DO  20  COUNT  =  0, STEPS 

RE.AD(1,*)R1(1),R1(2),B1(3),R1(4) 

BEAD(2,*)R2(1),R2(2),B2(3),R2(4) 

*  The  next  lines  define  the  visibility  function.  .According  to  this  * 

*  idea,  consider  a  plane  tangent  to  the  earth  as  the  line  below  which* 

*  satellite  to  satellite  visibility  is  not  possible.  The  visibility  * 

*  function  calculated  here  compares  angles  to  determine  whether  the  * 

*  satellites  are  above  or  below  the  tangent  plane.  In  this  instance  * 

*  only  two  satellites  are  involved.  * 

4[t«4:*****4:**’t************!*************************************i******** 

*  Satellites  ?1  emd  =2  ■  .  * 

PHI  =  ((R1(1)*R2(1))+(B1{2)*R2(2))+(R1(3)*R2(3)))/(R1(4)*R2(4)) 

IF  (PHI.GT.(l.ODO))  THEN 
PHI  =  l.ODO 

ELSE  IF  {PHI.LT.(-l.ODO))  THEN- 
PHI  =  -l.ODO 
ENDIF 

VIS  =  (DACOS(1.0DO/B1(4)))+(DACOS{1.0DO/R2(4)))-(DACOS(PHI)) 
WRITE' 10, 110 )COUST*SIZE, VIS 

20  CONTINUE 

CLOSE(UNITsl) 

CL0SE(UNITs2) 

_ CLOSE(UNITglO) _ 


no  FOmT(I10,5X,F12.5) 


END 


SttJKULilNt  IMKMb 


*  .  * 

*  PURPOSE:  This  subroutine  examines  the  visibility  opportunities  * 

*  between  three  satellites.  * 

*  ■  ’  ■  ♦ 

*  V.ARIABLES;  * 

»  INPUT:  » 

*  STEPS  -  number  of  steps  in  the  requested  interval  ♦ 

*  SIZE  -  step  size  (SEC)  ♦ 

*  ■  .  ♦ 

»  LOCAL  V.ARIABLES:  * 

*  COUNT  -  DO  loop  counter  *  , 

*  PHI  -  angle  between  position  vectors  to  two  satellites  * 

*  R1  -  position  vector  to  satellite  #1  * 

*  ,R2  -  position  vector  to  satellite  =2  '  * 

*  R3  -  position  vector  to  satellite  #3  * 

*  VIS  -  visibility  function  between  two  satellites  * 


'  tic*t******************************t*****************  >■  *t.***t**t*ie****** 

SUBROUTINE  THRVIS( STEPS, SIZE) 

INTEGER  COUNT, STEPS 

REAL*8  PHI,R1(4),R2(4),R3I4),SI2E,VIS 

0PEN(UNIT=1,FILE=’P0S1.D.AT’  ,ST.ATUS=’OLD’ ) 

OPEN  ( UNIT=2 ,  FI  LE=  ’  P0S2 .  D.AT  ’ ,  STATUS=’OLD’ ) 

OPEN { UNIT=3 , FI LE=  ’  P0S3 . D.AT ’ , ST.ATUS=  ’OLD’ ) 

t**tt************t***t****************t**************t******t********  . 


*  The  following  table  lists  filenames  and  their  content  * 

*  UNIT  =  10;  VISl.D.AT  -  visibility  between  1  and  2  * 

*  UNIT  =  11;  VIS2.DAT  -  visibility  between  1  and  3  * 

*  UNIT  =  12;  VIS3.D.AT  -  visibility  between  2  and  3  * 


*ft*****t********t*************************************t**t*****t**** 

OPEN(UN1T=10,FILE=’VIS1.D.AT’ ,STATUS=’NEW’ ) 
0PEN(UNIT=11,FILE=’VIS2.D.AT’ , ST.ATUS=  ’ NEK ’  ) 
OPEN(UNIT=12,FILE=’VIS3.D.AT’  ,STATUS=’NEK’ ) 

DO  20  COUNT  =  0, STEPS 

RE.AD(1,*)R1(1),R1(2),R1(3),R1{4) 

RE.AD(2,*)R2(1),R2(2),R2(3),R2(  I) 
RE.AD(3,»)R3{1),R3(2),R3(3),R3{4) 

■ti^t********»***t*************************************-t*********^l:******* 

*  The  following  lines  calculate  the  visibility  function  as  explained  * 

*  in  the  subroutine  TVOVIS.  Since  three  satellites  are  involved  in  '* 

*  this  case,  visibility  checks  must  be  made  between  satellites  1  and  * 

*  2,2  and  3,  and  1  and  3.  * 

ta*************************************************************^******^^ 

*  Satellites  ?1  and  f2  ,  » 

y*»»«««X«4i*4:***4i«***4E*4i.;*#**«******4!4!***t***:ti»»*****»**'***********«1ci*:tc 

PHI  =  ((R1(1)*R2(1))+(R1(2)»R2(2))+(R1(3)*R2(3)))/(R1(4)*R2(4)) 

IF  (PHI. GT.ll.ODO))  THEN- 
PHI  =  l.ODO 

ELSE  IF  (PHI.LT.(-l.ODO))  THEN- 
PHI  -  -l.ODO 
E.NDIF 

VIS  =  {DACOS(1.0D0/Rl(4)))  +  ID.ACOS(1.0D0/R2(4) ) )  -  ID.ACOS(PHI)) 
KRI TE ( 10 , 1 10 ) C0UNT*SI ZE , VIS 

«*««*««««««*»»*4i««««4t«««*««««««******«*«**««««********»«*t*«*««****«**« 

*  Satellites  #1  and  #3  * 

«««»*««««*«*»«»««»*»»a«*«a«*»*«*aaaa«»a«**********iE*a»**»a«»»**««***a*« 

PHI  s  ((B1(1)*R3(1))+(B1(2)*R3(2))+(R1('3)*B3(3J))/{R1(4)*R314)) 


ELSE  IF  (?HI.LT.(-1.0D0))  THEN 
PHI  =  -l.ODO 
ENDIF 

VIS  =  (DACOS(1.CDO/R1{ 1)))  +  (DACOS{ 1.0D0/R3( J) ) )  -  (DACOS(PHI)) 
VRITE*  ll,110)COm*SIZE,VIS 

*  Satellites  72  and  t3  * 

»**♦***♦**$************,, 

PHI  =  ((R2(1)»R3(1))+(R2(2)*R3(2))+{R2(3)*R3(3)))/(R2(4)*R3(4)) 

IF  (PHI.GT.il.ODO))  THEN' 

PHI  =  l.ODO 

ELSE  IF  (PHI.LT.(-l.ODO))  THEN- 
PHI  =  -l.ODO 
ENDIF 

VIS  =  (DACOS(1.0DO/R2(4)))+(DACOS(1.0DO/R3(4)))-(DACOS(PHI)) 
H-RITE(12,110)COLNT»SIZE,VIS 

20  CONTI  NTE 

CLOSE (UNIT  =1) 

CLOSE! UNIT  =2) 

CLOSE{UNIT  =  3) 

CLOSEdNIT  =10) 

CLOSE (UNIT  =11) 

CLOSEdNIT  =  12) 

RETURN 

no  FOR>UT!I10,2X,Fi2.3) 

END 


SUBfiOLTINE  FOURVIS 


t 


* 

* 

*  PURPOSE:  This  subroutine  examines  the  visibility  opportunities 

*  between  four  satellites. 

* 

*  V.VRIABLES: 

*  INPUT; 

*  STEPS  -  number  of  steps  in  the  requested  interval 

*  SIZE  -  step  size  (SEC) 

* 

»  LOCAL  VARIABLES: 


« 

COUNT 

-  DO  loop  counter 

* 

PHI 

-*  angle  between  position  vectors  to  two  satellites 

♦ 

R1 

-  position  vector  to  satellite  »1 

R2 

-  position  vector  to  satellite  #2 

* 

R3 

-  position  vector  to  satellite  #3 

* 

B4 

-  position  vector  to  satellite  k4 

* 

VIS 

-  visibility  function  between  two  satellites 

SUBROUTINE  FOURVIS(STEPS,SIZE) 


INTEGER  COUNT, STEPS 

BEAL«8  PH1,R1(4),R2(4).R3(4),R4(4).SIZE,VIS 


OPEN ( UNI T= 1, FI LE=’ P0S1.DAT’ ,STATUS=’OLD* ) 
OPEN { UNI T=2 , FI LE= ’ POS2 . DAT ’ , STATUS=  *OLD ’ ) 
OPEN ( UN I T=3, FI LE=’ P0S3.DAT’ ,STATUS=*OLD’ ) 
OPEN ( t  NI T= 4 , FI LE= ’ P0S4 . DAT ’ , STATUS= ’ OLD ’ ) 


»  Th?  following  table  lists  filenaaes  and  their  content 


*  UNIT  =  10;  VISl.D/J 
»  UNIT  s  11;  VIS2.DAT 
»  UNIT  =  12;  VIS3.DAT 
»  UNIT  =  13;  VIS4.DAT 

*  UNIT  =11;  VIS5.DAT 
»  UNIT  =  15;  YIS6.DAT 


-  visibility  between  1  and  2 

-  visibility  between  1  and  3 

-  visibility  between  1  and  4 

-  visibility  between  2  and  3 

-  visibility  between  2  and  1 

-  visibility  between  3  and  4 


OPEN(UNIT=10,FILE=’VIS1.DAT’ ,STATUS=’NEK’ ) 
OPEN{UNIT=ll,FILE=’VIS2.DAT’ ,STATUS=’NEW’ ) 
OPEN ( UN I T= 1 2 , FI LE= ’ V I S3 . DAT ’ , STATUS= ’ NEW ’ ) 
OPEN(UNIT=13,FILE=’VIS4.DAT’ ,STATUB=’NEW’ ) 
OPEN(UNIT=14 ,FILE=’VIS5,DAT’ ,STArjS=’NEW’ ) 
OPEN(UNIT=13,FILE=’VIS6.DAT’ ,STATIS=’NEW ) 


DO  20  COUNT  =  0, STEPS 

RE.A0(1,»}R1(1),B1(2],B1(3),R1(4) 

READ(2,«)R2(1),B2(2},R2(3).R2(4) 

READ(3,*)R3(1),B3(2),B3(3),R3(4) 

READ(4,*)R4(l).R4(2),R4(3),R>'Ul 


yy**««»*«****«*e*««««««**«*««*4i*««*«w»*«*««*»«*«**«***«*««4a«e»««»«*e« 


•  ■* 


« 


The  following  lines  calculate  U.ie  visibility  function  as  explained  * 
in  the  subroutine  TWOVIS.  F  .iice  four  satellites  are  involved  in  * 


*  .  this  case,  visibility  che'ks  mist  be  aade  between  satell:.tes  1  and  * 

*  2,1  and  3,  1  and  4,  2  and  3,  2  and  4,  and  3  and  4. 

*  Satellites  #1  and  fZ 

««***««*»**«*«***««*«»****««»«**«**«»«k««*«*««************«:i««k«*k*k*«« 

PHI  =  ((Ria)*R2(l))>(Rl(2)*R2(2))t(Bl(3)*R2(3)))/(Rl  4)»B2(4)) 


IF  (PHI.GT.d.ODO))  THEN 
PHI  =  l.ODO 

ELSE  IF  (PHI.LT.(-l.ODO))  THEN 


•  Satellites  «l  and  #3  * 

PHI  =  ((81(1)*83(1))+{R1(2)*R3{2))+(R1(3)*R3{3)))/IR1(4)*R3(4)) 

IF  IPHI.GT.d.ODO))  THEN 
PHI  3  l.ODO 

ELSE  IF  (PHI. LT. (-l.ODO))  THEN 
PHI  »  -l.ODO 
ENDIF 

VIS  3  (DACOS(1.0D0/81(4)))  *  (DACOS(1.0DO/83(4) ) )  -  (DACOS(PHI)) 
WBITE(ll,nO)(»UNT*SIZE,VIS 

«««*a«**««*3*«*«****«aa«s«*sa*«*««*««**«**«**«««**«a««***s«»««***a****« 

*  Satellites  #1  and  #4  * 

*s*ssse«s*»***«««««*ee*e«**«*«««*»«**«**«**************«*«**«******s«»* 

PHI  3  ((81(1)*B4(1>)>(B1(2)*B4(2))4(S1(3)*84(3)))/(R1(4)*S4(4)) 

IF  (PHI. GT. (l.ODO))  THEN 
PHI  3  l.ODO 

ELSE  IF  (PHI. LT. (-l.ODO))  THEN 
PHI  3  -l.ODO 
ENDIF 

VIS  3  (DACOS(l.ODO/Bl(4)))  ♦  (DACOS(1.0DO/B4(4) ) )  -  (DACOS(PHI)) 
V|-BITE(  12, 110)COUNT*SI2E,VIS 

*«**»««*«««*«*«««**««*««**««*««««««««**«*»««*«*•««*«*******««******«•** 

•  Satellites  »2  and  e3  * 

PHI  3  (|R2(1)*83(1))+(R2(2)*R3(2))+(R2(3)*R3(3)))/(R2(4)*R3(4)) 

IF  ( PHI ;GT.( l.ODO))  THEN 
PHI  3  l.ODO 

ELSE  IF  (PHI. LT. (-l.ODO))  THEN 
PHI  3  -l.ODO 
ENDIF 

VIS  3  (PaCCS(1.0DO/B2(4)))*(DACOS(1.0D0/B3(  l)))-(DACOS(PHI)) 
V8ITE(13,110)COCNT*SIZE,VIS 

*  Satellites  *2  and  rl  * 

***«e«*e««s***«««*«*«»«*««««««««««*««««««*»«**«*»»«*«*«*3**«**«*****«*« 

PHI  «  ((B2(l)«B4(l))*(R2(2)«84(2))f(B2(3)«B4(3)))/(R2(4)»B4(4)) 

IF  (PHI.GT.i l.ODO))  THEN 
PHI  3  l.ODO 

ELSE  IF  (FHI.LT.(-l.OOO))  THEN- 
PHI  3  -l.ODO 
ENDIF 

VIS  3  (DACOS(1.0DO/R2(4)))^(DACOS(1.0DO/B4(4)))-(DACOS(PHI)) 

VBI TE  ( 1 4 , 1 10  )(NX.-ST*SI  ZE ,  VI S 


*  Satellites  *3  and  #4  * 

tseeeeeeseeseseeeeeeeseeeesssseseeeeeeeeeeeeesesseeeeeeeeeeeesseeeeesee 

PHI  3  ((B3(l)«B4(l))f(B3(2)*B4(2))t(B3(3)*B4(3)))/(R3(4)*B4(4)) 

IF  (PHI. GT.( l.ODO))  THEN 
PHI  3  l.ODO 

ELSE  IF  (PHI. LT.(-l.OOO))  THEN 
P!1I  3  -l.ODO 
ENDIF 


KKiTE(15,110)COUNT*SIZE,VIS 

20  CONTINUE 

CLOSE (UNIT  =  1) 

CLOSE (UN IT  =  2) 

CLOSE (UNIT  =  3) 

CLOSE! UNIT  =4) 

CLOSE! UNIT  =  10) 

CLOSE (UNIT  =11) 

CLOSE (UN IT  =12) 

CLOSE! UNIT  =13) 

CLOSE! UNIT  =14) 

CLOSE ( UNIT  =  15) 

RETURN 

110  FOBMAT(I10,5X,F12.5) 

END 


aCDHUt-’l  1»C.  OLAKU 


m 

* 


PURPOSE:  The  purpose  of  this  subroutine  is  to  parabolically  blend  * 
the  points  of  the  visibility  function.  The  equation  of  ♦ 
the  blended  curve  will  then  be  solved  for  real  roots.  If  ♦ 
a  real  root  exists,  it  represents  either  a  rise  or  set  » 
tiae  for  the  objects  of  interest.  * 

* 


V.\RI.4BLES: 

INPUT: 

STATE 

STEPS 

T 

VIS 


specifies  a  file  unit  nuaber 
nuaber  of  steps 

chronological  tiae  over  tiae  span 
value  of  the  visibility  function 


*  LOCAL  VARI.ABLES:  ♦ 

*  A,B  -  variables  used  in  the  solution  of  the  closed  fora  * 

*  cubic  equation  (Escobal,  432-434)  * 

*  ALPHA  -  coefficients  of  the  cubic  ,  * 

*  CHECK  -  slightly  later  tiae  than  cross.  The  value  of  check  ♦ 

*  is  used  to  deteraine  is  a  horizon  crossing  is  a  * 

*  rise  or  set.  ♦ 

*  '  COSPHI  -  cosine  of  the  angle  phi  used  in  solution  of  the  * 

*  cubic  equation  * 

*  COUNT  -  do  loop  counter 

*  CROSS  -  tiae  of  a  rise  or  set  * 

*  D120  -  120  DEG  changed  to  radians  * 

*  D240  -  240  DEG  changed  to  radians  * 

*  DELTA  -  see  closed  fora  solution  of  cubic  equation  * 

*  EVENT  -  1  if  the  crossing  is  a  rise;  0  if  the  crossing  is  a  * 

*  set  * 

*  EO  -  see  closed  fora  solution  of  cubic  equation  * 

*  INC  -  increaent  of  time  T  ♦ 

*  0  -  coefficient  of  cubed  variable  in  cubic  equation  * 

*  P  -  coefficient  of  squared  variable  in  cubic  equation  * 

*  FKI  -  angle  used  in  calculation  of  roots  of  cubic  equation  * 

*  PI  -  aatheaatical  constant  pi  * 

*  Q  -  coefficient  of  the  first  degree  variable  in  the  * 

*  cubic  equation  ♦ 

*  QTNT  -  portion  of  the  quadratic  equation  under  the  radical  .* 

*  R  -  constant  value  in  the  cubic  equation  • 

*  ROOTl,ROOT2,ROOT3  -  possible  roots  of  the  cubic  equation  * 

*  ,  SINPHI  -  sine  of  the  angle  phi  used  in  solution  of  the. cubic  * 

*  equation  • 

*  Z1,Z2,Z3  -  variables  used  in  detersining  roots  of  the  cubic  * 

*  .  equation  ,  * 

*  '  '  ♦ 

*  REFERENCES:  Adaas  and  Rogers,  278-284  * 

*  Alfano,  Negron,  and  Moore,  4-5  * 

*  Escobal,  431-433  * 

««»»*«««*«»*«aa**a*«««*«»*****»«s**»«*«*y««*«*»*«*«**»«»***«*«*»«*»«*«* 

SUBROUTINE  BLEND (STATE, STEPS,!, VIS) 


INTEGER  COUNT, STATE, STEPS 

REALMS  A , ARGl , ABG2 , B, ALP2 (0:3), ALPHA( 0:3), CHECK , COSPHI .  CROSS 
REAL*8  D120,D240, DELTA, EO, EVENT, INC, P,PHI, PI, Q, QTNT, VAU.VALB 
REAL*8  B,ROOT1,ROOT2, BOOTS, SINPHI, T(0;STEPS+2) 

REAL*8  V1S(0:STEPS>2),Z1,Z2,Z3 

*«**««*sss***«s*«*s*«*s«***««*s*s«»««s*««*s*««*«*********s**«s**«s*ss*« 

*  CONSTANTS  * 

sataas**ss«sas*ss*«s>^ss*****s«a*«*s«s«s*««*«s******»****«««ss«*sw-3s**«* 

PI  =  DACOS(-l.ODO) 

D120  =  120.0D0»(PI/180.0DO) 


* 

t 


The  folloKin?  lines  assicn  values  to  the  first -and  last  points  in  ^ 
in  the  visit. iiitv  curve.  Thest  points  fail  imEediatciy  befere  and  * 
'  iranu-diateiv  after  the  endpoints  of  interest.  This  insures  nc  * 

portion  of  the  curve  will  be  onatted  »-nen  testing  for  hcriccn  * 

*  crossings. 


T(0)  =  1:1) 

visto)  =  vis;i! 

T £TEPS+‘' * '  =  T '  STE!"-- 1 ' 

VIS' STEF'S*!  ;  =  VIStSTEI'S*!  ! 


*  In  the  fcllcKin?  do  loop,  the  coefficients  of  tiio  cubic  ociiiaticn  * 

*  are  deterEined  usinc  parabolic  blending.  * 

DO  10  COIKT  =  1, STEPS 

ALPHA ( 3 )  =  ( - .  50D0 ;  ♦V i  S  ( COUNT- 1 )  +  { 1 . 50D0 ! *VI S ( COUNT )  - 1 . 5CD0  '  ^ 
+  VIS(COUNT+l)+( .50D0)*VIS(COUKT*2) 

ALPHA ; 2 }  =  VIS'  COUNT- 1 ) t ( -2 . 50D0 ) *VI S { COUNT )  + ! 2 . ODO ! * 

+  VIS! COUNT+ 1 ) - ( - . 50D0 i *VI S { COUKT+2 ) 

.tLPHA (1 !  =  ! - . 50D0 )  * V I S (COUNT- 1 ) + ( . 5D0 1 *VI £ ( COUNT* 1 ) 

ALPHA! C)  =  VIS! COUNT) 

zz*z*»#»******z*»**z»x»»*z»Z*  ********  zz**^!***»»**»*z**z*rz*r*«»-z»»*l:*»z 

In  the  event  that  the  ec,uation  is  a  quadratic  equation  instead  of  * 

*  a  cubic  equation,  the  following:  lines  determine  the,  roots  of  the  * 


*  quadratic. 

******************** »»*********>:*******»********** ********************* 


IF  iDAESlALPHAOl  ).LE.  !  0.  OOOOOIDO)  .AND.  DABS!  ALPHA!  2 1 ).  LE. 
:  0 .  OOOCO IDO ; .  and .  D.IBS  !  alpha  1 1 )  l .  LT .  ( .  OOOOO IDO  ) )  then 
ROOT!  =  -i.CDO 
do:":  =  -i.ODC 
ROOTc  =  -l.ODO. 


ELSE  IF  !  DABS 'ALPHA!. 3  i  ).LE.  !0.  OO'OOOIDO) .  AND.DAES!  ALP.HA!  2  )  I .  L£. 
(O.OOOCClDOl  •;  THEN 

R'DOTl  =  -ALPHAIOl/.'.LF'HUl ) 

ROCTl:  =  -l.ODO 
.ROOT3  =  -l.ODO 


GOTO  5 

ELSE  IF  !D\BS!ALrHA;3j;.LE.{0.000001D0))  THEN 

QTNT  =  'ALPHA!  1 1*»2.0D0)-!  1.0, DO»ALPKA(2)*ALrHA!C'! ) 

IF  iOINT.GE. !0.0D0) !  THEN  • 

FOOT 1 '  =  ( -ALPHA ( 1 ) +DSQRT ! QTNT ) ) / ' 2 . ODO* ALPHA ( 2 ) ) 
f.'OOT2  =  !  -ALPHA !  1 )  -DSQRT !  QTNT )  1  / ( 2 .  ■0DC*.'.Lr'HA !  2  )  .' 

'  ROOT .3  =  -l.ODO 
GOTO  5 
ELSE 

.ROOTl  =  -l.ODO. 

ROOT2  s  -l.ODO 
ROOT3  =  -l.ODO 

GOTO  5  ' 

END!  F 
ENDIF 

x**************************************************'* ******************* 

♦  The  following  is  a  closed  form  solution  to  the  cubic  equation..  * 

«*******************x**»******x*************«»****»*«****************** 

P  =  .\LPHA(2)/ALPHA(3) 

Q  =  ALPHA(1)/ALPH.A!3) 

R  =  .ALPHA{0)/.ALPHAI3) 

A  =  (1.0D0/3.0D0)*({3.0D0*O)  -  (P»*2.0D0)) 

B  =  ( 1 . ODO/2  7 . ODO )* ( ( 2 . ODO* ( F**3 . ODO ) ) - ( 9 . 0Dd*P*Q )  +  ( 27 . 0D0*R ) ) 
DELTA  =  ((A*»3.0D0)/27.0D0)  *  ((B*»2.0D0)/r.0D0) 


IF  (ARGl.LT. (O.ODO; )  THEN 

VALA  =  -((DABS(ARG1))**(.33C32DC)^ 

ELSE 

VALA  =  ARG1**(.33333D01 
ENDIF 

IF  ; ARG2.lt. (O.ODO)  )  THEN 

VALB  =  - ( { DABS ( ARG2  ))**'.  33233DO i 
ELSE 

VALE  =  ARG2**(. 33332D0: 

ENDIF 

Z1  =  VALA  .+  VALB 
R(X)T1  =  Z1  -  {P/3.0D0) 

ROOT2  =  -l.ODO 
ROOT3  =  -l.ODO 

ELSE  IF  (DELTA. EG. (O.ODO))  THEN 
ARGl  =  -B/2.0D0 
ARG2  =  B/2.0D0 

IF  ('JIG1.LT.( O.ODO))  THEN 

Z1  =  -2.0D0*((DABS(.ARG!;)*»(.33333D0)) 

ELSE  IF  (ARGl. GT. (O.ODO))  THEN 

Z1  =  2.0D0»((DABS(ARG1 !)*»{. 33232DO)! 

ENDIF 

IF  {ARG2.LT. (O.ODO) )  THEN 

22  =  -( (D.ABS{.\RGL  "*( . 32323DO i  ) 

ELSE  IF  (AiiG2.GT.  (O.tOO) )  T.HEN 
Z2  =  .ARGL**!  .33333D0; 

ENDI” 

ROCTl  =  21  -  fP/3.0L'0; 

ROOT2  =  Z2  -  (P/3.0D01 
.ROOT  3  =  -l.ODO 

ELSE  IF  (DELTA.lt. (O.ODO))  THEN 
.ARGl  =  -(.A»*3.0D0)/27.0D0 
ARG2  =  -(A/3.0D0) 

IF  (  ARC  1  ■.  IT.  i  O.ODO).  OR.. ARGG.LT.i  O.ODO)  THEN 
.ROOTl  =  -l.ODO 
root:  -  -l.ODO 
ROOTS  =  -l.ODO 
ELSE 

EO  =  2.0D0*DSQRT(ARG2) 

COSPHI  =  -B/{2.0D0*DS0RT(.ARG1)) 

SI.NPHl  ,=  DS(3RT(  1. ODO-( COSPHI *-»2.0D0)) 

PHI  =  D.\T.A.N2  ( SI  NPHI,  COSPHI) 

Z1  =  EO»DCOS(PHI/3.0D0) 

Z2  =  EO*DCOS((PHI/3.0D01+D120) 

Z3  =  EO*DC’OS((PHI/3.0DO)+D2IO) 

ROOTl  =  ZT  -  {P/3.0D0) 

ROOT2  =  Z2  -  (P/3.0D0) 

ROOT3  =  Z3  -  IP/3.0D0) 

ENDIF  . 

E.NDIF 

5  CONTI.NUE 

*  If  a  positive  root  is  found,  parabolic  blending  is  used  to  deter-  * 

*  a  cubic  equation  based  on  the  tines  corresponding  to  the  visibility* 

*  values  of  interest.  The  value  of  this  cubic  equation  evaluated  at  ♦ 

*  the  root  of  the  equation  above  ( ROOTl , KOOT2 ,  or  ROOT3  above)  to  * 

*  find  tlie  e:;art  time  of  crossinc.'  The  same  equation  is  evaluated  * 


IF  (ROOT1.GE.0..4ND.ROOT1.LE.1)  THEN 

ALP2(3)  =  (“.SODOntCOL'NT'D  +  d.SODOjnt  COUNT  )  +  (-1.50D0)* 
f  T ( COUNT+ 1 ) + ( . 50D0 ) *T { COUNT+2 ) 

ALP2(2)  =  T(COUNT-l)+(-2.50D0)*T(COUNT)+(2.0D0)« 

T ( C0UNT+ 1 ) + ( - . 50D0 ) ♦T ( COUNT+2 ) 

ALP2(1)  =  {-.50DO)*T( COUNT-1 )+(.5DO)*T(COUNT+l) 

ALP2(0)  =  T( COUNT) 

CROSS  =  ({ROOT1**3.0DO)*ALP2(3))+((ROOT1*«2.0DO)*ALP’(^>))+ 
(ROOT1*ALP2(1))+ALP2(0) 

ROOTl  s  ROOTl+INC 

CHECK  s  ((ROOT1**3.0DO)*ALPH.A(3))+((ROOT1«2.0DO)*ALPH4{2))  + 
(ROOTl*ALPH.4(  1 )  )+ALPHA(0) 

IF  (CROSS. GE.(T(COUNT)). AND. CROSS.LE.(T(COUNT+l)))  THEN 
IF  (CHECK. GT. lO.ODO) )  THEN 
WRITE(25,100)CBOSS 
EVENT  =1.000 
WRI TE ( 60 ,*) CROSS , EVENT 
WR I TE ( STATE ,*) CROSS , EITNT 
ELSE 

V(R1TE(25,110)CBOSS 
EVENT  =  O.ODO 
WRI TE ( 60 , t ) CROSS , EVENT 
WR I TE ( STATE ,♦) CROSS , EVENT 
ENDIF 
END  IF 
ENDIF 

IF  { ROOT2.GE.O. AND. ROOT’ -LE.l)  THEN 

ALP2 ( 3 )  =  ( - . jODO ) *1 { COUNT- 1 ) + { 1 . oODO ) *T { COUNT ) + { - 1 . 50DO ) * 

T ( COUNT+ 1 )+(. 50DO ) *T (COUNT+2 ) 

ALP2(2)  s  T(COUNT-l)+(-2.50D0)*T{COUNT)+l2.0D0)*, 

T ( cot NT+ 1 )+(-. oODO ) *T ( COUNT+2 ) 

Ar.P2(l)  =  {-.50D0)*T( COUNT-1 )+(.5D0)*T(C0UNt+lJ 
ALP2(0)  =  T( COUNT) 


CROSS  =  ({R00T2**3)*ALP2(3))+{(E00T2»*2)*ALP2(2))+(R00T2» 

+  ALP2(1))+ALP2(0) 

R0OT2  =  R00T2+INC 

CHECK  =  ((ROOT2**3)*ALPH.A(3))+((ROOT2*»2)»ALPHA(2))+(ROOT’» 

+  ALPH.A(1))+ALPH.A(0) 

IF  (CROSS.GE.(T(COUNT)).AND.CROSS.LE.  (KCOUNT+l)))  THEN 
IF  (CHECK.GT. (O.ODO) j  THEN 
WRITE(25,100)CROSS 
El'ENT  =  l.ODO 
KR I TE ( 60 ,*) CROSS , EVENT 
WRI TE ( STATE ,*) CROSS , EVENT 
ELSE 

WRITE! 25, 110 )CROSS 
EVENT  =  O.ODO 
WRITE! 60 ,» )CROSS , EVENT 
WHI TE ( STATE ,*) CROSS , EVENT 
ENDIF 
ENDIF 
ENDIF 


IF  (ROOT3.GE.O.AND.ROOT3.LE.1)  THEN 

ALP2(3)  s  (-.50DO)*T(COUNT-l)+(1.50DO)*T(COUNT)+(-1.50DO)* 
♦  T(COUNT+1)+(.50DO)*T{ COUNT+2) 

ALP2(2)  =  T(COUNT-1)+(-2.50DO)*T(COUNT)+(2.0DO)* 

+  T(C0UNT+1)+(-.50D0)»T(C0UNT+2) 

- - .UP2  { 1 )  =  ( - .  50D0 )  «T ( COUNT-1 ) + ( .  5D0  >*Tt  rnnKT* 


+  ALP2(1))+ALP2(0) 

R00T3  =  R00T3+INC 

CHECK  =  ((ROOT3**3)*ALPHA(3J)+{(ROOT3**2)*ALPHA(2))+(ROOT3* 

+  ALPHA(1))+ALPHA(0) 

IF  (CROSS.GE.(T(COUNT)).AND.CROSS.LE.(TtCOLNT+l)))  THEK 
IF  (CHECK.GT.(O.ODO))  THEN 
WRITE! 25, 100 )CROSS 
El'ENT  =  l.ODO 
VfRI  TE  ( 60  ,*  )CROSS ,  EVENT 
TO I TE { STATE , ♦ ) CROSS , EVENT 
ELSE 

TOITE(25,110)CROSS 
EVENT  =  O.ODO 
WR I TE ( 60 ,*) CROSS , EVENT 
KH I TE ( STATE ,*) CROSS , EVENT 
ENDIF 
ENDIF 
ENDIF 


10  CONTiNLt 


100  FORMAT (r 12. 5,'  SEC  -  RISE’) 
no  FORMAT! F12. 5,’  SEC  -  SET’) 

RETTON 

END 


swonuciiAC. 

• 

* 

PURPOSE:  The 

purpose  of  this  subroutine  is  to  accomplish  a 

* 

* 

check 

for  constant  visibility  between  a  ground  station 

* 

* 

and 

2 

satellites. 

* 

* 

VARIABLES: 

* 

* 

IXaT: 

' 

* 

* 

¥ 

STATE 

- 

File  number  specifier 

* 

» 

LOCAL  VARIABLES: 

* 

« 

CHANGE 

- 

specifies  when  coaauncation  is  broken  or  re¬ 

* 

* 

established 

* 

« 

COUNT 

- 

do  loop  counter 

* 

* 

E1,E2 

- 

Specifies  rise/set  time  for  risibility  function 

* 

* 

a  ground  station  and  2  satellites 

* 

* 

EVENT 

- 

Array  of  all  rise/set  times 

* 

* 

IS 

- 

do  loop  counter . 

t 

* 

RS1,RS2 

- 

1  if  the  event  is  a  rise;  0  if  the  event  is  a  set* 

* 

SLT1 

- 

Sun  of  RSI  and  RS2  at  a  given  event 

* 

T 

- 

tipe  zero  of  any  selected  visibility  function 

* 

-  in  this  subroutine  T  is  a  place  holder 

* 

TEMP 

- 

Temporary  value  used  in  sort  loop 

t 

TOTAL 

- 

Total  number  of  events 

* 

TOTl 

- 

Number  of  events  between  Ground  Station  f?  and 

* 

Satellite  #1 

* 

T0T2 

- 

Number  of  events  between  Ground  Station  9?  and 

* 

Satellite  ffZ 

* 

TT 

- 

Holds  all  event  times  and  the  correspcndi:ig 

* 

status  of  the  visibility  function  for  each 

s 

viewing  path 

* 

V 

- 

the  first  value  of  the  selected  v^ibility 

* 

function 

* 

X 

- 

do  loop  counter 

* 

¥*««««»X««»*»******»«*««««*«**«««*«»«*««***«**«*»*»»***»*»»»»**»«*«** 

SUBROUTINE  GSSCVC2 { STATE ) 

REAL*8  El(0:100), £2(0:100), E\*ENT{0:500) 

REAL*8  RS1|0:100),RS2(0:100),SUM,T,TEMP,TT;0: j00,2),V 
lNTEGER  COUNT,  CH.ANGE,  CHECK,  IX,  STATE,  TOTAL,  TOT  l,TOT2,X 

OPEN ( US I T= I 1 . Fi LE* ’ STATl . DAT ’ , STATUS^ ’OLD ’ ) 
OPEN(UNIT=42,FILEs’STAT2.DAT’,STATUS=’OLD’ ) 

OPEN { LSI T=60 , FI LE* ’ TRACKER . DAT ’ , STATUS^ ’ OLD ’ ) 

REWIND{UNIT=60) 

OPEN(USlT=10,nL£s’VISl.DAT’,STATUS=’OLD’) 

OPLN  ( UN I T= 1 1 , FI LE= ’ V I S2 . DAT ’ , STATUSs ’ OLD ’ ) 

VRlTE(*,*)’The  coaputer  is  in  the  process  of  deteraining  whether’ 
'^■RITE(*,*)’or  not  continuous  coaaunicaticn?  is  possible,  Khen  ’ 
KRITE{*,*)’the  prograa  is  finished  runni  "  results  can, be  found’ 
WRITE! *,*)’ in  the  file  RISESET.DAT.’ 

*  GROUND  STATION  AND  SAVL-LITE?!  * 

«**»«a***s»*«**«s*«««««««s««*«««»«««««*«***««««***»***«s»«****s»«»»** 

TOTl  *0 

DO  20  COUNT  s  1,100 

BEAD ( STATE ESD= 2 5 ) E 1 { COUNT J ,RS 1 { COUNT ) 

TOTl  s  COUNT 
20  COfTTINUE 


READ(10,*)T,V 
IF  (V.GT.IO.ODO))  IKES 
RSIIO)  =  l.ODO 
ELSE 

RS1{0)  =  O.ODO 
ENDIF 

ELSE  IF  (RSl(l).EQ. (O.ODO))  THEN 
RSl(O)  =  l.ODO 
ELSE 

RSKO)  =  O.ODO 
ENDIF 

*  GROUND  STATION  .AND  S.ATELLITE  »2  ♦ 

TOT2  =  0 

DO  30  COUNT  =  1,100 

8LAD( STATE+1 , * , END=35 ) E2 1  COUNT ) , ES2 ( COUNT ) 

TdT2  =  COUNT 
30  CO.NTINUE 
35  E2(0)  =  O.ODO 

IF  (TOT2.EQ.0)  THEN 
RE.ADI11,*)T,V 

IF  (V.GT. (O.ODO))  THEN 
RS2(0)  s  l.ODO 
ELSE 

RS2(0)  s  O.ODO 
ENDIF 

ELSE  IF  (RS2(1).EQ. (O.ODO))  THEN 
RS2(0)  s  l.ODO 
ELSE 

8S2(0)  s  O.ODO 
ENDIF 

*  The  follovin?  calculates  the  total  nuaber  of  events  between  the  ♦ 

*  ground  station  and  two  satellites.  » 

»»*»**»t*«**»«****a«*««:«**««*««>«4'«***«««»**»»*«»**«»»4«**«»»*»**«**»» 

TOTAL  =  TOTUTOT2 
CHECK  =0 

DO  15  COUNT  =1,500 

READ ( 60 END= 18 ) EVENT ( COUNT ) 

CHECK  s  COUNT 
15  CONTINUE 
18  CONTINUE 

IF  (CHECK. EO.O)  THEN 

WRITE(25,*) 'There  were  no  visibility  opportunities  with  this 
configuration.* 

GOTO  no 
ENDIF 

*  The  following  orders  all  the  events  froa  the  saallest  to  the  * 

*  largest.  * 

sssssssasssssssssssssssssssssssssssssssssssssstsssssssssssssss******* 

DO  60  COUNT  »  1, TOTAL-1 
DO  50  X  *  1, TOTAL-1 

IF  (E\'E.NT(X*l).LT.EVEh*T(X))  THEN 
TEMP  s  EA’ENT(X+1) 

EVENTlX+1)  »  EVENT(X) 

EVEKT(X)  =  TEMP 
E.NDIF 

50  CONTINUE _ 


*  The  following  deteraines  the  status  of  coaouni cat ions  for  each  * 

*  possible  visibility  path  at  each  event.  * 

****«*aa****#»*««**:^««*«««*««««a«»**t:***«««*«««*****»****««*a*****4:»* 

DO  80  COtNT  =  0, TOTAL 

IF  (TOTl.EQ.O)  THES 
TT(COUNT,:)  =  RSKO) 

GOTO  80 
ENDIF 

DO  70  IS  =  O.TOTl  -  1 

I F  ( EVENT { COUNT ) . GE . El ( I S ) . AND . EVENT ( COL  NT ) . LT . El ( I S+ 1 ) )  THEN 
TTICOUNT.IJ  =  RSI (IS) 

GOTO  80 

ELSE  IF  iEVENT( COUNT ).EQ. El (IN+D)  THEN 
TT(COUNT,l)  =  BSKIS+l) 

GOTO  80 

ELSE  IF  ( EVENT ( COUNT ).GT, El (IN+1).. AND. I N.EQ.(TOTl-l))  THEN 
TT((^’KT,1)  =  8S1(IN>1) 

GOTO  80 

ENDIF 
70  CONTINUE 
80  CONTINUE 

DO  120  COUNT  =  0, TOTAL 

IF  (TOT2.EO.O)  THEN 
n{COUNT,2)  =  RS2(0) 

GOTO  120 
ENDIF 

DO  110  IN  =  0,TOT2  -  1 

i F  { EVENT i cot  NT  1 . CE . E2 ( I N ; . AND . EVENT . COUNT ) . LT . E2 ( I N+ 1 ) )  THEN 
TT1C0UN:,2)  =  RS2(IM 
GOTO  120 

ELSE  IF  ( EVENT!  COUNT  ).EQ.E2 1  IN+D)  THEN 
TT(COUNT,2)  =  RS2(IN+D 
GOTO  120 

ELSE  IF  {EVEN'T(COUNT).GT.E2(IN+l).AND.IS'.Eg.(TOT2-D)  THEN 
TT(COUNT,2)  =  RG2(1N+1) 

GOTO  120 

ENDIF 
■10  CONTINUE 
120  (ONTINUE 

CHANGE  =  1  . 


*  The  riext  DO  loop  deteraines  if  constant  coaauni  cat  ion  is  .* 

♦  possible.  If  the  sua  is  eiqUri  to  zero,  coaaunication  is  not  * 

*  possible.  If  the  sua  is  greater  than  or  equal  to  one,  coaauni-  * 

♦  cation  is  possible.  * 


DO  130  COUNT  =  O.TOTAI. 

S171  =  TT(COUNT,l)  +  TT(COUNT,2) 

IF  (SL^.EQ.(0.0D0)..AND.CH.ANGE.FQ.1)  THEN 

KRITE(23,*)’CoBBunications  broken  at  ’ .EVENT (CODT) 

CHANGE  =0 
ENDIF 

IF  (SUM.GE.(l.ODO).ASD.CHANqE.EQ.O)  THEN 

WHITE(25,*)’CoBBunication8  re-established  at  ’ ,E\*EST( COUNT) 
CHANGE  ?  1 
ENDIF 
130  CONTINUE 
140  CONTINUE 


CLOSE M'SIT=10 ,  STATl’S-  ’ DELETE ’ ) 
CLOSE I  IS I T= 1 1 , STATLS= ’ DELETE ’ ) 

RETURN 

END 


tift^t^^-^i-IHH^t*************^**************^*****:*****  '.l'5i!*i?;*i*******v*» 

*  SUBROUTINE  CVCi  * 

*  '  * 

*  PURPOSE:  The  purpose  of  this  subroutine  is  to  accoBpllsh  a  * 

■*  check  for  constant  visibility  between  3  satellites.  * 

*  For  visibility  to  be  maintained,  at  leapt  tso  * 

*  paths  visibility  muuc  be  maintained  at  a’l  times.  * 


V.aRIaBLES: 

INPUT. 

KIND 


STATE 


-  spijcifies  whether  3  satellites  are  invol\*d  (1) 
or  a  ground  station  and  3  satellites  are 
involved  (0) 

•  File  number  specifier 


LOCAL  VARIABLES: 
CHANGE  -  soi 


FVF.NT 

IN- 


TEMP 

TOTAL 

TOTl 


*  CHANGE  -  specifies  when  conauncation  is  broken  or  re-  .* 

*  established  * 

*  COUNT  -  do  loop  counter  * 

*  E3  -  Specifies  rise/net  tine  for  visibility  function  * 

'*  between  3  satellites  or  a  ground  station  and  * 

*  3  satellites  ♦ 

*  FVF.NT  -  Array  of  ail  rise/set  tines  * 

*  IN  „  -  do  loop  counter  * 

*  RS1,..RS3  -  1  if  the  event  is  a  rise;  0  if  the  event  is  a  set* 

*  '-UM  -  Sun  of  RSI . . .  RS3  at  a  given  event  * 

*  r  -  tine  zero  read  fron  file  of  visibility  data  * 

*  -  in  this  subroutine,  1  serves  only  as  a  place  * 

*  holder  * 

*  TEMP  -  Temporary  value  used  in  sort  loop  * 

*  TOTAL  -  Total  number  of  events  * 

*  TOTl  -  Number  of  events  between  Satellite  ?1  and  * 

*  Satellite  *2  or  Ground  Station  ??  and  * 

*  Satellite  * 

»  T0T2  -  Number  of  events  between  Satellite  =2  and  * 

*  Satellite  *3  or  Ground  Station  ??  and  * 

*  Satellite  '  * 

*  T0T3  -  Number  of  events  between  Satellite  #1  and  * 

*  Satellite  =3  or  Ground  Station  ??  and  * 

*  '  Satellite  =3  ♦ 

*  TT  -  Holds  all  event  times  and  the  corresponding  ♦ 

*  status  of  the  visibility  function  for  each  ♦ 

*  '  viewing  path  '  , 

*  ,  V  -  value  of  the  chosen  visibility  function  at  "time  * 

*  equal  to  zero  * 

*  X  -  do  loop  counter  * 

**********S******4?******»*****«**«<*»«*«**S****************4:****#*4:*** 

SUBROUTINE  CVC3( KIND, STATE) 

REAL*8  El  { 0 : 100 ) ,  E2  ( 0 : 100  j , E3  { 0 ;  100 ) ,  El^ENT ( 0 : 500 ) . RSI  ( 0 ;  100 ) 
REAL*8  RS2 { 0 : 100 ) , RS3 { 0 : 100 ) , SUM, T .TEMP , TT { 0 : 500 , 3 ) , V 
INTEGER  CHANGE,CCUNT,IS,KIND,STATE,TOTU,TOT1,TOT2,TOT3,X 

IF  {KIND.EQ.il  THEN 
KRITEI*,*) 

HRITfil*,*)*The  computer  is  in  the  process  of  determining* 

;  WBITE(*,*) 'whether  or  not  constant  communications  is  * 
WRITE(*,*)’ maintained.  When  the  program  has  finished*. 
WBlT£(*,*)*ruhning,  results  may  be  found  in  RISESET.DAT.’ 


,0PEN{UKITs41.FILE=*STAT1.DAT''  .STATUSs’OLD’  ) 

.A-—;-..!,.  A.. — ^  Ma  imi_:i  g  mul  t  ' 


JcL'-A-  .-vV-. 


OPEN  ( UNiTrbO ,  FI  IE=  ’  TkAC KEk .  DAI  ’ , STaTI  S>=  ’  OLD  ’ ) . 

IF  (KIND.EQ.O)  THEN 
REVIND(L’N1T=60) 

ENDIF 

OPEN(USIT=10,FILE='VIS1.DAT’,STATIS=’OLD’ ) 
OPEN(USIT=n,FILE=’VlS2.DAT’  ,STATUS=‘OLD’ ) 

OPEN ( UN I T= 1 2 , F I LE= ’ V I S3 . DAT ‘ , STATUS= ’ OLD ’ ) 

*  The  folloKing  DO  loops  read  the  rise  and  set  tiaes  for  each  vis-  * 

*  ibility  path.  .  * 

***************t*t*t*******************************t****9t****t****** 

*  SATELLITE  ?1  AND  SATELLITE  »2  * 

*  GROUND  STATION  AND  SATELLITE  «1  * 

TOTl  =  0 

DO  20  COUNT  =1,100 

READ! STATE , * , ESD=25 ) El ( COUNT ) , RSI ( COUNT ) 

TOTl  =  COUNT 
20  CONTINUE 
25  EKO)  =  O.ODO 

*  If  there  were  no  rise  or  set  tines,  the  file  containing  the  vis-  * 

*  ibility  data  aust  be  referenced  to  deteraine  if  the  two  objects  * 

*.  of  interest  are  always  in  view  or  never  in  view.  * 

***«a**«s;>**«s***s***s*s***««*s«s«s**s*s*«s««s*s«sa***ss*s»*s**sss**s 

IF  (TOTl.EQ.O)  THEN 
READ(10,*)T,V 
IF  (V.GT. (O.ODO))  THLN 
BS1(0)=1.0DO 
ELSE 

RS1(0)=0.0D0 

ENDIF 

****4E«*«**x**«*««««««««4;«»«««*««X********S*««*S«**a**S*******4:*»***»* 

*  The  following  line  checks  the  first  condition  of  the  first  rise  * 

*  or  set  time  then  assigns  the  appropriate  condition  to  time  zero.  * 

ELSE,  IF  lESKD.EQ. (O.ODO))  THEN 
BSKO)  =  l.ODO 
ELSE 

BSKO)  =  O.ODO 

ENDIF 

tS^:**************************************************************’**** 

*  '  SATELLITE  »1  AND  SATELLITE  f3  * 

*  ground  station  AND  SATELLITE  »2  ♦ 

TOT2  =0 

DO  30  COUNT  =  1,100 

READ(  STATE+1 ,  * ,  END=35 )  E2  ( COLTiT ) ,  BS2  ( COLUT) 

TOT2  =  COUNT 
30  CONTINUE 
35  E2(0)  =  0.01)0 

IF  (TOT2.EQ.O)  THEN 
READ(11,*)T,V 
IF  (V.Gl. (O.ODO))  THEN 
BS2(0)=1.0Dn 
ri 

BS2C0)=0.0D0 


*  ThE  following' line  checlcs  the  first  condition  of  the  first  rise  * 

*  pr  S(^t.  ti^  a^ropriste  condition  to  tine  zero.  *, 


Rsi:(0)  =  i.oDo 
ELSE 

RS2{0)  =:  O.ODO 

ENDIF 

*  SATELLITE  ?2  ASD  SATELLITE  #3  ,♦ 

*  GROUKD  STATION  AND  SATELLITE  ?3  * 

I  «***4;*i4;»*«««t*««4t«»#*»«*»4!«*4!*««»»*«»«***«4:««**»*****»**«****4:**:t:4!*** 

TOTS  =0 

DO  40  COUNT  =  1,100 

READ{  STATE+2 ,  * ,  ESD=45  )E3  ( COUNT )  .RS**  { COUNT ) 

TOT3  =  COUNT 
40  CONTINUE 
45  E3(0)  =  O.ODO 

IF  (TOT3.EQ.O)  THEN 
READ(12,*)T,V 
IF  (V.GT. (O.ODO))  THEN 
RS3(0)=1.0D0 

FT 

RS3(0)-0.0D0 

ENDIF 

t*******************^************************************************ 

*  The  folloKing  liite  checks  the  first  condition  of  the  first  rise  * 

*  or  set  tiae  then  assigns  the  appropriate  condition  to  time  zero.  * 
*»»***e«***»e»«»«es**«*«*«s««****««ss****»«*«***«**s«e**s**es«:i;*e**«* 

ELSE  IF  (RS3(1).EQ.(0.0D0))  THEN 
RS3(0)  s  I.ODO 
ELSE 

RS3(0)  s  O.ODO 

ENDIF 

T:««««***»*«*«***»*t*««*******««***»**4c****»«e«**«*e*»*******4:****«*** 

*  The  folloving  calculates  the  total  number  of  rise/set  events  * 

*  between  the  ground  station  and  the  three  satellites  or  between  * 

*  the  three  satellites.  ♦ 

TOTAL  =  TOT1+TOT2+TOT3 

DC  47  COUNT  =1,500 

read; 60 END=48 ) EVENT ( COUNT ) 

47  CONTINUE 

48  CONTINUE 

t***:^****************************************)************************ 

*  If  there  were  no  rise/set  tines  at  all,  the  following  message  is  * 

*'  written  to  the  final  output  file  and  the  prograa  ends.  * 

.  *««»i«S***«**e**«««****«4:«***««*««******S**S«****S«*S*«»**^«**«***««*« 

’  IF  (C0UNT.EQ.1)  THEN 

WRITE! 25,*) ’There  were  no  opportunities  for  viewing 
+  wtih  this  configuration. ’ 

GOTO  215 
ENDIF 

*  The  following  performs  a  sort  to  place  the  list  of  all  events  in  * 

*  order  from  smallest  number  to  largest.  * 

«««*««*««s«**«s*»****»*****s«s*s*s**»**sss«*«»««»**«»*««*sss«««**««** 

DO  60  COUNT  s  l,T0tAL-l  . 

DO  50  X  s  1,T0TAL-1 

.  :  IF  (EVENT(X+11.LT.EVEKT{X))  THEN^  ^  ^ 

-S' EVENT(X).  .  ■ 

EVENTCX)  *  TEMP 

nroip 

50  COSTINUE 
60  CONTINUE 


*  The  following  loops  determin'?  the  interaediate  status  of  v^^ibil-* 

♦  ity  at  each  event.  ♦ 

*****t************:$*****t************************t**t***t************ 

DO  80  COUNT  =  0, TOTAL 

IF  (TOTl.EQ.O)  THEN 
TT(C0UNT,1)  =  BSl(O) 

,  GOTO  80 
ENDIF 

DO  70  IN  =  O.TOTl  -  1 

IF  (EVENT ( COUNT ) . GE . E 1 ( I N I . AND . EVENT ( COUNT ) . LT . El { I N+ 1 ) )  THEN 
TT(COUNT,l)  =  RSKIN) 

GOTO  80 

ELSE  I F  ( EVENTI COUNT ) . EQ. El ( IN+1 ) )  THEN 
TT(CO''NT,l)  =  RS1(1N+1) 

GOTO  80 

ELSE  IF  (EVENT(C0UNT).GT.E1(IN+1).A.ND.IN.EQ.(T0T1-1))  THEN 
.  TT(COUNT,l)  =  RSKTOTl) 

GOTO  80 
ENDIF 
70  CONTINUE 
80  CONTINUE 


DO  100  COUNT  =  0, TOTAL 

IF  (TOT2.EQ.0)  THEN 
TT{COUNT,2)  =  BS2(0) 

GOTO  100 
ENDIF 

DO  90  IN  =  0,TOT2  -  1 

I F  ( EVENT  ( COUNT ) .  GE .  E2  ( I N ) .  .AND .  EVENT  ( COUNT ) .  LT .  E2  ( I N+ 1 ) )  THEN 
TT(COUNT,2)  s  RS2(IN) 

GOTO  100 

ELSE  IF  I  EVENT ( COUNT ).EQ.E2( IN+1 ))  THEN 
TT(COUNT,2)  =,RS2(1N+1) 

GOTO  100 

ELSE  IF  (EVENT(COUNT).GT.E2{IN+l)..AND.IN.EQ.(TOT2-l))  THEN 
TT(COUNT,2)  =  RS2(TOT2) 

GOTO  100 

ENDIF 
90  CONTINUE 
100  CONTINUE 

DO  120  COUNT  =  0, TOTAL 
r  Do  110  IN  =  0,TOT3  -  1 

IF  (EVENT(COUNT).GE.E3(IN)..AND.EVENT(COUNT).LT.E3(IN+l) )  THEN 
TT(COUNT,3)  =  fiS3(IN) 

GOTO  120 

ELSE  IF  (EVENT{COUNT).EQ.E3fIN+l))  THEN 
TT{COUNT,3)  s  BS3{IN+1) 

GOTO  120 

ELSE  IF  (EVE!JT{COUNT).GT.E3(IN+l).A.ND.IN.EQ.{TOT3-l))  THEN 


TT(COUNT,3)  *  RS3(TOT3) 

•  GOTO  120 

ENDIF 

no  CONTINUE 

120  CONTINUE 

CHANGE*! 

*  The  following  detcraines  if  continuous  co«Bunicsti<m  is  possible.* 

*  For  3  satellites t  the  sua  sust  be  greater  thsa  or  equal  to  2.  .  * 

*  For  a  ground  station  and  three  satellites,  the  sua  aast  be  * 

*  greater  than  or  equal  to  1.  • 

eeeeeeeeeeeeeeeeeeseeesseeesaeseeaessessssesassaakctMtMteeeeeseeeeee 

IF  (UND.EQ.l)  THEN 


ir  \  sun«  LiC.*  \  1  •  uuu  j  •nni/^v.iuuiuc.* £.<•<  A  /  mbit 

WRI TE { 25 , 200 ) EVENT ( COUNT ) 

CHANGE  =  0 
ENDIF 

.  IF  (SUM.GE.(2.0D0).AND.CHANGE.EQ.O)  THEN 

»  V»TIITE(25,210)EV'ENT(COUNT) 

CHANGE  =  1 
ENDIF 

130  CONTINUE 

ELSE  IF  {KIND.r.3.0)  THEN 
DO  140  COUNT  =  0, TOTAL 

SUM  r  TT(COUNT,l)  +  TT(COUNT,2)  +  TT(C0UNT,3) 

IF  (SW!.LT.I1.0DO).AND.CH.ANGE.E0.1)  THEN 
WHITE (25 , 200 ) EVENT ( COUNT ) 

CH.ANGE  =  0 
ENDIF 

IF  ( SUM. GE.(l.ODO). AND. CHANGE. EQ.O)  THEN 
VRI TE ( 25 , 210 ) EIHNT  { COUNT ) 

CHANGE  =1 
ENDIF 

140  CONTINUE 
ENDIF 

215  CONTINUE 

CLOSE ( UN I T=STATE , STATUS= ’ DELETE ’ ) 

CLOSE ( UNI T=STATE+ 1 , STATUS= ’ DELETE ’ ) 

CLOSE ( UNI T=STATE+2 , STATUSs ’ DELETE’ ) 

CLOSE! UNI T=60) 

CLOSE ( I NI T= 1 0 , STATUSs ’ DELETE ’ ) 

CLOSE ( UNI T= 1 1 , STATUSs ’ DELETE ’ ) 

CLOSE ( UN I T= 1 2 , STATUS= ’ DELETE ’ ) 

200  FORMAT! ’Couunications  broken  at  ’,F12.5,’  SEC’  ) 

210  FORMAT! ’CoBiunications  re-established  at  ’,F12.5,’  SEC’) 

RETURN 

END 


SUBROUTINE  CVC4 


* 

* 

*  PURPOSE: 

* 

* 

* 

*  VARIABLES: 

*  INPUT: 

*  STATE  -  File  nuaber  specifier  * 

*  i 

*  LOC.AL  V.ARI.ABLES:  < 

«  CHANGE  -  specifies  when  couuncation  is  broken  or  re*  * 

*  established  * 

*  COUNT  -  do  loop  counter  < 

*  E1,..,E6  -  S'^ecifies  rise/set  tiae  for  visibility  function  ^ 

*  '  bi.'  een  6  satellites  ^ 

♦,  EVENT  -  Array  of  all  rise/set  tiaes  < 

*  IN  -  do  loop  counter  < 

*  RS1,..RS6  -  1  if  the  event  is  a  rise;  0  if  the  event  is  a  set* 

«  SUM  -  Sun  of  RS1...RS6  at  a  given  event  * 

*  T  .  -  tiae  zero  read  froa  file  of  visibility  data  * 

*  -  in  this  subroutine,  T  serves  only  as  a  place  * 

*  holder  * 

*  TEMP  -  Tenporary  value  used  in  sort  loop  * 

*  TOTAL  -*  Total  nuaber  of  events  > 

*  TOTl  -  Nuaber  of  events  between  Satellite  #1  and  * 

*  Satellite  f2  , 

*  T0T2  >  Nuaber  of  events  between  Satellite  #l  and 

«  Satellite  v3 

*  T0T3  -  Nuaber  of  events  between  Satellite  ?!  and 

*  Satellite  #4 

*  TOTi  -  Nuaber  cf  events  between  Satellite  «2  and 

*  Satellite  v3 

*  T0T5  -  Nuaber  of  events  between  Satellite  #2  and 

*  Satellite  #4 

*  TOTS  -  Nuaber  of  events  between  Satellite  73  and 

*  Satellite  74 

*  TT  -  Holds  all  event  tiaes  and  the  corresponding 

*  status  of  the  visibility  function  for  each 

*  viewing  path 

*  V  -  value  of  the  chosen  visibility  function  at  tiae 

*  equal  to  zero 

*  X  -  do  loop  counter 

SUBROITINE  CVCKSTAH) 

SEAL*8  £1(0:50), E2(0:50),E3(0:50), £4(0:50), E5(0:50) 

SEAL*8  £6(0:50), E1^ENT(0:200),BS1(C:50),RS2(0:50),BS3(0:50) 
REAL*8  RS4(0:50),BS5(0:50),RS6(0:50),SUM,T,T£MP,TT(0:200,6),V 
INTEGER  CHANGE,COUST,IN,STATE,TOTAL,TOT1,TOT2,IOT3,TOT4,TOT5 
INTEGER  T0T6,X 

WHITE!*,*) 

WRITE!*,*) 'T)m  coapQter  is  ia  the  process  of  detersining* 

>  WRITE!*,*) ’whether  or  not  constant  coasuniclitiws  Is  * 

WRITS!*,*) ’asintained.  Wbea  the  prograa  lias  flaished* 
WRITE!*,*) ’mnniaf,  results  aay  be  found  la  RISISET.DiAT. * 

OPn(UHIT»41,riLEa’STAT1.0AT*  tSTATVSa’OLD’ ) 
0PES!UllIts42,riL£s»STAT2.DAr  ,STATUSa’OLD» ) 

OPfR!UIIlts43,  riLEa'STATS.  DAT*  .STATUS* ’OLD* ) 

OPEW(UNIT«44. FILES *  STAT4 . DAT' . STATUS* ’ OLD ’ ) 


The  purpose  of  this  subroutine  is  to  accoaplish  a 
check  for  constemt  visibility  between  4  satellites. 
For  visibility  to  be  aaintained,  at  least  3 
paths  of  visibility  aust  be  aaintained  at  all  tiaes. 


OPES{UN*IT=10,FILE=’VIS1.DAT’ ,STATL'S=’OLD’ ) 

OPEN  ( U!S  I  T=  1 1 ,  F I  LE=  •  V I S2 .  DAT  ’ ,  ST ATLSi  •  OLD  ’  J 
OPEN{UNIT=12,FILE=’VIS3.DAT’ ,STATUS=’OLD’ ) 
OPEN(UNIT=13,FILF='VIS4.DAT’ ,STATUS=’OLD’ ) 
OPES(USIT=14,HLE=’VIS5.DAT’,STATUS=’OLD* ) 
OPEN(L'SIT=15,FILE=*VIS6.DAT’  ,STATUS=’OLD’ ) 
»»*»»♦*»»♦***♦*♦**♦♦*»*»**»»»»**♦»*♦»*»*»•*****»***»***»***♦**♦***♦♦* 

*  SATELLITE  H  AND  SATELLITE  fZ  * 

***♦**»****«»***»****»*«*********»*»♦*»*»»*♦♦**♦**♦*»**♦**♦*♦♦*»♦♦♦* 

TOTl  =0 

DO  20  connr  =  i.ioo 

READ ( STATE ESD=25 ) El ( COUNT ), RSI ( COUNT ) 

TOTl  =  COLTT 
20  CONTINUE 
25  EKO)  =  O.ODO 

IF  (TOTl.EQ.O)  THEN 
READ(10,*)T,V 

IF  (V.GT. (O.ODO))  THEN 
RS1(0)=1.0D0 

FI 

RS1(0)=O.ODO 

ENDIF 

ELSE  IF  (BSKD.EQ.  (O.ODO))  THEN 
RSKC)  =  l.ODO 
ELSE 

RSKO)  s  O.ODO 
ENDIF 

*«*«««««««***««»«««««*«**«**S«*««<*«X«***««««*«i»«»*«**«**«***»*«*«*«» 

*  SATELLITE  *1  .AND  SATULITE  f3  * 

T0T2  *0 

DO  30  COUNT  =1,100 

RLAD  ( ST.ATE+ 1  END=35 )  E2  { COUNT ),  RS2  ( COUNT ) 

T0T2  =  COUNT 
30  CONTINUE 
35  E2(0)  s  O.ODO 

IF  (TOT2.EQ.O)  THEN 
RLAD(11,*)T,V 
IF  (v.<rr. (O.ODO))  then 
RS2(0)=1.0D0 

ELSE  _ 

BS2(0)=O.ODO 

ENDIF 

ELSE  IF  (RS2(1).EQ. (O.ODO))  THEN 
BS2(0)  «  l.ODO 
ELSE 

R82(0)  s  O.ODO 
ENDIF 

x******************************************************************** 

«  SATELLITE  »1  AND  SATELLITE  #4  * 

**«««**««*«««*««««««««**«*««*«««*««*««**«*«*«X«XX»***«**«**X***«X*«*« 

T0T3  *  0 

DO  40  COUNT  *  1,100 

BEA0(STATE^2,«,£ND«43)E3(C0UST),RS3(C0U!iT) 

TOTS  *  COUNT 
40  COSFilNUI 
49  £3(0)  »  O.ODO 

IP  (TOT3.tQ.O)  THIN 

>iADn2.*)T.y  .  . .  ■ ,  _ _ ; : .  ,, 


LNDIF 

t**********9*************t*****************t***t*t*****************%* 

*  The  following  line  checks  the  first  condition  of  the  first  rise  * 

*  or  set  time  then  assies  the  appropriate  condition  to  tiae  zero.  * 

ELSE  IF  (BS3(1).EQ.{0.0D0))  THEN 
RS3(0)  -  l.ODO 
rt  <;r 

ES3(0)  =  O.ODO 
EXDIF 

t*************t********t*******t*****************Jt********t********** 

*  SATELLITE  *2  AXD  SATELLITE  #3  * 

«**««****»««sss*«*a»«aa««««*a«**«»«***«***s**$¥«**«**s«*«*a««aa**aa** 

T0T4  =  0 

DO  50  COUNT  =  1,100 

READ( STATE+3 , * , ENOroa ) E4 { COUNT ) ,RS4 { COUNT ) 

TOT4  =  COUNT 
30  CONTINUE 
55  £4(0)  s  0.000 

IF  (TOT4.EQ.0)  THEN 
READ(13,*)T,V 
IF  (V.GT. (O.ODO))  THEN 
BS4(0)=1.0D0 
FT  9F 

RS1(0)=0.0D0 

ENDIF 

»  The  following  line  checks  the  first  condition  of  the  first  rise  * 

*  cr  set  tiae  then  assigns  the  appropriate  condition  to  tine  zero.  * 

ELSE  IF  (RS4(l).EQ.(0.0DO))  THEN 
RS4(0)  3  1.000 
ELSE 

QS4(0)  3  O.vIDO 
ENDIF 

S*********************************************’#********************** 

*  SATELLITE  »2  AND  SATELLITE  f4  * 

TOTS  »  0 

DO  GO  COUNT  a  1,100 

BEAD! STATE>4 ,*, ENOaGS ) E5 ( COUNT ), RS5 ( COUNT ) 

TOTS  *  COINT 
60  CONTINUE 
65  E5(0)  *  O.ODO 

IF  (TOT5.EQ.O)  THEN 
READ(14,«)T,V 
IF  (V.(7r. (O.ODO))  THEN 
BS5(0)*1.0D0 
ELSE 

B85(0)«0.0D0 

ENDIF 

aaaaaaaaaaasaaaasaaaaaaaaaaaaaaaaaaaacaaaaaaaaasaaaaaaaaaaaaa******** 

*  Tha  fqlloainf  lina  cfeacks  tba  first  coadition  of  :tha  first  riso  * 

*  or  sot  tiao  tboa  asslfas  tho  sppropriat*  conditioo  to  tlao  zoro.  * 
soaaaoaooaoaosssssaaaasoaaaaassaaasssssasaaaasasaaaaaosaosasat 

nn  IF  (ISS(l).IQ. (O.ODO) )  TID 
185(0)  •  1.000 
DLSf  . 

185(0)  a  O.ODO 
„JW01F 


DO  70  COUNT  =  1,100 

RL4D{  STATE>5 ,  * ,  E}{D=73 )  E6(  COUNT )  .BS6  ( COUNT ) 

TOT6  =  COUNT 
70  CONTINUE 
73  E6(0)  =  O.ODO 

IF  (TOT6.EQ.O)  THEN 
HLAD(15,*)T,V 
IF  (V.GT. (O.ODO))  THEN 
RS6(0)=1.0D0 
ELSE 

RS6(0)=0.0D0 

ENDI? 

«4!«»«*«»»«*4:««*«***«*«*4:«*«««***«««***»*««**«***»«»«*««t««***«»«««««»« 

*  The  folloving  line  checks  the  first  condition  of  the  first  lise  * 

*  or  set  tiae  then  assigns  the  appropriate  condition  to  tine  zero.  * 
*»sea*«**»**»*«««*«**«*«***«****««««*»««*a«a*«*«**«*«*»*»s«*«s«««»a«* 

ELSE  IF  (RS6(1).EQ. (O.ODO))  THEN 
RS6(0)  ::  l.ODO 
ELSE 

RS6(0)  =  O.ODO 
ENDIF 

«***«*a******««««««*«*a«»«*««««*«««**a««e«»«**e***»«*»*y»«*»««»»««««« 

*  Tlie  following  calculates  the  total  nuaber  of  rise/set  events  *' 

*  between  the  ground  station  and  the  three  satellites.  * 

TOTAL  s  TOTltTOT2TTOT3+TOT4+TOT54TOT6 

DO  75  COUNT  *  1,500 

RE-4D  ( 60 ,  ♦ ,  END*  7  6 )  EIXNT  ( COUNT ) 

75  CONTINTE 
75  CONTINUE 

«fs***a«»«««*«ti««««««««««*«*««e***«*««*«*$*«««***«**»»«****«**«*»**e« 

*  If  there  were  no  rise/set  tiaes  at  all,  the  following  aessage  is  * 

*  written  to  the  final  output  file  and  the  prograa  ends.  * 

IF  { COUNT. EQ.l)  THEN 

VRITE( 23,*) 'There  were  no  opportunities  for  viewing  with  this 
4  configuration.' 

GOTO  215 
ENDIF 

*********e««««***««««*«»«««**««s*«*««*««***«**«***«**««*****»**e***««. 

*  The  following  perforas  a  sort  to  place  the  list  of  all  events  in  * 

*  order  froa  saallest  nuaber  to  largest.  * 

«s«*s««*«*«*«*s«*s««*«s«««*«s*«ss*«s«««**sa*«s*«a**e*«*«««s********** 

DO  80  COUNT  «  1, TOTAL-1 
DO  78  X  »  l,T0TAL-l 

IF  (EVZST(X4l).LT.EVENT(X))  THEN 
TEJC*  a  EVEST(X41» 

EVCNTdtl)  «  Emr(X) 

El’ENT(X)  a  TEMP 
ENDIF 

78  CONTINUE 
80  CONTINUE 


EVENT(O)  a  O.ODO 

Mwsaaaaaasaaaaaa* 

*  Tba  foilowlttf  loepa  dataraiaa  tha  iataraadiata  atataa  of  Tisibil-* 

*  itr  at  each  avant.  a 


DO  too  COUNT  a  0, TOTAL 

IF  (TOTl.EO.Ol  THEN 


TT(C0UNT,1)  a  HSKISI 
GOTO  100 

ELSE  IF  (FVLNT{C0LST).EQ.E1(IS+1))  THES 
TTtCOUST.l)  =  RS1(1S+1) 

GOTO  100 

ELSE  IF  {EVENT(C0UST).GT.E1(IS).AND.TS.EQ.{T0T1-1))  THES 
TTCCOl'KT.l)  a  RSKTOTl) 

GOTO  100 
ESDIF 
90  COSTISIE 
100  COSTISUE 


DO  120  COUNT  a  O.TOT.U 

IF  (TOT2.EQ.O)  THES 
TT(COUNT,2)  a  RS2(0) 

GOTO  120 
EMOIF 

DO  110  IN  a  0,TOT2  -  1 

IF  (EVENT(COU!iT).GE.E2(IS).AND.EVENT(COLST).LT.E2(IN>in  THEN 
n(COUNT,2)  a  BS2(IN) 

GOTO  120 

ELSE  IF  (EYENT(COLNT).EQ.E2(IN+l))  THEN 
n(COUNT,2)  a  HS2(IN4l) 

GOTO  120 

ELSE  IF  (EVENT(C0LNT).GT.E2{IN)..4XD.IS.EQ.(T0T2-1))  THEN 
n( COUNT, 2)  a  RS2(T0T2J 
GOTO  120 

ESDIF 
110  COSTISUE 
120  CONTINUE 

DO  no  COUNT  a  0, TOTAL 

IF  (TOT3.EQ,0)  THEN 
TT{COUNT,3I  a  BS3(0J 
GOTO  140 
ESDIF 

DO  130  IN  a  0,TOT3  -  1 

IF  (EVENT{COUNT).GE.E3(IS)..4SD.EYENT(COUST).LT.E3(IN+l))  THEN 
niCOUNT,3)  a  RS31IN) 

GOTO  110 

ELSE  IF  (EVENK COUNT ).EQ.E3(IS+1))  THEN 
TT(COUNT,3)  =  RS3(IS+1) 

GOTO  140 

ELSE  IF  (EVEltT(COUNT).GT.E3(lN).AND.IS.EQ.(TOT3-in  THES 
TKca'st.a;  *  Bssdors) 

GOTO  140 

ESDIF 
130  CONTINUE 
140  CONTINUE 

DO  160  COUNT  «  O.TOTAL 

IF  (TOT4.EO.O)  THEN 
TT(COUNT,4)  a  BS4(0) 

GOTO  160 
IHDIF 

DO  150  IN  •  0,TOT4  -  1 

ir  (IVIHT(OOinrr).GI.I4(lH).AND.EV'ENT(COUST).LT.E4(ISn))  THIN 
TT(00U1IT,4)  ■  BS4(IK) 

0010  160 

ELSI  IF  (ETnT(C0CNT).n.I4n54l))  THEN 
TT(C0Uirr.4)  ■  tS4(IS4l) 


160  CONTINUE 

DO  180  COUNT  =  0, TOTAL 

IF  <TOT5.EQ.0)  THEN 
TT{COUNT,5)  =  RS5(0) 

GOTO  180 
ENDIF 

DO  170  IS  s  0,TOT5  -  1 

I F  ( EVENT  I  COUNT ) .  GE .  E5  ( I S ) .  AND.  EVENTl  COUNT ) .  LT  .  E5 '  I S+ 1 )  *,  THEN 
TT(COUNT,5)  s  RSodN) 

GOTO  180 

ELSE  IF  (E\’EST(COnfT).EQ.E5(IN+l))  THEN 
TT(COUNT,5)  =  RS5(IS+1) 

GOTO  ISO 

ELSE  IF  {EVENTIC01NT).GT.E5(IS).ASD.IS.EQ.!T0T?-1))  THEN 
TT(COUNT,5)  =  RSj(TOTo) 

GOTO  18C 

ENPIF 
170  CONTINUE 
180  CONTINUE 

DO  200  uWST  =  0, TOTAL 

IF  (TOT6.EO.O)  THEN 
TT(COUNT,6)  =  RS6(0) 

GOTO  200 
ENDIF 

DO  ISO  IN  s  0,TOT6  -  1 

IF  { EVENT { COUNT ) .GE. EC < IN) .AND ^ EVENTl COUNT ) .LT. ESI IN>1 ) )  THEN 
TT(COUSr,6)  =  RSSIIN) 

GOTO  200 

ELSE  IF  IEYENT(C0UN*T).EQ.E5{IN>1))  THEN 
TT I COUNT, 6)  a  8S6(IN+1) 

GOTO  200 

ELSE  IP  (EVE.NT(C«;NT).GT.E6(IS).AND.IS,EQ.(T0T6-1))  THEN 
TT(COUNT,5)  s  R35(T0T6) 

GOTO  200 

ENDIF 
ISO  CONTINUE 
200  CONTINUE 

CHANGE  *  1 

«««*«»«*«««««««•«**««•««««!»*«*«««**«««***«*«««*«»*««*****«**«****«**« 

*  The  foilowintf  detcrainee  if  cotitinuous  coaaunicatlbn . is  possible.*. 

*  For  4  satellites,  the  stui  sust  be  freater  than  or  e^jual  to  3.  '  * 

efeeeeeeeeeeeeseeeseeeseseeeesseeeeeeeseeeesesefeeeeseeeeeees******** 

DO  210  COUNT  •  0, TOTAL 

SUM  •  TT(C0UNT,1)  ♦  TT(COU5T,2)  ♦  TTIC0UNT,3)  ♦  n(C0UNT,4)  ♦ 

♦  TT(COlirr,51  ♦  TT(C0UKT,6) 

IF  (SUN.LT. (3.000). A.ND.CliANGE.EQ.1)  THE.N 
WRITE! 25 ,220 )E?EHT( COUNT ) 

C11A.NGE  «  0 

ENDIF 

IF  (SUN.aE. (3.000). AND.CHANGE.EQ.O)  THEN 
WRITE!  25 , 230)EVE.TT(C0UNT ) 

CHAMR  si. 

ENDIF 

210  CONTINUE 

215  CONTINUE  . _ _  _ ' 


CLOSE ( LTCI  T=STATE+5 ,  STATUS*  *  DELETE  ’ ) 
CLOSE(LNIT=60) 

CLOSE(UKIT=10,STATIS=’DELETE’ ) 

CLOSE ( UN I T= 1 1 , STATUS* ' DELETE ’ ) 

CLOSE ( US I T= 1 2 , STATUS* ’ DELETE  * ) 

CLOSE ( US I T= 1 3 , STATUS* ’ DELETE ‘ ) 

CLOSE (US IT* 14, STATUS* ’DELETE’ ) 

CLOSE ( US I T* 1 5 , STATUS* ’ DELETE ’ ) 


220  FORMAT! ’Couunicatiom  broken  at  ’,F12.5,*  SEC) 

230  FORMAT! ’Coaauni  cat  ions  re-established  at  ',F12.3,'  SEC) 

RETLHS 

END 


) 


SLB80UTISE  GSSCVC4  * 

(Ground  Station  Satellite  Continuous  Visibility.  Check  (i))  * 

n'RPOSE:  The  purpose  of  this  subroutine  is  to  accomplish  a  ,  * 

check  for  constant  visibility  between  a  ground  station  * 
and  3  satellites.  * 


VARIABLES: 

ISPIT: 

STATE 


-  File  number  specifier 


L(X:aL  VARIABLES:  * 

CHANGE  -  specifies  when  communcation  is  broken  or  re-  * 

established  * 

COUNT  -  do  loop  counter 

Elt..|E4  -  Specifies  rise/set  time  for  visibility  function 

a  ground  station  and  4  satellites  * 

EVENT  -  Array  of  all  rise/set  times  * 

IS  -  do  loop  counter  ♦ 

RSI,., RSI  -  1  if  the  event  is  a  rise;  0  if  the  event  is  a  set* 

SUM  -  Sun  of  RS1,..,RS4  at  a  given  event  ,* 

T  -  tine  zero  of  any  selected  visibility  function  * 

-  in  this  subroutine  T  is  a  place  holder  * 

TD!P  -  Temporary  value  used  in  sort  loop  * 

TOTAL  -  Total  number  of  events  » 

TOTl  -  Number  of  events  between  Ground  Station  #?  and  * 

Satellite  *1  * 

T0T2  -  Number  of  events  between  Ground  Station  #?  and  * 
Satellite  *2  * 

TCT3  -  Number  of  events  between  Ground  Station  5?  and  * 

*  Satellite  s3  * 

*  T0T4  -  Number  of  events  between  Ground  Station  =?  and  * 

*  Satellite  *4  * 

*  TT  -  Holds  all  event  times. and  the  corresonding  * 

*  status  of  the  visibility  function  for  each  * 

*  '  \ iewing  path  * 

*  V  -  the  first  value  of  the  selected  visibility  * 

*  function  * 

*  X  -  do  loop  counter  * 

»«»**«»«**s«*««t«***««4t*aS«*:i[«*«$«*«#t*e**4E*4r:^i(c*4c«»***»***4c!|:*4i**t4:*** 

SUBROUTINE  GSSCVCKSTATE) 

RE.AL*8  El  ( 0: 100  > ,  E2  ( 0: 100 )  ,E3(  0: 100 ) ,  E4  ( 0 :  iqO ) ,  EVENT  ( 0:500 ) 

RLAL*8  RSI  1 0 : 100 ) , RS2 ( 0 : 100 ) , RS3 ( 0 : 100 ) , RS4 ( 0 : 100 ) , T , TEMP 

REAL*8  TT(0:500,4),V 

INTEGER  CHANGE, CHECK, COUNT, IS, ST ATE, SUM 

INTEGER  TOTAL, TOTl, T0T2, TOTS, TOT4,X 

0PEN(UNIT=41,FILE=’STAT1.DAT’ ,STATUS=’OLD’ ) 
0PEN(USIT=42,F1LE=’STAT2.DAT' ,STATUS=’OLD’ ) 

OPEN ( UNI T= 4 3 , FI LE= ’ STAT3 . DAT * , STATUS* ’ OLD ’ ) 

OPEN { UNI T=4 4 , FI LE= ’ STAT4 . DAT* , STATUS* ’ OLD * ) 

OPES { UNI Ta60 , FI LE* ’ TRACKER . DAT STATUS* ’ OLD ’ ) 


REWIXD(UNITs60) 

OPEX(USIT*10,FILE»’VIS1.DAT» .STATUSs’CLD’ ) 
0PEX(USIT*ll,FILEs’VIS2.DATSSTATUSs’0LD’) 
0PES(UNITal2,FILE**VIS3.DAT’ .STATUSs’OLD* ) 
OPES(niITal3,FILZ*’VlS4.DAT* .SIAIUSs’OLD’ ) 

sssmsmsmssmssssmmmsf****************************’********************* 


IHJ  iU  tULAl  =.i,lUU 

RE.4D ( STATE ,  ♦ ,  EN'D=25 )  El  I COUKT ) , R31  ( COUNT ) 

TOT 1  =  COUNT 
20  CONTINUE 
25  EKO)  =  O.ODO 

IF  (TOTl.EQ.O)  THEN 
READ(10,*)T,V 
IF  (V.GT. (O.ODO))  THEN 
RSKO)  =  l.ODO 
ELSE 

RSKO)  =  O.ODO 
ENTIF 

ELSE  IF  (RSl(l).EQ. (O.ODO))  THEN 
RSKO)  =  l.ODO 
ELSE 

RSKO)  =  O.ODO 
ENDIF 

tit******t*t**t*t***********tt*****************t**t***t*t***t********* 

*  GROUND  STATION  .AND  SATELLITE  #2  * 

*««4i**4E*»4:*4:*#:*«**4:«*«*»««*:»*««*»«*»**»»4;«»**»****«««»«»««*4(*tt«»«*** 

T0T2  =  0 

DO  30  COUNT  s  1,100 

READ(  STATE+1 ,  * ,  END=35 )  E2  ( COUNT )  ,RS2  ( COUNT ) 

T0T2  =  COUNT 
30  CONTINUE 
33  E2(0)  s  O.ODO 

IF  (TOT2.EQ.O)  THEN 
EE.AD(11,*)T,V 
IF  (V.GT. (O.ODO))  THEN 
RS2(0)  =  l.ODO 
ELSE 

RS2{0)  =  O.ODO 
ENDIF 

ELSE  IF  (RS2(1).EQ. (O.ODO))  THEN 
RS2(0)  =  l.ODO 
ELSE 

RS2(0)  =  O.ODO 
ENDIF 

t***********t**************************************************±***** 

*  GROUND  STATION  .AND  SATELLITE  =3  * 

T0T3  =0 

DO  40  COUNT  =  1,100 

RE.AD  ( ST.ATE+2 ,  * ,  ENT-4  5 )  E3  ( COUNT ) ,  RS3  ( COUNT ) 

T0T3  =  COUNT 
40  CONTINUE 
45  E3(0)  =  O.ODO 

IF  (TOT3.EQ.O)  THEN 
READ(12,*)T,V 
IF. (V.GT. (O.ODO))  THEN 
RS3(0)  ^  l.ODO 

FT  'nF 

RS3(0)  s  O.ODO 
LNDIF 

ELSE  IF  (RS3(1).EQ. (O.ODO))  THEN 
RS3(0)  3  l.ODO 
ELSE 

BS3(0)  s  O.ODO 
ENDIF 

«««««*««*«««*4i«««*«*«*««««**«««*»*«««*««*««»«**««*«**«***«*******«*** 


DO  50  COUNT  =  1,100 

8LAD{ STATE+3 , * , END=53 ) El ( COUNT ) , RSI  ( COUNT ) 

TOT4  =  COUNT 
50  CONTINUE 
53  E4(0)  =  O.ODO 

IF  (TOT4.EQ.O)  THEN 
RL4D(13,»iT,V 
IF  (V.GT. (O.ODO))  THEN 
RS4(0)  =  l.ODO 
ELSE 

RSKO)  =  O.ODO 
ENDIF 

ELSE  IF  (RSKD.EQ.  (O.ODO))  THEN 
RS4(0)  =  l.ODO 
ELSE 

RSKO)  =  O.ODO 
ENDIF 

■^tt***********t***********t***t***************t********************** 

*  The  following  calculates  the  total  number  of  events  between  the  * 

*  ground  station  and  four  satellites.  * 

TOTAL  =  TOT1+TOT2+TOT3+TOT4 
CHECK  =  0 

DO  55  COUNT  s  1,500 

READ  ( 60 , ’T,  END=58 )  EVENT  ( COUNT ) 

CHECK  s  COUNT  , 

55  CONTINUE 
58  CONTINUE 

IF  ( CHECK. EQ.O)  THEN 

W8ITE(25,*) ’There  were  no  visibility  opportunities  with  this 
+  configuration.’ 

GOTO  160 
ENDIF 

**Xt****!t*t*****************************************************1t***t 

*  The  following  orders  all  the  events  from  the  smallest  to  the  * 

*  largest.  * 

**»«4!«««t4:*****«’t:»********«****«***4:**«t***********»***4:**¥4:*4:*4:*:|::t:4:t 

DO  65  COUNT  =  1, TOTAL-1 
DO  60  X  =  1, TOTAL-1 

IF  (EVENT(X+1).LT;EVENT(X))  THEN 
TEMP  =  EVENT (X+1) 

EVENT(X+1)  =  EVENT! X) 

EVENT! X)  =  TEMP 

ENDIF 

60  CONTINUE 
65  CONTINUE 

EVENT! 0)=  O.ODO 

*  The  following  determines  the  status  of  communications  for  each  *' 

*  possible  visibility  path  at  each  event.  * 

**4!*«*«***««4!«««*a«**»s««««**«*«**a«***«*******«*a*******«***«»«*«*»* 

DO  80  COUNT  -  0, TOTAL 

IF  (TOTl.EQ.O)  THEN 
TT(COU.NT,l)  X  BSl(O) 

GOTO  80 
■  .  ENDIF 

DO  70  IN  X  O.TOTl  -  1 

IP  (EVENT(COUNT).GE.El(IN).ASD.EVEST!COljNT).LT.El{IS+l))  THEN 
TT(COUNT,l)  X  eSKIK) 

GOTO  80 


UUlU  ou 

ELSE  IF  (EVEXT(COLNT).GT.E1(IX+1).AXD.IN.EQ.(TOT1-1)>  THEN 
TT(COUST,l)  s  RS1{IN+1) 

GOTO  80 

ENDIF 
70  CONTINUE 
80  CONTINUE 

DO  100  COUNT  =  0, TOTAL 

IF  (T0T2.EQ.0)  THEN 
TT(COUNT,2)  =  RS2(0) 

GOTO  100 
ENDIF 

DO  90  IN  =  0,T0T2  -  1 

I F  ( EVE-NT  { COUNT ) .  GE .  E2  ( I N ) .  AND .  EVENT  I  COUNT ) .  LT .  E2  ( I N+ 1 ) )  THEN 
TT(COLNT,2)  =  HS2(IN) 

GOTO  100 

ELSE  IF  (EVENT(COUNT).EQ.E2{IN+l))  THEN 
TT(COUNT,2)  s  RS2(IN+1) 

GOTO  100 

ELSE  IF  (EVENT(COUNT).GT.E2(IN+l).AND.IN.EO.(TOT2-l))  THEN 
TT(COUNT,2)  =  RS2IIN+1) 

GOTO  100 

ENDIF 

90  CONTINUE 
100  CONTINUE 


DO  120  COUNT  =  0, TOTAL 

IF  (TOT3.EO.O)  THEN 
TT{COUNT,3)  =  RS3(0) 

GOTO  120 
ENDIF 

DO  110  IN  =  0,TOT3  -  1 

I F  ( EVENT ( COUNT ) . GE . E3 (I N ) . AND . EVENT ( COUNT ) . LT . E3 ( I N+ 1 ) )  THEN 
TT(COUNT,3)  =  RS3(IN) 

GOTO  120 

ELSE  I F  { EVENT { COUNT ) . EQ . E3 ( I N+ 1 ) )  THEN 
TT(COUNT,3)  =  RS3(IN+1) 

GOTO  120 

ELSE  IF  { EVENT ( COUNT ) . GT . E3 { I N+ 1 ) . AND . I N . EQ . { TOT3- 1 ) )  THEN 
n(COUNT,3)  =  RS3(IN+1) 

GOTO  120 

ENDIF 
110  CONTINUE 
120  CONTINUE 

DO  140  COUNT  =  0, TOTAL 

IF  (TOT4.EQ.O)  THEN 
TT(COUNT,4)  =  RSKOl 
GOTO  140 
ENDIF 

DO  130, IN  s  0,TOT4  -  1 

I F  ( EVENT { COU-Vr ) .  GE .  E4  ( I N ) .  AND.  El^NT  ( COUNT ) .  LT .  E4  ( I N+ 1 ) )  THEN 
TTlCOnfT,4)  =  RS4(1N) 

GOTO  140 

ELSE  IF  { EVENT (COU'NT).EQ.E4( mi J)  THEN 
TT(COUNT,4)  s  RS4(I.N>1) 

GOTO  140 

ELSE  IP  (EVENT(COn{T).CT.E4(mi).AND.IN.EQ.(TOT4-l))  THEN 
TT{COUNT,4)  a  RS4(IS+1) 

GOTO  140 

ENDIF 

130  CONTINUE  ,  ' 

140  CONTINUE 


*  The  next  DO  loop  determines  if  constant  communication  is  * 

*  possible.  If  the  sum  is  equal  to  zero«  communication  is  not  * 

*  possible.  If  the  sum  is  greater  than  or  equal  to  one,  communi-  ♦ 

*  cation  is  possible.  * 


DO  150  COUNT  ,  =  0, TOTAL 

SUM  =  TTICOUNT,:)  +  TT(C0UKT,2)  +  TT(COUST,3J  +  TTJCOUNT,!) 

IF  (SUM.EQ.{0.0D0)..AN-D.CH.AXGE.EQ.1)  THEN 

KRITE(25,*) ’Communications  broken  at  EVENT {COUNT) 

CHANGE  =0 
ENDIF 

IF  (SUM.GE.(1.CD0).AND.CE\NGE.EQ.0)  THEN 

VfRITEI 25,*) ’Communications  re-established  at  ’, EVENT  (COUNT) 
CHANGE  =  1 
ENDIF 


150  CONTINUE 
160  CONTINUE 


CLOSE (UNI T=STATE , STATUS'  ’ DELETE ’ ) 
CLOSE(UNIT=STATE+l ,STATUS=’DELETE’ ) 
CLOSE ( UN I T=STATE+ 2 , STATUS= ’ DELETE ’ ) 
CLOSE ( UNIT=STATE+3 ,STATUS= ’ DELETE’ ) 
CLOSE (UNIT=60) 

CLOSE { UN I T= 10 , STATUSs ’ DELETE ’ ) 
CL0SE(UNIT=11 ,STATUS=’ DELETE’ ) 

CLOSE ( UN I T= 1 2 , STATUS= ’ DELETE ’ ) 

CLOSE ( L N I T= 13 , STATUS=  *  DELETE ’ ) 

RETURN- 

END 


I 


